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ABSTRACT 


A structural synthesis methodology for the minimum mass design of 
three dimensional frame-truss structures under multiple static loading 
conditions and subject to limits on displacements, rotations, stresses, 
local buckling and element cross sectional dimensions is presented. A 
variety of approximation concept options are employed to yield near 
optimum designs after no more than 10 structural analyses. Available 
options include: (A) formulation of the nonlinear mathematical program- 
ming problem in either reciprocal section property (RSP) or cross sec- 
tional dimension (CSD) space; (B) two alternative approximate problem 
structures in each design space; and (C) three distinct assumptions 
about element end-force variations. Fixed element, design element link- 
ing and temporary constraint deletion features are also included. The 
solution of each approximate problem, in either its primal or dual form, 
is obtained using CONMIN, a feasible directions program (n.b. , dual for- 
mulation not available for all options). 

The frame-truss synthesis methodology is implemented in the COM- 
PASS computer program and is used to solve a variety of problems. These 
problems were chosen so that, in addition to exercising the various 
approximation concepts options, the results could be compared with pre- 
viously published work. The types of problems solved include both 
planar and three dimensional frame-truss structures and contain frame 
members having various cross sectional shapes including: (1) a thin 
walled tube; (2) thin walled box sections; (3) an I section; and (4) a 
solid square section. Finally, the collection of numerical examples are 
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used to form guidelines for the solution of future problems. 
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CHAPTER I 


Introduction 


i*i Introduction 

During the past decade optimization via general nonlinear 
mathematical programming techniques has become widely accepted as a 
viable methodology for engineering design* This has been particularly 
true in the structural engineering field (Refs. 1-2). Here, mathemati- 
cal programming methods have been coupled with finite element based 
structural analysis methods to yield a potentially powerful design tool; 
leading to the emergence of a number of rather general structural syn- 
thesis capabilities (Refs. 3-10). The success of many of the methods is 
due, in large part, to computational efficiencies gained through the 
application of various approximation concepts pioneered by Schmit, et 
al. (Refs. 11-12). 

Vhile the basic methodology for structural synthesis is in place 
for a large class of problems, the majority of the reported computa- 
tional experience has focused on truss and membrane type structures. 
The fact that many practical structures are of this type, or can be ade- 
quately approximated as such, has certainly contributed to this situa- 
tion. There is, however, a significant class of problems for which a 
combined bending-membrane element representation must be used to ade- 
quately capture the essential structural behavior (e.g. frame-truss 
structures). The extension of synthesis methodology to the design of 
these types of structures has been slow and has met with only limited 
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success. The principal difficulty encountered in this case has been 
that of choosing an appropriate set of design variables for which accu- 
rate behavior constraint approximations can be constructed while simul- 
taneously maintaining adequate design freedom. 

The objective of the work reported here is twofold. First, modif- 
ications to the current structural synthesis methodology are suggested 
which will enhance its generality and allow for the more efficient solu- 
tion of bending-membrane structural design problems. Secondly, numerous 
example problems, selected to be representative of larger problems of 
more practical interest, are solved to illustrate the effectiveness of 
the structural synthesis technique described and to provide a body of 
computational experience upon which the solution of future problems may 
be based. 


1 .2 Background 


Much of the early work in the area of optimum design of frame 
structures was motivated by civil engineering applications. Since many 
of the structural systems typical of such applications are built using 
standard section members (e.g. wide flange I-beams) it became popular to 
use assumed size-inertia relationships of the form 


I = cZ P 


( 1 - 1 ) 


where I is the cross sectional moment of inertia, c and p are constants 
and Z is some element sizing variable (e.g, cross sectional area, 
characteristic cross sectional dimension) • These relationships 
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represented assumptions governing the geometry of the element cross sec- 
tion during the re-design process and were frequently based on interpo- 
lation and/or extrapolation of tabulated values for standard sections. 
This approach has the advantage of representing a structural element by 
a single variable and consequently leads to optimum design problems hav- 
ing relatively few design variables. The popularity of this technique 
is illustrated by its extensive coverage in the literature (e.g. Refs. 
13-19) . 

While the use of assumed size-inertia relationships has the advan- 
tage of reducing the design problem size, it also greatly restricts the 
amount of design freedom. Although this lack of design freedom is not a 
serious disadvantage in the design of many civil structures, it can be a 
severe limitation for more weight critical design applications, such as 
in the aerospace and automotive industries where the structural elements 
are usually custom fabrications. As a result, a second approach emerged 
in which some or all of the element cross sectional dimensions (CSD's) 
were selected as the structural design variables (Refs. 20-24) . The 
increase in design freedom and the generality of structural elements 
afforded by this technique lead to its application to increasingly com- 
plex problems. However, as in the case of earlier work on truss syn- 
thesis, it again became apparent that the implementation of approxima- 
tion concepts would be required in order that the method be computation- 
ally viable for large structural systems. 

The integration of approximation concepts into the frame design 
methodology does not, in itself, present any conceptual difficulties. 
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However, as in the case of the truss design problem, the implementation 
requires careful attention to the selection of the intermediate design 
variables so that accurate approximate expressions for the structural 
behavior can be generated. It has been demonstrated that high quality 
first order approximations for nodal displacements and element stresses 
can be constructed using compliance variables (i.e. reciprocal truss 
areas and membrane thicknesses) for moderately redundant truss and mem- 
brane structures (Ref. 25). Indeed, for the statically determinate case 
these behavior approximations are exact when formed using the compliance 
variables. For the frame design problem, the element stresses are, in 
general, complex nonlinear functions of both the element reciprocal sec- 
tion properties (RSP's) and the element CSD's. As a result there is no 
particular choice of intermediate design variable which will yield gen- 
erally high quality approximations for element stresses. However, the 
nodal displacements are well approximated in terms of the element RSP's. 
Therefore it might be expected that the compliance variables (RSP's) 
will yield the best overall behavior constraint approximations for many 
synthesis problems. Several innovative approaches to the frame design 
problem have emerged which are based on this concept. 

One of the most successful of these approaches is based on the 
observation that in the case of thin walled beam sections having fixed 
external dimensions and uniform wall thickness the element RSP's are 
nearly linear in the reciprocal of the wall thickness (Ref. 5). As a 
result, high quality approximations of the structural behavior are 
obtained by selecting the design variables to be the reciprocal wall 
thicknesses of the design elements. This approach suffers somewhat in 
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that the design freedom is obviously limited. The adverse effects of 
this limitation are minimized for cases where design element external 
dimensions are fixed by other considerations such as packaging or 
attachment requirements. This technique has been applied quite success- 
fully to the preliminary design of automotive frame structures (Refs. 
26-27). Extension of the method to cases where the external dimensions 
are also included as design variables has been explored for several 
alternative choices of intermediate CSD variables with moderate success 
(Ref. 28). However, for a general mul ti-variable design element, 
behavior approximations based on first order expansions generated 
directly in terms of the element CSD's or their reciprocals may lead to 
slow convergence and require an excessive number of structural analyses. 

An alternative approach to the frame synthesis problem, which has 
received somewhat less attention for design elements of general cross 
sectional geometry, is to perform the structural design directly in 
terms of the element RSP's. The advantage of such a formulation lies in 
its ability to capitalize on the high quality behavior approximation for 
nodal displacements in terms of the element RSP's. However, since the 
RSP's are treated as independent design variables, a fundamental con- 
sideration must be how the actual physical dimensions of the design ele- 
ment cross section are to be recovered from the RSP's. In general, 
explicit relations for the recovery of the element CSD's are not avail- 
able. In Ref. 29 a technique which coupled the use of RSP's as design 
variables and an approximate CSD recovery method was suggested. In this 
case an approximate linear relationship between the element RSP's and 
CSD's was constructed and used in the CSD recovery process. As 
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originally presented this technique was applicable only when the number 
of CSD's equaled the number of RSP's. This restriction was subsequently 
removed thereby making the method more generally applicable (Ref, 30), 
Unfortunately, even with the available high quality behavior approxima- 
tions, the initial numerical experience with the method of Ref. 29 indi- 
cated the need for larger than expected numbers of structural analyses 
and some convergence difficulties. This may have been due to the 
adverse effect that the linear approximation between element CSD's and 
RSP's has on the net behavior approximation. 
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CHAPTER II 


The Structural Synthesis Problem 


2 .1 Introduction 

Structural synthesis is, by its very nature, a complex, iterative 
process. Fundamentally, this process consists of the generation and 
evaluation of a sequence of trial designs. Each successive design 
represents an attempt to improve some measure of structural performance. 
Historically, design modifications were based on the experience and 
insight of the design engineer. Acceptable designs were frequently 
obtained only after a considerable number of trial designs had been 
evaluated. This was particularly true for complex design problems such 
as those encountered in the design of aerospace structures. This situa- 
tion, together with increased interest in generating designs which were, 
in some sense, optimal, subsequently lead to the development of several 
formal design methodologies based on assumptions as to the number and 
types of critical failure modes (e.g. structural index, fully stressed 
design and optimality criterion methods). A more general method based 
on nonlinear inequality constrained mathematical programming was pro- 
posed by Schmit in 1960 (Ref. 31) and forms the basis for this work. 
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2.2 Problem Formulation 


A significant class of structural synthesis problems may be stated 
as follows: seek a minimum mass design such that all structural behavior 
quantities and design variables remain within specified limits. 
Mathematical ly, this can be written in the form of a nonlinear inequal- 
ity constrained mathematical programming problem as 

min M(D) 

D 

s.t. G(D) 1 0 (2-1) 


5 L i 


1 ? 


where the objective function M is the structural mass, D is a vector of 
design variables, G is a vector of constraints on the structural 


behavior (e.g. nodal displacements, element stresses) and and D L are 
the vectors of upper and lower bounds on D. If it is assumed that the 
structural topology, configuration, materials and loading conditions are 
prescribed, then the design variables D represent element sizing vari- 
ables. For frame-truss structures the element sizing variables are typi- 
cally the element cross sectional dimensions (CSD's) and/or element 
reciprocal section properties (RSP's). The mathematical program 
represented by Eq. (2-1), then, can be rewritten for the frame-truss 
synthesis problem as 
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min M(X,Y) 

X,Y 

s.t. G(X, Y) i 0 

H(X, Y) = 0 

X L i X < X U (2-2) 

Y L i Y i Y° 

where X is the vector of RSP's, Y is the vector of CSD's and X^, X^, Y^ 

and Y^ are their corresponding vectors of upper and lower bounds. The 
equality constraints H(X,Y) have been introduced to account for any 
inter-dependence in the set of sizing variables {X,Y}. Since these con- 
straints are, in general, nonlinear, the solution of the mathematical 
program represented by Eq. (2-2) may be computationally burdensome. It 
is therefore useful to rewrite Eq. (2-2) in terms of a vector of 
independent generalized design variables Z as follows: 

min M(Z) 

Z 

s.t. G(Z) 1 0 (2-3) 

X L 1 X(Z) 1 x u 
Y^ i Y(Z) 1 Y° 

The generalized design variables are, in general, some subset of the 
element CSD's and RSP's. This design problem can be solved for Z, with 
the element CSD's and RSP's being subsequently determined via a recovery 
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transformation of the form 


(X,Y) = T(Z) 


(2-4) 


The stractoral synthesis problem represented by Eqs. (2-3) and 
(2-4) is, in general, a complex, implicit nonlinear problem in terms of 
the generalized design variables. As a result, the direct solution of 
Eqs. (2-3) and (2-4) is computationally impractical even for relatively 
small structures. A more tractable approach to the solution is to 
replace this implicit, nonlinear problem with an explicit approximate 
problem of reduced dimensionality having the following form: 


min MZ^ 
Z. 


s.t. g q (Z L ) <0 ; q€Q R (2-5) 

~i ~ _ ~ n 

X i X(Z L ) < x u 
I 1 Y(Z l ) < Y° 

(X, Y) = T(Z) (2-6) 

where, for the general case, M is an explicit approximation of the 

objective function; the g" are explicit approximations of subset Q n of 

q it 

the original constraints G; is a vector of linked generalized design 

cr U “L 

variables; X , X , Y and Y are upper and lower bounds on the RSP's and 
CSD's chosen to insure the validity of the approximations; and T 
represents some approximate recovery transformation. The solution to 
the original problem (Eqs. (2-3) and (2-4)) is obtained via the itera- 
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tive construction and solution of approximate problems baying the form 
of Eqs. (2-5) and (2-6). Hence, tbe solution to the implicit, nonlinear 
design problem is obtained through the solution of a sequence of expli- 
cit approximate problems. The generation and solution of each approxi- 
mate problem consists of the following four phases: 1) structural 
analysis, 2) approximate problem generation, 3) optimization and 4) 
detail design recovery. Each of these four phases is described in detail 
in Chapters III- VI. 
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CHAPTER III 


Structural Analysis 


3_.l Introduction 

Structural analysis is an integral part of the structural syn- 
thesis problem. The solution of the analysis problem yields the struc- 
tural response quantities (e.g. nodal displacements and element forces) 
required for the evaluation of the design constraints. Various tech- 
niques are available for the linear analysis of frame-truss structures. 
One of the most widely used techniques, and the one chosen here, is the 
well known finite element displacement method. This method is particu- 
larly attractive in the structural synthesis context because 1) a 
variety of different structures and loading can be treated in a unified 
manner, 2) the method is relatively efficient and easy to implement and 
3) the method is well suited for subsequent response quantity sensi- 
tivity calculations (as will be shown in Chapter IV), 

While the finite element method is quite general, the class of 
problems considered here is frame-truss structures subject to multiple 
static loading conditions (including discrete nodal loads and loads uni- 
formly distributed along the element) and homogeneous displacement boun- 
dary conditions. The underlying analysis equations are described in 
detail in the next section. 
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3..2 Static Analysis 


The equations governing the response of a linear structural system 
subject to multiple static loading conditions are of the form 


[KHu) k = {PJ k ; k = 1,2,...K 


(3-1) 


where LK] is the system stiffness matrix, {u> k and {PJ k are the vectors 
of unknown displacements and known applied nodal loads (corresponding to 
the k-th loading condition), and E is the total number of loading condi- 
tions. Eqs. (3-1) represent a set of linear simultaneous equations 
which can be generated from the element level stiffness matrices [K^] e 
and load vectors {P^}* using an assembly technique known as the direct 
stiffness method (Ref. 32). The stiffness matrices and work equivalent 
load vectors (for uniformly distributed loading) for the space frame and 
truss elements are given in Appendix A. 


Prior to the actual assembly of the system stiffness matrix and 
load vectors the element level quantities [K^]° and {P^}® must be 
expressed in terms of a common system level or global coordinate system. 
This is accomplished by using the following transformation equations 

[K^ 8 = [T 1 ] T [K i l e [T i ] 


{P iJ® = [T i 3 T {P i }® 

where [K^] 8 and {P^}£ are the element level stiffness matrix and 
vector, in global coordinates, for the i-th structural element, 
orthogonal transformation matrix [T^] has the general form 


(3-2) 

load 

The 
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[Tj] = 


[til 


[t.] 


[ ti ] 




where (dropping the subscript for convenience) 


and 


It] = [R o ][R e ][R p ] 


[R o ] = 


10 0 
0 cosa sina 

[_0 -sina cosa_ 


l V - 


cosO sinO l" 

-sin© cosO 0 

0 0 0 


[ y - 


cosp 0 sinp 

0 10 
[_-sinp 0 cosp_ 


(3-3) 


(3-4) 


(3-5) 


The angles a, 0 and p, between the local and global coordinate systems, 
are shown in Fig. 1. It should be noted that the matrix [t] for the 
space truss element reduces to the form 


It] = [R e ][R p l 

by virtue of the fact that a may be arbitrarily set to zero making 
an identity matrix. 


(3-6) 

IR a ] 
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Once Eqs. (3-1) have been assembled the homogeneous displacement 
boundary conditions may be applied. Conceptually this is done by elim- 
inating those equations associated with the boundary degrees of freedom 
(in actuality these equations are never assembled). With the appropri- 
ate boundary conditions imposed Eqs. (3-1) represent a positive definite 
system of equations which can be solved for the unknown displacement 
vectors {u}^. The solution method used here is based on a modified 
Cholesky decomposition technique which replaces [K] by a factorization 
of the form 


[K] = [L][D][L] T 


(3-7) 


where [L] is a lower triangular matrix and [D] is a nonsingular diagonal 
matrix. Once [K] has been factorized the solution vectors {u}^ are 
obtained through the usual series of forward and backward substitutions. 
It is important to recognize that significant computational and computer 
storage savings can be realized by taking advantage of the banded struc- 
ture of Eqs. (3-1) . Therefore, in this study, the solution method 
described above is implemented for a compact "skyline" storage arrange- 
ment of [K] as described in Ref. 33, 


Having calculated the nodal displacement vector {u}^; 

k = 1,2, ...E, the end forces for the i-th structural element are given 
by 

{F }® = [K ] ® {u }f + {FEF }® 

(3-8) 

where {F^} a , (u^} a an< * are forces, displacements and fixed 

end forces (corresponding to the uniformly distributed loading) 
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associated with the i-th element for the k-th loading condition, written 
in the local coordinate system. The local displacements are cal- 
culated from the global displacement vector Cn^}^ via the transformation 


<»i>I- IT i ]{ Vk 


(3-9) 


where it is understood that is the subset of the global displace- 
ment vector {u}^ associated with the i-th element. The fixed end forces 
(FEF^}£ for the uniformly distributed loading are given by 


{FEF i )J 



(3-10) 


where is the work equivalent loading vector as defined in Eqs. 

(A-ll) and (A-19) for the frame and truss elements, respectively. 


3_.3. Multiple Boundary Conditions 

The consideration of multiple sets of boundary conditions during 
the analysis of a structural system is quite common in engineering 
design. These boundary conditions may represent actual physical res- 
traints corresponding to varying service environments or they may 
represent "artif icial" boundary conditions created as part of the model- 
ling process (e.g. using a half model for a symmetric structure). The 
application of multiple sets of boundary conditions can significantly 
increase the analysis solution time if the complete system stiffness 
matrix is decomposed for each boundary condition set. This computa- 
tional burden can be greatly reduced by recognizing that, for many 
structures, a significant portion of the system stiffness matrix is 
unaffected by the changes in the boundary conditions (Ref. 5). 
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Conceptually, the system stiffness matrix may be partitioned into por- 

' *1^ the 

boundary condition changes as shovn in Fig. 2. The free portion, IT , 
need be decomposed only once for all boundary condition sets. The 
decomposition of the entire matrix is then completed separately for each 
boundary condition set. This represents a considerable computational 
sayings for structures in which the number of degrees of freedom associ- 
ated with the changing boundary conditions is small compared to the 
total number of system degrees of freedom. It should be noted that the 
degrees of freedom associated with the changing boundary conditions do 
not actually have to be positioned together in the lower portion of the 
matrix as depicted in Fig. 2. This is important since it eliminates 
the need to perform row and column interchanges on the matrix and 
preserves the original matrix bandwidth. 


tions which are independent of (K^p) and dependent on (Kgg 
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CHAPTER IV 


Approximate Problem Generation 


4.1. Introdnction 

The key to a tractable structural synthesis formulation lies in 
the replacement of the original implicit nonlinear design problem with a 
sequence of explicit approximate problems of reduced dimensionality. 
The generation of these approximate problems is accomplished through the 
application of a variety of techniques commonly referred to as approxi- 
mation concepts (Refs. 11-12). Primarily, these techniques serve to 1) 
reduce the numbers of design variables and constraints in the design 
problem and 2) reduce the required number of detailed (exact) constraint 
and objective function evaluations. There are various methods available 
for this purpose. Those implemented here include design variable link- 
ing, temporary constraint deletion and explicit first order constraint 
approximations. These techniques form the foundation of the approximate 
problem generation procedure which consists of the following steps: 1) 
design variable selection, 2) design variable linking, 3) constraint 
evaluation, 4) constraint deletion and 5) objective function and con- 
straint approximation. This procedure is described in detail in the 
following sections. 
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£.2 Design Variable Selection 


In Chapter II it was shown that the frame-truss structural syn- 
thesis problem could be stated, in a general manner, in terms of an 
independent set of generalized design variables Z (Eq. (2-3)). Concep- 
tually, any independent combination of element cross sectional dimen- 
sions (CSD's) and element reciprocal section properties (RSP's) may be 
selected as the design variables as long as the changes in the dependent 
variables can be determined from the changes in the independent vari- 
ables. Two such design variable selection schemes are implemented in 
this study. 

Probably the most natural approach to design variable selection is 
to simply choose the element CSD's as the design variables (CSD design 
space). This has been popular in much of the reported literature (e.g. 
Refs. 20-24 and 26-28) primarily due to the fact that changes in the 
element RSP's (AX) are easily related to given changes in the element 
CSD's (AY) for any cross section shape. Although it is possible to com- 
pute these changes exactly it is useful (in the construction of the 
explicit objective function and constraint approximations) to use the 
following approximate linear relationship to reflect AX in terms of AY 
for each element: 

m - «»> - [8 ] {<« - <v} (4 . 1) 

where can bo determined either through differentiation of analytic 
expressions for the RSP's in terms of the CSD's or from finite differ- 
ence calculations. 
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The primary disadvantage in choosing the element CSD's as the 
design variables is that high quality behavior constraint approximations 
(in terms of the element CSD's) can not be constructed for a variety of 
cross sections. This difficulty subsequently lead to a second approach 
in which the element RSP’s are choosen as the design variables (Ref. 
29) . While this technique offers the potential for the construction of 
high quality behavior constraint approximations (particularily for dis- 
placement constraints), the following inherent difficulties must be 
addressed: 1) changes in the element CSD's are generally not easily 
determined from changes in the element RSP's, 2) the element RSP's may 
not represent an independent set of variables (e.g. the number of RSP's 
is greater than the number of CSD's for a particular cross section) and 
3) the element RSP's may not adequately represent the design freedom 

associated with the CSD's (e.g. the number of RSP's is less than the 

.* 

number of CSD's for a particular cross section). 

The difficulty associated with determining AY in terms of AX is 
due to the nonlinearity of the relationship between the CSD's and RSP's 
which: 1) limits the range of AX changes where linearized approximations 
are useful and 2) admits the possibility that multiple sets of CSD's can 
be found which will yield the same values for the RSP's. As a result, 
an exact representation of AY in terms of AX, for the general case, 
would be difficult if not impossible to determine. 

To overcome the foregoing problem approximate linear relationships 
are constructed between AX and AY for each design element. These rela- 
tionships are obtained by first rewriting Eq, (4-1) as 
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{AX} = [||] {AY} = [ J} {AY} 


(4-2) 


where 

{AX} = (X) - {X Q } 


{AY} = {Y} - {Y q } 


(4-3) 


The desired approximate relationship may now be obtained if Eq. (4-2) 
can be solved for {AY}. For the special case where [J] is square (i.e. 
the numbers of CSD's and RSP's are equal) and non-singular (i.e. the 
RSP's are linearly independent) Eq. (4-3) can be solved directly to 
yield 


{AY} = [J) 1 {AX} 


(4-4) 


and the element RSP's are, indeed, selected to be the generalized design 
variables {Z} . However, this is clearly not the general case and, 
therefore, some alternative method must be employed. Such a method has 
been suggested in Ref. 30 and is described in detail in Appendix 6. 
This technique automatically selects a set of linearly independent 
(free) variables 


m = [{Xj,} {y f j) t 


(4-5) 


and constructs a linear transformation between these variables and the 

T 

remaining dependent (basic) variables [ {X^} {Yg}] of the form 


[ {AXg} {AY b )] T = [HHAZ} 


(4-6) 


where 
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[H] = 


ax 


B 


az 

!Is 

az 


(4-7) 


is calculated from [JJ as shown in Appendix B. An important featnre of 
this technique is not only that it allows for the changes in the depen- 
dent variables to be written in terms of the changes in the independent 
variables, but that the number of independent design variables selected 
is always equal to the number of element CSD's. Therefore all three of 
the difficulties associated with the RSP design space, as described pre- 
viously, are overcome simultaneously. It should also be noted that the 
design variables are selected such that the independent element RSP's 
are chosen first and any additional variables that may be required (to 
make the total number of design variables equal to the number of element 
CSD's) are chosen from the element CSD's. Hence, this design variable 
selection scheme is referred to as the RSP design space. 


Expressions relating the changes in both the element CSD's and 
RSP's to changes in the generalized design variables can now be con- 
structed using Eqs. (4-2) and (4-6). In the CSD design space the fol- 
lowing relationship can be written 

{“} - [Mi nr <“* - 11 1 i|T (az » - iwaz) (4 _ g) 


Similarity, in the RSP design space we may write 


22 



(4-9) 


I ““| * [iF «* ' T* 7?] (4Z) * lH l » T ,4Z) - lDI I4Z > 

These relationships will prove to be useful in the subsequent construc- 
tion of the objective function and constraint approximations. It should 
be noted that frequent updating of these relationships can be obtained 
at low computational cost because they do not involve finite element 
analysis or behavior sensitivity analysis. This will prove to be an 
important observation as will be shown in Section 4.7. 

£.3. Design Variable Linking 

Once the generalized design variables Z have been selected, as 
described in the previous section, design variable linking concepts may 
be employed, thereby reducing the dimensionality of the synthesis prob- 
lem. Topically, design variable linking is used to reflect actual 
design requirements and/or to reduce the problem size enough to make its 
solution tractable. In this latter case the designer is forced to 
approximate the actual design problem in terms of a reduced number of 
design variables in much the same manner that the structural analyst 
must approximate an analysis problem with a limited number of degrees of 
freedom in order that the analysis problem can be solved. Various forms 
of linking are conceivable including the fixing of the relative sizes of 
a given set of design elements of the same type (total linking) and the 
linking of a single design variable of one element to that of another 
element (partial linking). In this work, only total linking between 
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elements of the same type, such that they are identical, has been imple- 


mented. 


±.± Behavior Constraint Evaluation 


The definitions of acceptable structural behavior are central to 
the structural synthesis problem statement. These definitions are 
included in the mathematical problem statement in the form of behavior 
constraints. Two basic types of behavior constraints are included here: 
1) constraints on overall structural stiffness (in the form of nodal 
displacement/rotation constraints) and 2) constraints on local element 
strength (e.g. stress and local buckling constraints). These con- 
straints may be written in terms of the structural response quantities 
(nodal displacements (u) and element forces (F)) and the element ESP's 
(X) and CSD's (I) as follows: 


g q = R q (u(X) ) -110 ; q = 1,2,...^ 

for the displacement constraints, and 


(4-10) 


g = R (F(u(X)),X,Y) - 1 < 0 ; a » 1,2, ...CL 
q q ^ (4-11) 

for the strength constraints, where and are the numbers of dis- 
placement and strength constraints, respectively, and where the response 
ratio represents the ratio of the behavior value to the associated 
allowable and approaches unity as the behavior constraint becomes criti- 
cal. Evaluation of Eqs. (4-10) is clearly straightforward, given the 
values of a particular nodal displacement and its allowable. However, 
the evaluation of the strength constraints (Eqs. (4-11)), in general. 
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requires the evaluation of a different expression for each type of 
design element. The strength constraint formulations for the various 
design element types are given in Appendix C. 

£.5 Constraint Deletion 

Proper design of a structural system usually requires the con- 
sideration of a substantial number of possible failure modes since, in 
general, the critical failure modes are not known at the outset of the 
design process. As a result, the structural synthesis problem statement 
may contain a large number of inequality constraints. In order to 
reduce the number of constraints, and the associated computational bur- 
den. it is possible to temporarily ignore certain constraints which are 
not expected to currently participate in the design. In effect, this 
process reduces the number of constraints by approximating the poten- 
tially critical constraint set. 

The criteria by which particular constraints are judged to be par- 
ticipating (active) or non-participating (passive) forms the basis of 
the constraint deletion technique. Various criteria are conceivable, 
however a relatively simple but effective strategy consists of deleting 
all constraints with response ratios (R ) less than a specified con- 
straint truncation parameter CTP. The value of CTP may, in general, be 
chosen separately for each constraint type and may change during the 
design process. In this work, a single value for CTP is used for all 
behavior constraints. Simply stated, the value of CTP is selected so 
that: 1) constraints with R^ »7 are always retained, 2) constraints 
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with R^ < .3 are always deleted and 3) constraints with .3 £ R < .7 are 

retained or deleted depending on the value of the response ratio cntoff 

value R . This criteria can be written as 
c 


CTP = min 


{ 


max{R ,.3}, 
c 



( 4 - 12 ) 


where R^ is the maximum response ratio rounded down to the nearest tenth 

(e.g. if max R = .65 then R = .6). The value of CTP is updated for 
q€Q q C 

each approximate problem. 


4 ..6 Objective Function and Constraint Approx imat ions 

A key element to the efficient solution of the structural syn- 
thesis problem lies in the construction of accurate explicit function 
approximations. This is particularity true in the case of the behavior 
constraint functions because, in general, exact evaluation of these con- 
straints requires that the structural analysis problem be solved. Vari- 
ous methods are available for the construction of these approximations, 
with the most commonly used techniques requiring only first order 
derivatives of the functions to be approximated (Refs. 12 and 34-35) . 
Two types of first order approximations are used here. 


The first type of approximation consists of expanding the function 
in a linear first order Taylor series of the form 


- „ ~ - _ B 3f(Z ) 

f(Z) = f L (z> = f ( z 0 ) + (Z,, - Z^) 


(4-13) 


where the expansion variables Z are chosen so that the resulting approx- 
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imation is of the highest possible quality. In many cases, however, no 
single set of expansion variables may be chosen such that all function 
approximations are of sufficient quality. In this case it has been sug- 
gested (Ref. 36) that a hybrid or mixed variable approximation might be 
a useful alternative. This approximation can be constructed by the com- 
parison of Eq. (4-13) with a first order Taylor series expansion of the 
form 


B 

f(Z) = fj(Z) = f(Z Q ) + I 


3f(Z 0 ) 
3(1/ Z^ 




(4-14) 


or, equivalently. 


- B 

f(Z) s f (Z) = f(Z ) + I 
• ” 1 « 


3 f ( Z Q ) 
^b 




Subtracting Eq. (4-13) from Eq. (4-15) gives 


fj(Z) - f L (Z) 


B 3 f ( Z_ ) 

I — 

b=l dZ b 


(z b " V 

P 


(4-16) 


For the case where f(Z) represents an objective function to be minimized 
or a constraint function of the form f(Z) _< 0 Eq. (4-16) indicates that 
f j is more conservative than when 


\ 9Z b 


< 0 


(4-17) 


or, if represents some physical variable known to be positive in 
sign, when 
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(4-18) 


3f(Z 0 ) 


< 0 


Consequently, comparison of f^ and on a term by term basis leads to 
the following first order mixed variable approximation: 


_ _ _ B 3f(Z n ) 

f(Z) s f M (Z) = f(Z 0 ) + B b 


where 


3f(Z n ) 

(Z b - z 0 b > if “izf - > 0 


B b- 


-A( 


) 


3f(Z 0 ) 

d/z b ) - d/z 0 ) j if < 0 


(4-19) 


This mixed variable approximation (f u (Z)) is more conservative than 
either the pure linear approximation (f^(Z), see Eq. (4-13)) or the pure 
inverse approximation (fj(Z), see Eq. (4-14) or (4“15)). In this work 
two types of approximations (Eq. (4-13) and Eq. (4-19)) form the basis 
for the objective function and constraint approximations described in 
the following subsections. 


Objective Function Approximations 

The objective function (structural mass) can be written explicitly 
in terms of the element RSP's as 
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(4-20) 


M(X) = 


I 

I 

i=l 



where p^, L i and X^ are the mass density, length, and reciprocal cross 

sectional area, respectively, for the i-th design element. Clearly the 
objective function is easy to evaluate (exactly) when the design is car- 
ried out in the RSP design space and, therefore, no approximation is 
required. However, in the CSD design space Eq. (4-20) cannot be 
evaluated directly to yield the exact value for M. It would first be 
necessary to calculate the design element areas from the element CSD's. 
While this computation is certainly not as burdensome as the detailed 
evaluation of the behavior constraints, it is, never the less, useful to 
replace Eq. (4-20) with an explicit approximation in this case. The 
linear and mixed variable (hybrid) approximations for the objective 
function can be written in terms of the I HSP's as 


M(X) = Mj^X) = M(Xq) 


I 

+ I 
i=l 


3M(X 0 ) 

ax. 

l 



and 


M(X) £ M M (X) = M(X Q ) + 


I 

X 

i=l 


3M(X 0 ) 

3X i 


B. 


where 


(4-21) 


(4-22) 
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9M(X ft ) 

lI i - V lf lif- > 0 


B i- 


- x 

A/1 


^(1/Xi) - (l/x 0 


9M(X n ) 

if -7x^ <0 

1 


Using Eq. (4-8), the approximations given by Eqs. (4-21) and (4-22) can 
be rewritten in terms of the B generalized design variables (Z^) as 


M(X) = M^Z) = M(X Q ) 


B dM(X 0 ) 

+ A <z > ' V 


(4-23) 


and 


_ _ _ B 3M(X 0 ) 

M(X) £ M m (Z) = M(X 0 ) + Z - r y— ■ 
M 0 b=l 3Z b 


(4-24) 


where 


3M(X 0 ) I 3M(X 0 ) 9X i (Z 0 ) 

' A 8i i 8z b 


tz b - v i£ 

D 


B b = 


9M(X 0 ) 

flZ b 


> 0 


- z \{ 


(1/Zb) - (1/z 0 b > ] if 


') 


9M(X 0 ) 


< 0 
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4. .<>.2 Behavior Constraint Approx imations 


The linear and mixed variable approximations for the structural 
displacement constraints represented by Eq. (4-10) may be written in 
terms of the I RSP's as follows: 


I 

g(X) S g L (X) = g(u 0 ) + I 

i=l 


9g(» 0 ) 


du 

dX 


* <x * - \ 


> 


(4-25) 


and 


_ _ _ I 3g(0 

g(X) = g M (X) = g(u Q ) + X Z 2 - 

i=l du 



where 


(4-26) 


- V lf 


a<( °o ) ii , 0 


B. = 

i 


- xj ((1/X.) - (1/X 0 )j if 


d 8 (n 0 ) an 
du 3X i 


< 0 


Similarity, approximations for the element strength constraints (e.g. 
stress and local buckling) represented by Eq. (4-11) may be written in 
terms of the element RSP's and CSD's as 
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g(X,Y) g L (X,Y) - g(F 0 .X 0 ,Y 0 ) + ^ 3X F . + ^ (X i X O i ) 


+ y (y - Y ) 

A aij (Y J Y o.> 


(4-27) 


and 


g(X,Y) 


r M «.D = g(F 0 .i 0 .Y 0 > + [|f wi * ^] B i 


where 


B. 


+ I -^C. 

J-l 3Y J J 


Hi - x 0 1 if ^ it + > 0 

1 °i dF 0X i aX i 


- ^ ( a/i,, - a/i„ 4 ,) a ^ ^ < » 


(4-28) 


[I j - *o 1 “ it > 0 


c. 

J 


- T 


ij ( (i/ v - ci/T 0j >) 


it < 0 

"j 


and where it is understood that the summation over J includes only those 
CSD's corresponding to the design element in which the strength con- 
straint is located. Both the displacement and element strength 
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constraint approximations can be written in terms of the generalized 
design variables (using either Eq. (4-8) or (4-9) depending on the 
selection of the design space) as follows: 


B 


8(Z) = g L (Z) = g(X 0 ,Y 0 ) + ^ <V" Z Qb ) 


(4-29) 


where 


B 


g(Z) - g M (Z) = g(X 0 ,Y Q ) + Z g/ B J) 


b=l 


dZ b 


[Z b " 


3Z b 


> 0 


(4-30) 


- 7 2 


(<*'V - (l/Zo b >) « 


and where 


Jig- = 
3Z b 


8Z b 


< 0 




9Y 


3Y j 3Z b 


qCQj 


Construction of the approximations given by Eqs. (4-29) and (4-30) 
clearly requires the computation of the partial derivatives of the 
structural response quantities with respect to the element RSP's (i.e 

). Various methods are available for the computation of these 

ax ax 

quantities (Refs. 37-38). The technique used in this worh is based on 
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the direct implicit differentiation of the equations governing the 
structural response and it is commonly referred to as the pseudo-load 
method. This method is particularily easy to implement and it rela- 
tively efficient (particularily when implemented using a partial inverse 
technique as described in Ref. 12) when only first order sensitivities 
are required. A detailed formulation of the method is contained in 
Appendix D. 


£.6,. 3. Side Constraint Approximations 

Side constraints on the element RSP's and CSD's of the form 


x L 1 X < x c 


x 1 i X i I c (4-31) 

where jP, X^, and are the upper and lower bounds on the element 
RSP's and CSD's, respectively, play two important roles in the solution 
of the structural synthesis problem. Primarily, the side constraints 
represent bounds on the design element sizing variables corresponding to 
physical design requirements (e.g. packaging limitations, manufactura- 
bility) * Secondly, the side constraints can be used to limit design 
changes during the solution of each approximate problem so as to protect 
the accuracy of the objective function and behavior constraints approxi- 
mations. In this latter case the global upper and lower bounds 

— U — T — n — f. — tj — L “U — L 

(X , X , Y , x ) are replaced by the stepwise bounds (X , X , Y , i ) 

and Eqs. (4-31) become 
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(4-32) 


X L <. X < X U 

f 1, i X < 1° 

The stepwise bounds are calcnlated from move limits (d^, d^) supplied by 
the designer as follows: 


X^ = max[Xj, X £ - d^] 

X? = min[X?, X. + djX.] 

= max [Y^ , Y. - d,Y.] 

J J J J (4-33) 

+ Vj' 

where X i and Y^ are the values of the i-th RSP and j-th CSD at the 
beginning of each approximate problem stage. 

Evaluation of these constraints (Eqs. (4-32)) during the approxi- 
mate problem solution requires that they be rewritten in terms of the 
generalized design variables. For the case in which this solution is 
performed in CSD design space we may write 


Y = Z 


X = X(Y) = X(Z) 


and Eq. (4-32) becomes 


(4-34) 
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Under the additional assumption that the primary element sizing vari- 
ables are the element CSD's, the side constraints on the RSP's may be 
ignored yielding (without approximation) the following form for the side 
constraints in CSD design space: 

Y 1 < Z 

(4-36) 

In the RSP design space the side constraints given by Eq. (4-32) 
must be rewritten in terms of the dependent (basic) and independent 
(free) variables as follows: 

% * VV V 1 


^ ^B^^F’^F* — ^B 


— L — 


(4-37) 


S 1 1 X 


F ^ ~F “ F 


Using Eq, (4-9) the following first order approximation of the side con- 


straints can be written: 



(4-38) 



where {Xg Yg Xg YpJ^f contains the valnes of the element RSP's and CSD's 
at the beginning of the design stage. 

£,£. 4 . Selective Constraint Dependen ce 

A significant amount of the computational effort associated with 
the generation of the approximate design problem is expended during the 
calculation of the partial derivatives of the structural response quan- 
tities with respect to the element RSP's (i.e. ), which are 

dX ax 

required for the construction of the behavior constraint approximations 
(Eqs. (4-29) and (4-30)). For some structural synthesis problems signi- 
ficant reductions in the computational effort can be realized by using a 
selective constraint dependence technique (Ref. 5). This technique is 
based on the observation that in many practical design problems a single 

f 

behavior constraint may be strongly dependent on only a relatively few 
design elements. If these design elements can be identified, then the 
partial derivatives of the constraint with respect to the other design 
variables (for elements upon which the constraint is weakly dependent) 
may be ignored. This technique often leads to dramatic reductions in 
the required number of derivative calculations. Conceptually, selective 
constraint dependence may be applied to both the system displacement and 
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the element strength constraints. However, selection of the design ele- 
ments which strongly influence the displacement constraints may be dif- 
ficult, in general, and it is certainly problem dependent. On the other 
hand, the element strength constraints are well suited to the applica- 
tion of three special cases of the selective constraint dependence con- 
cept. 


Examination of the element strength constraint approximations 
given by Eqs. (4-27) and (4-28) clearly shows that the strength con- 
straints for a given design element are coupled to the design variables 
associated with all other design elements only through the element force 

At? 

derivatives ( _) . It is therefore possible to apply the selective con- 

3X 

straint dependence technique via assumptions made as to the expected 
nature of any element force redistribution which may occur during the 
design process. The following three assumptions are considered here: 
1) the element forces are invariant during the current stage in the 
design process, 2) changes in the forces on a given design element are 
primarily dependent on the design variables associated with that element 
and 3) the element forces are strongly dependent on all of the struc- 
tural design variables. These assumptions lead to the following hierar- 
chy of element force sensitivity calculations: 1) no element force 

derivatives are calculated, 2) element force derivatives are calculated 
only with respect to the RSP's associated with that element and 3) ele- 
ment force derivatives are calculated with respect to the RSP's of all 
elements. 
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i.7 Updating the Approximate Problem 

The extent to which it is possible to minimize the number of 
structural analyses and response quantity sensitivity calculations 

~) required during the design process clearly depends on the qual- 
ax 3X 

ity of the approximate design problems. Approximate problem statements 
for which the underlying constraint approximations are of high quality 
are valid over larger changes in the design variables and, as a result, 
fewer such problems are required to obtain the solution to the actual 
synthesis problem. As discussed previously, the generation of the 
response quantity sensitivities in terms of the element RSP's yields 

values of ~ and ^ which are relatively accurate over large changes in 

ax ax 

X. This, in turn, tends to improve the quality of the behavior con- 
straint approximations, particularily for displacement constraints in 
the RSP design space. However, the overall or net quality of the con- 
straint approximations depend not only on the response quantity sensi- 
tivities but also on the approximated relationships between the element 
CSD's/RSP's and the generalized design variables (e.g. see Eqs. (4-29) 
and (4-30)). Frequently these latter quantities are accurate only for 
small changes in X and Y because of the highly nonlinear relationships 
between the element CSD's and RSP's. Hence, the net approximations may 
be accurate over smaller than desired changes in the design variables. 
To a degree, this problem is less severe when the mixed variable con- 
straint approximations are employed. However, additional computational 
savings may be realized via a procedure which periodically updates the 
approximate problem without recourse to structural re-analysis or 
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response quantity sensitivity calculations. 

The motivation for the approximate problem update procedure 

described in this section lies in the desire to utilize, to the fullest 

extent possible, the quality of the response quantity sensitivities; 

since it is these derivatives which are computationally burdensome to 

generate. To this end the partial derivatives of the constraints with 

respect to the element RSP's and CSD's, (^, are assumed to be of 

3X dY 

high quality and are saved during the approximate problem generation for 
use in the subsequent approximate problem update procedure. This pro- 
cedure consists of the following steps: 


1. calculate new objective function derivatives with respect to 
element RSP's ^ 

ax 


2. update the approximate behavior constraint values (to com- 
pensate for the approximate relationship between the element 
CSD's and RSP's) nsing the equation 


S 


»u A > 


I 

I 

i=l 


JLs. 

ax. 


( V 


x A ) 


(4-39) 


where Xg and X^ are the exact and approximate values of the 
RSP's corresponding to the current values of the CSD's. 


3. 


calculate new values for ^ and ^ (see Eq. (4-8) or Eq. 


dZ 


az 


(4-9)). 
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4. form new side constraint approximations (Eq. (4-36) or (4- 
38)) 

5. calculate new objective function and behavior constraint 
derivatives using 

9M = 3M 91 

0Z 9X 9Z 

a g _ a g ax ^ a g ay 

az ax az ai az (4 _ 40) 

6. form new objective function and behavior constraint approxi- 
mations (Eqs. (4-23 , 4-24) and (4-29, 4-30)). 

Using this procedure it is now possible to update and solve the 
approximate problem repeatedly wihtout recourse to structural re- 
analysis or response quantity sensitivity generation. It should be 
noted that, in practice, the number of times which this update procedure 
may be performed depends on the quality of the reponse quantity sensi- 
tivities. Therefore, it is of paramount importance that the structural 
response quantity sensitivities be generated directly in terms of those 
variables which yield response gradients of the highest possible qual- 
ity, irrespective of the final choice of design variables. 
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CHAPTER V 


Optimization 


5. .A Introduction 

The approximate problem generation techniques discussed in Chapter 
IV make it possible to replace the implicit nonlinear frame-truss syn- 
thesis problem (Eq. (2-2)) with a sequence of explicit approximate 
design problems, each having the form 

min M(Z) 

Z 

s.t g (Z) <0 ; q€Q p 

q K 

~ — ~TT 

X < X(Z) < X 


1 Y(Z) < Y U 


( 5-1) 


where M, g^, X and Y are, in general, explicit approximations of the 
structural mass, retained behavior constraints, design element recipro- 
cal section properties (RSP's) and design element cross sectional dimen- 
sions (CSD's) in terms of the generalized design variables (Z) . When 
these approximate problems are constructed using the objective function 
and constraint approximations described in Chapter IV, Eq. (5-1) 
represents an explicit, separable, convex inequality constrained 
mathematical programming problem. As such, Eq. (5-1) can be solved via 
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any number of veil known, nonlinear constrained minimization techniques 
(Ref, 39), Each of these techniques can be classified as either a pri- 
mal or a dual method, depending on whether the solution is carried out 
in terms of the primal variables (Z) or the dual variables (X). Two 
methods, one of each type, have been implemented here and they are 
described in the following sections, 

£•2 A Primal Solution Method 

Numerous primal methods are available for the solution of the 
mathematical programming problem represented by Eq. (3-1), including 
both direct and transformation (e.g. penalty, barrier) methods. The 
method chosen here is based on the feasible directions method of Zouten- 
dijk (Ref, 40-41) with modifications to improve numerical stability and 
efficiently solve initially infeasible problems (Ref, 42); as imple- 
mented in the CONMIN (Ref, 43) optimization program. This technique was 
selected for the following reasons: 1) the method is applicable to the 
rather general class of problems represented by Eq. (5-1) and 2) the 
implementation of the method, in the form of the CONMIN program, is 
reliable and relatively efficient for the class of problems considered 
here. 

The feasible directions method serves as the primary solution 
technique and can be used to solve any of the approximate problem formu- 
lations shown in Figs, 3 and 4, During the CONMIN solution process some 
computational efficiencies are realized by identifying all side con- 
straint approximations as being linear. Similarily, additional computa- 
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tional savings may be gained by identifying the behavior constraint 
approximation as being linear for problem formulations 1-3 and 7-9. 
This is. however, optional and the default case is to treat the behavior 
constraint approximations as being nonlinear, thereby causing the solu- 
tion to be ’^pushed off" somewhat from the constraint surfaces. 

5,.3 A Dual Solution Method 

An alternative procedure for solving the mathematical programming 
problem represented by Eq. (5-1) consists of replacing this primal prob- 
lem by its dual mathematical programming statement and solving the 
resulting problem in terms of the dual variables (X). This may be done 
by first rewriting Eq. (5-1) in the following slightly more general 
form: 


min H(Z) 
Z 


s.t h^(Z) 1 0 ; q^Q 


Z L 1 Z < Z U 

(5-2) 


where h 

q 


q€Q represent the approximated behavior and side constraints 


and where and are the upper and lower bounds on the generalized 

design variables. The dual of Eq. (5-2) may now be written as 


max 

UO 


min 

Z L i z < z D 



(5-3) 
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where 


L(Z,X) = M(Z) + I X h 

q€Q q q 

Alternatively, Eq. (5-3) may be written as 

max l(X) 

X 2 0 

where 


(5-4) 


(5-5) 


/(X) = m i n L(Z»X) 

(5-6) 

is defined as the dual function. This procedure is viable if it can be 
demonstrated that the dual maximization problem represented by Eq. (5-5) 
has a unique solution (saddle point). It is well known that if the pri- 
mal problem (Eq. (5-2)) is a convex program (i.e. M(Z) and 
h q (Z) ; q£Q are convex functions and Z is contained in a convex subset 
of E n ) and has at least one strictly feasible solution (i.e. there 
exists some Z s.t. h^(Z) < 0 ; q€Q) then the dual problem has a 

♦ — ■ ♦ ■ ♦ 

unique saddle point (Z X ) . If this saddle point can be found* Z is 

the solution to the primal problem (Ref. 44). The existence of a saddle 
point can be demonstrated when Eq. (5-2) represents one of the approxi- 
mate design problems described previously* under the assumption that a 
strictly feasible solution exists. Therefore* in principle* the dual 
solution method may be applied to any of the approximate problems shown 
in Figs. 3 and 4. 
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While the question of saddle point existence is certainly crucial 
in determining the applicability of the dual solution method, another 
important consideration concerns the computational efficiency of solving 
the dual problem (Eq. (5-5)). Clearly, for the general case, maximiza- 
tion of the dual function is considerably complicated by the imbedded 
Lagrangian minimization represented by Eq. (5-6). However, if L(Z, X) 
is additively separable then the solvability of Eq. (5-5) is enhanced by 
the fact that the Lagrangian minimization can be performed as a sequence 
of smaller minimization problems (Ref. 44). The attractiveness of the 
dual solution method is further enhanced by the recognition that for the 
types of approximate problems constructed here the minimization of the 
Lagrangian simply consists of solving a sequence of explicit single 
variable minimization problems. This important observation was first 
made in Ref. 45 in the context of a generalized optimality criteria 
method and subsequently coupled with approximation concepts in Refs. 46 
and 47. 


A final consideration in the application of the dual solution 
technique concerns the method by which the dual function maximization is 
to be performed. Gradient methods are particularily attractive since it 
is well known that the first derivatives of the dual function with 
respect to the dual variables are immediately available from the primal 
constraint values, i.e.: 


iilAi _ v 

a \ q 


/N 

q€Q 


(5-7) 


However, if l(X) does not possess continuous first derivatives such gra- 
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dient techniques may exhibit slov or nonconvergent behavior unless spe- 
cial precautions are taken. The dnal function can be shown to be con- 
tinuously differentiable under the following conditions; 1) Z is con- 
tained in a closed and bounded subset of E n (S), 2) M and q€Q are 

continuous on S and 3) L(Z, X) is minimized over S at a unique point 
Z(X) for all X 2. 0. It can be shown that these conditions are satisfied 
for the case in which the design element CSD's are selected to be the 
generalized design variables and the objective function is approximated 
via a mixed variable (hybrid) approximation. Therefore, an explicit 
mixed variable dual problem can be formulated and solved, via an exist- 
ing first order technique, for approximate problem options 10-12 shown 
in Fig. 4. The mixed variable dual formulation given in the following 
was originally presented in Ref. 48. 

The dual problem can be constructed by first writing the approxi- 
mate primal problems (10-12) as 

min M(Z) 

Z 


s.t g q (Z) < 0 ; qeC^ 


(5-8) 


zJ* < z < z u 

where it is recognized that ? q (Z) ; qCQ^ are the approximations of the 

retained behavior constraints and Z^ and Z^ are the stepwise upper and 
lower bounds on the design variables (element CSD's). Introducing the 
mixed variable approximations for M and g q ; q€Q^ (see Eqs. (4-26) and 
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(4-30)) Eq. (5-8) can be rewritten as 


m 

min I bl Z. - Z t Z? + M 

2 «b>0 ^ b m b <0 ^ °b 


1 C bq Z b" 1 + »q i0 8 q ■ 


°bq>0 * °bq<0 ® b 


4 ^* *5 


b H 1|2|* i*B 


(5-9) 


where 


8M 


J*a 

'bq = 32^ 


N - M(Z 0 ) - I «„ Z * z ^ z 

m. >0 b m. <0 b 


g„ - g (Z n ) - I c. Z n + I c. Z A 

1 « 0 c > 0 b ’ °b c t <0 b<1 °b 

bq bq 


The dual problem (Eq* (5-3)) may now be written as 


max / min L(Z,X) 

^ 0 I Z L ( Z i z D 


(5-10) 


where 
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L(Z,X) = I m_Z, 


"b 


>0 


Vb 


- * ^ Z* + M 

^<0 \ D 


I K 


q€QR 


! I c.Z. - 

c bq b 

C bq>0 


I >z? 


bq<0 


^ ~°b + ^ 


Interchanging the order of the double summation in the fourth and 
terms, Eq. (5-11) can be rewritten as 


UZ,X) = I aZ.- I + 

m^O ^ b m^O h °b 


I < I K o h „ \ Z. - 

C bq>° U eQ R 4 4 ‘ 


C bq <0 


Ij« e Vb 4 z > 


+ I kt + m 


b q€Q, 


qq 


R 


Letting 


% - - °b \ 1 V° 
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and substituting into Eq. (5-12) yields 
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Recognizing that the last two terns of Eq. (5-14) are constant and that 
the remaining terns are additively separable, the ninimization of L(Z,X) 
can be performed via B single variable minimizations, i.e.: 


min L(Z»X) 

z 1 i z < z D 



min 




(5-15) 


where 
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(5-16) 


The solution to the b-th single variable minimization, (temporarily 
ignoring the side constraints on Z. ), is given by 
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( 5 - 17 ) 
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Taking the side constraints into consideration, the solution to the b-th 

single variable minimization becomes 
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Note also the following special cases: 
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Equation (5-20) represents a relatively unconstrained maximization 
problem of a differentiable concave function and, as such, may be solved 
using a gradient based maximization algorithm. The method used here is 
the feasible directions method described previously, where the only con- 
straints are the non-negativity constraints on X. Since the solution is 
carried out in terms of the full set of dual variables (X) the dimen- 
sionality, of the optimization problem is n x n where n is the number of 
retained behavior constraints. Therefore, the dual solution method is 
generally more efficient than the primal method if the number of 
retained constraints is less than the number of primal design variables 
(Z). However, it should be noted that specialized solution schemes for 
the dual problem (e.g. Ref. 46 where the dimensionality of the dual 
space is gradually increased but does not exceed the number of truly 
critical behavior constraints) can make the dual solution method more 
efficient than the primal method even when the number of retained con- 
straints exceeds the number of primal design variables. 
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CHAPTER VI 


Detail Design Recovery 


6..1 Introduction 

A fundamental consideration in structural synthesis is that of how 
the actual detail design quantities (sizing variables) are to be deter- 
mined from the structural design variables. When the relationships 
between the design variables and sizing variables are simple and expli- 
cit (as in the case of a truss design element when selecting A or 1/A as 
the design variables, or a frame design element when selecting the cross 
sectional dimensions (CSD's) or their reciprocals as the design vari- 
ables) then the detailed design recovery process is, of course, trivial 
and is rarely mentioned as a distinct part of the structural design 
problem. However, when the relationship between the design variables 
and sizing variables is not explicit, as in the case where the design 
variables are the element reciprocal section properties (RSP's) and the 
sizing variables are the element CSD's, then the detail design recovery 
process can be quite complex and must be treated as a separate phase of 
the structural synthesis methodology. Two basic types of detail design 
recovery techniques are discussed in the following section for use in 
conjunction with the RSP design space option described in Chapter IV. 
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6.2 The Recovery Process 

The detail design recovery process for the frame-truss synthesis 
problem (in which the element sizing variables are the design element 
CSD's) can be viewed as a procedure for calculating the element CSD's 

(Y) from the vector of optimal generalized design variables Z* 
corresponding to the solution of each approximate problem. In general , 
such a procedure seeks the solution to the set of nonlinear equations 


Z(Y) = Z 

subject to the following restrictions on the element CSD's 


( 6 - 1 ) 


I 1 i I i y d 

( 6 - 2 ) 

For the case where the element CSD's are selected as the generalized 
design variables (CSD design space option) the solution to Eq. (6-1) is 
immediately given by 


Y = 


— * 

Z 


(6-3) 

where it is recognized that the restrictions on the element CSD's 
represented by Eq. (6-2) have already been accounted for in the form of 
bounds on Z (see Eq. (4-35)). Clearly, in this case, the recovery pro- 
cess is computationally trivial and may be carried out without recourse 
to iterative or approximate techniques. This is not, however, the case 
when the synthesis procedure is carried out using the RSP design space 
option. In this case the recovery process must attempt to solve Eqs. 
(6-1) (subject to the constraints represented by Eq. (6-2)) directly. 
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It should be recognized that Eqs. (6-1) may not possess a solution 
within the acceptable domain defined by Eq. (6-2) and, therefore, any 
potential solution procedure must be capable of dealing with this possi- 
bility. As a result, the design recovery procedure for the RSP design 
space option will, in general, be approximate. Two such recovery pro- 
cedures are discussed below. 

Possibly the most natural method for recovering the element sizing 
variables from the generalized design variables is to formulate the 
recovery process as an element level optimization problem of the form 

min tj 
n 2 o 

s.t. -1) i Z* - Z i (T) £ r\ ; i = 1,2, ...M 


yJ* i Yj < yJ ; j = 1.2....N 


(6-4) 


where H and N are, respectively, the numbers of generalized design vari- 
ables and cross sectional dimensions associated with a given design ele- 
ment. This procedure seeks the solution to Eqs. (6-1), one design ele- 
ment at a time, so as to minimize the maximum error in any one equation 
while forcing the solution (T) to lie within the allowable upper and 
lower bounds (a similar method is suggested in Ref. 23). This procedure 
will tend to yield the "best" approximate solution to Eqs. (6-1) and, as 
a result, has the advantage of preserving, to the extent possible, the 
quality of the behavior constraint approximations. Unfortunately, the 
recovery scheme summarized by Eqs. (6-4) would require the solution of 
many (equal to the number of design elements for each approximate 
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problem stage) nonlinear mathematical programming snbproblems. This 
would be computationally burdensome and therefore the recovery scheme 
represented by Eqs. (6-4) has not been implemented in this work. 


An alternative recovery procedure, which has been implemented 

here, was first used in Ref. 29. This procedure makes direct use of the 

previously constructed approximate linear relationships between the 

changes in the element CSD's {AY} and the changes in the generalized 

design variables {AZ} (see Eq. (4-9)). This relationship can be written 
1 

as 


(AT) - gfJtAZ) 

where is constructed as described in Chapter IV. Using Eq. 

the actual recovered values for the element CSD's are given by 


(6-5) 

(6-5) 


{Y} = (Y 0 ) + (AY) 


( 6 - 6 ) 


where {Y^} contains the values of the element CSD's at the beginning of 
the design stage. Clearly, this procedure requires few additional com- 
putations and, therefore can be applied efficiently to large problems. 
The main disadvantage of the method lies in the fact that the linear 
approximation (Eq. (6-5)) may be valid for only relatively small changes 
in the generalized design variables. This can require the use of tight 
move limits which, in turn, can make it necessary to construct and solve 
an excessively large number of approximate design problems. Fortunately, 
this difficulty can be effectively overcome by using the approximate 
problem update technique (described in Chapter IV) in which the rela- 
tionships given by Eq. (6-5) are periodically updated during the 
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solution of an approximate design problem while holding the results of 
the structural analysis and behavior constraint sensitivity analysis 
invariant (until the beginning of the next stage) . 


57 



CHAPTER VII 


Program Description 


2.1 Introduction 

The frame-truss synthesis methodology described in Chapters II-VI 
has been implemented in the COMPASS ( Comp uter Program for Analysis and 
Synthesis of Space-frames) computer program. This program is intended 
to serve as a research code for the study and development of practical 
and efficient synthesis techniques for structural systems whose essen- 
tial structural behavior requires a bending-membrane element representa- 
tion. While primarily a research tool, the program is capable of solv- 
ing problems large enough to be of some practical interest. Although 
its primary function is structural design, the program can be used for 
basic structural analysis, with or without design modelling data. A 
command oriented input data structure makes the program relatively easy 
to use. Also, the programs modular organization and in-core storage 
management system serve to facilitate future expansion and development 
efforts. 

2.2 Scope of Program 

The COMPASS program is currently capable of determining the 
minimum mass design of three dimensional frame-truss structures subject 
to multiple static loading conditions. The structural topology, confi- 
guration, material and loading information is supplied in the form of a 
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finite element-analysis model and is assumed to be invariant during the 
design process. The structural behavior (nodal displacements and member 
end forces) is determined via a linear displacement method finite ele- 
ment technique, using a combination of space frame and truss elements 
( see Appendix A) . 

The structural synthesis problem is solved using the sequence of 
approximate problems approach pioneered by Schmit et al. (Refs. 11-12). 
Both a first order Taylor's series approximation and a mixed variable 
approximation (Ref. 36) are available for construction of the explicit 
behavior constraint functions. Each approximate problem is constructed 
and solved in a generalized design variable space, consisting of either 
cross sectional dimensions (CSD's) or a combination of reciprocal sec- 
tion properties (RSP's) and CSD's, vith the ultimate goal being that of 
determining optimum values for the element sizing variables. These siz- 
ing variables are associated with the various design element cross sec- 
tion shapes described in Appendix C. User specified bounds on the ele- 
ment sizing variables prevent the design from assuming unrealistic 
dimensions. Move limits can also be applied to the design variables to 
ensure that the behavior of all candidate designs is well represented by 
the approximate problem. A design element linking capability which 
links all sizing variables between selected design elements is also 
available. 

The COMPASS program allows the user to design against a variety of 
failure modes. Limits on nodal displacements and rotations and on 
design element strength (e.g. stresses and local buckling) can be 
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treated. Differences in the design element cross section shapes require 
that, in general, the strength failure criteria be tailored to the 
specific design element. The failure criteria for the various design 
elements are described in Appendix C. 

The approximate design problem may be solved using either a primal 
or dual mathematical programming algorithm depending on the users choice 
of design space and constraint approximation technique. The CONMIN 
(Ref. 43) computer program, based on a feasible directions method (Ref. 
43), is used to solve the primal form of the approximate design problem. 
The solution to the dual form of the design problem is based on the 
development described in Ref. 48; with the CONMIN program being used to 
perform the actual dual function maximization. The available combina- 
tions of optimization method, design space and constraint approximation 
techniques are illustrated in Fig. 5* 

2.3 Organization 

The basic program organization is shown in Fig. 6. Pre and post 
processing routines are provided to perform one time input and output 
data processing functions. The design control routine directs the exe- 
cution of the following four primary design functions; (1) structural 
analysis, (2) approximate problem generation, (3) optimization and (4) 
detail design recovery. 

The design process begins with the structural analysis phase from 
which exact values for the structural behavior quantities (nodal dis- 
placements and element end forces) are obtained based on the initial 


60 



design input. The approximate problem generation routine then performs 
the following operations; (1) design variable selection and linking (2) 
constraint evaluation, (3) constraint deletion, and (4) objective func- 
tion and retained constraint approximation. The approximate design 
problem is then passed on to the optimization routine where the re- 
design function is performed. Finally, the design element sizing vari- 
ables are calculated from the design variables in the detail design 
recovery phase. The new design is then passed back to the design con- 
trol block where the entire process is repeated until design convergence 
is achieved. 

Under certain circumstances the user may wish to periodically 
bypass the structural analysis phase, proceeding directly to the approx- 
imate problem generator as depicted by the dashed line in Fig. 6. In 
this case the approximate problem is updated without recourse to struc- 
tural analysis or response quantity sensitivity calculations as outlined 
in Chapter IV. The approximate problem updating procedure allows the 
high quality displacement derivatives to be utilized over larger changes 
in the design variables than would otherwise be possible* It is impor- 
tant to note that considerable gains in solution efficiency can be real- 
ized by bypassing the structural analysis and displacement derivative 
calculations in this way. 
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2..£ Storage Management 


Program data storage management is one of the most important con- 
siderations in program development. Unfortunately, in research program- 
ming it is often times ignored. It is generally agreed that some type 
of centralized data base management system (Ref. 49-50) is essential to 
the development of modern large scale engineering analysis/synthesis 
programs. However, it is less clear that the implementation of such 
data base management schemes in a research code is time and/or cost 
effective. The COMPASS program utilizes an in-core storage management 
facility which attempts to provide some of the benefits of general data 
base management within the constraints of the research environment. 

The basic concept-behind the storage management system implemented 
here is shown in Fig. 7. The storage manager consists of three parts; 
(1) a data vector, (2) a dictionary and (3) a void area table. All pro- 
gram data is stored within the data vector while the associated descrip- 
tions, locations and lengths are stored in the dictionary. The void 
area table contains the locations and lengths of unused portions of the 
data vector that may appear as the program data storage is altered dur- 
ing execution. While the storing of program data in a single data vec- 
tor is quite common in engineering programming, the associated descrip- 
tive information is usually maintained outside of the program. Here, 
the programmer is allowed complete control over the access, creation, 
deletion and alteration of data storage from within the program through 
the use of a variety of storage management commands (Fig. 8) . Such a 
capability leads to increased programming flexibility and facilitates 
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program maintenance and development. Finally, although the use of such 
a storage management system does increase solution time, computational 
experience indicates that these increases are small when the system is 
efficiently applied. 

2.5 Implementation 

The COMPASS program is operational as a stand alone program on the 
IBM 3033 computer using the MVS/SP operating system. The program con- 
sists of approximately 14,000 Fortran, 1000 PLI and 50 Assembler state- 
ments. The current version of the program requires approximately 650 K 
bytes of memory, excluding data storage. Program storage requirements 
may be considerably reduced, however, by taking advantage of the 
program's structure through the use of overlay or segmentation tech- 
niques. The detailed flow diagrams shown in Figs. 9-17 are provided as 
a guide for the overlay process. It should be noted, however, that 
memory savings realized through the overlay process may not be signifi- 
cant for problems in which the data storage is large compared to the 
program storage. 

The standard I/O device unit designations (FT05 and FT06 for For- 
tran, SYSIN and SYSPRINT for PLI) are used for program input and output. 
Auxiliary external files are also required for certain program options. 
File numbers 10 and 11 are required for all problems for which multiple 
boundary condition sets are specified. File number 12 is required if 
the approximate problem generation update procedure is enabled. For 
problems requesting program data checkpoint or restart options the user 
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must define the files associated with the checkpoint or restart file 
numbers. 

The implementation of the COMPASS program on other computer sys- 
tems would require some program modifications. For computer systems on 
which the PLI programming language is available the required modifica- 
tions are relatively minor and confined to the storage management sub- 
routines. On systems where the PLI programming language is not avail- 
able the input data subroutines would have to be rewritten in Fortran or 
some other suitable language. The small amount of Assembler code is 
used only for CPU timing and is easily replaced by any equivalent system 
CPU timing routine. 

2.6 Input Data Commands 

The COMPASS program input data format is designed to be easily 
used and highly flexible. A problem oriented free format command 
language is coupled with a data scanning feature to provide a data entry 
method which is essentially free of organizational and formatting res- 
trictions. For convenience, the data input stream is divided into three 
sections; (1) analysis data, (2) design data and (3) control data. Each 
section is headed by a data block command. Commands which describe the 
structure and its loading are supplied in the analysis data block. This 
data is similar to that provided to most finite element structural 
analysis programs. The design data block is used to supply the informa- 
tion associated with the structural design problem (e.g. initial design 
data, side constraints, and behavior allowables). Finally, all of the 
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program control information is supplied via tlie control data block. 

Three input data command forms are shown in Fig. 18. All of the 
program data commands appear in one of these three basic forms (or 
slight variations thereof). In form 1 the command is followed on the 
same line by its associated data. A command followed by several 
separate lines of data is shown in form 2. Finally, form 3 shows a com- 
mand followed by several lines of sub-commands and data. In all cases, 
only the underlined portion of the commands or sub-commands need to be 
given. However, inclusion of the full command phrase is allowed. Some 
commands contain optional data (denoted by a quantity enclosed in brack- 
ets, e.g. [data]) which may or may not be specified at the user's dis- 
cretion as outlined in the command description. In this case a default 
value is assigned for the missing data item. Comment cards are allowed 
in the data stream and are designated by placing a dollar sign (f) in 
column 1 of the data card. Data may also be continued from one card to 
the next by placing a continuation character (-) in column 72 of the 
card which is to be continued. Data entity (element, node, load set, 
etc.) numbering is also unrestricted. For example, node point numbering 
does not have to begin with the number 1 or be in ascending order. The 
available data commands are described in Appendix E. A set of sample 
data input is shown in Fig. 19. 
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2.2 Restrictions and Limitations 


The COMPASS program has relatively few operational restrictions. 
The major restriction is that of problem size. Like most programs which 
make extensive nse of an in-core data structure, the maximum solvable 
problem size is dependent on the amount of computer memory available. 
The use of the storage management system described in Section 7.4 helps 
to alleviate this problem but does not eliminate it. 

The dynamic nature of the synthesis problem makes data storage 
requirements difficult to estimate. In general, however, the problem 
data requirements are most effected by the number of degrees of freedom 
and bandwidth in the analysis model and the number of design variables 
and constraints in the design model. Careful attention to node number- 
ing will help reduce the analysis model bandwidth. To the extent that 
it is possible, design variable linking can help to reduce the dimen- 
sionality of the design space. 
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CHAPTER VIII 


Numerical Examples 


8.1 Introduction 

In this chapter, the detailed results for numerous structural syn- 
thesis problems are presented. Each problem has been solved using the 
previously described frame-truss synthesis methodology as implemented in 
the COMPASS computer program on the IBM 3033 at UCLA. The example prob- 
lems were selected so that, in addition to exercising the various solu- 
tion options outlined in Table 1, the results could be compared with 
previously published work. The types of problems solved include both 
planar and three dimensional frame-truss structures subject to multiple 
static loading conditions with constraints on nodal 
displacements/rotations and on element strength (e.g. stress and local 
buckling). These problems contain frame members having various cross 
sectional shapes including: 1) a thin walled box beam with 4 cross sec- 
tions dimensions (CSD's) (B, H, t^, t^) , 2) an I beam symmetric about 
the x-y plane with 6 CSD's (B^, B^, H, t^, t 2 , t^), 3) a solid square 
beam with 1 CSD (B) , 4) a thin walled tubular beam with 2 CSD's (R,t) 
and 5) a thin walled box beam with 3 CSD's (B, H, t) (see Appendix C, 
Figs. C3, C5-C8). 


67 


8., 2 Tied Cantilevered Beam ( Problem 1) 

Figure 20 depicts a tied cantilevered beam subject to two indepen- 
dent loading conditions. The structure is modelled using one frame ele- 
ment (member 1) and one truss element (member 2). The members are made 
of materials having the same modulus of elasticity and weight density 
but having different yield stress allowables. The truss member (design 
element type 1) is described by its cross sectional area while the frame 
member (design element type 13) has a square cross section with one siz- 
ing variable (B) . This structure is designed for minimum weight subject 
to element stress constraints at the ends of each member and side con- 
straints on the element sizing variables. A summary of the material 
properties, loading conditions, constraint allowables, initial design 
and bounds on the element sizing variables is given in Table 2. The 
design element descriptions are given in Appendix C. 

This problem was solved using four different solution options. 
The iteration history data for these runs is given in Table 3. The 
corresponding iteration history plots are shown in Figs. 21-24. In the 
first three runs (options 1(P), 2(P) and 3(P)) the design is carried out 
in the RSP design space using linear constraint approximations and 40% 
move limits on the ESP's (d^ - 0., ~ .4). Each approximate problem 
is solved via a primal solution method (CONMIN) • The differences 
between runs 1, 2 and 3 lie in the assumptions made regarding the ele- 
ment end force variations during the solution of each approximate prob- 
lem. Comparison of the results for runs 1, 2 and 3 indicates little 
difference in the number of analyses required for convergence as the 
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amount of element end force sensitivity information included in the 
stress constraint approximations is increased. However, the final 
design weight for run 3 is approximately 9% less than that of runs 1 and 
2. This result is not unexpected since this problem is known to possess 
several local minimum solutions (Sef. 51). It is, however, interesting 
to note that this solution was obtained as a result of an improvement in 
the quality of the constraint approximations. 

The fourth run (option 6(P)) for this problem is the same as run 3 
except that the element stress constraints are approximated via a mixed 
variable (hybrid) approximation. There is no improvement in the conver- 
gence rate in this case; indeed the results are identical to those in 
run 3. This is due to the simple form of the stress constraints which 
are well approximated in terms of the design element RSP's. 

The final designs and critical* constraints for runs 1-4 are 
shown, along with those of the reference solution (Ref. 51), in Tables 4 
and 5. The three reference solutions, designated as Method I, Method 
II-B and Method IV-B, were obtained using the same assumptions regarding 
the element end force variations as runs 1, 2 and 3, respectively. The 
comparisons between the final designs for runs 1-3 and the corresponding 
reference solutions is quite good in terms of both final weight and 
material distribution. The largest differences occur for run 2 where 
the reference solution is approximately 3.3 lb. lighter. However, in 
this case the stress constraints for the reference solution are slightly 


* A critical constraint is defined here as any constraint having 
an associated response ratio (ratio of the response value 
to its allowable) greater than or equal to .95. 
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violated (see Ref. 51). Also, it is interesting to observe that the 
solution obtained in runs 3 and 4 is distinctly different than that of 
runs 1 and 2 both in terms of material distribution and critical con- 
straints. Again, this is not unexpected since the structure clearly has 
two competing primary load paths. 

8.3 Two Member Frame ( Problem 2 ) 

A two member plane frame subject to a single out of plane load is 

shown in Fig. 25. This structure is modelled using two thin walled box 

section frame elements (type 11) having four sizing variables (B, H, 

t, , t.) . Both members are made of the same material. The two member 
d n 

frame is designed for minimum mass subject to two independent sets of 
constraints. The first set (Case A) includes stress constraints on both 
members and side constraints on the element sizing variables. The 
second set (Case B) consists of the constraints included in Case A with 
the addition of constraints against local wall budding of the members. 
A summary of the material properties, loading conditions, constraint 
allowables, initial design and bounds on the element sizing variables is 
given in Table 6. The design element is described in Appendix C. 

8..3..1 Case A: Stress and Side Constraints 

This case was solved using eight different solution options. The 
iteration history data for these runs is given in Table 7. The itera- 
tion history plots are shown in Figs. 26-33. The first three runs 
(options 1(P), 2(P) and 3(P)) are made using the RSP design space option 
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and linear constraint approximations. Each approximate problem is 
solved via a primal solution method (CONMIN) with 40% move limits on the 
RSP's (d^ = 0., dj = .4). The differences between runs 1, 2 and 3 lie 
in the assumptions made regarding the element end force variations dur- 
ing the solution of each approximate problem. Comparison of the results 
for runs 1-3 shows little change in terms of the number of analyses 
required for convergence. However, the results do improve in terms of 
the maximum constraint violation for intermediate designs as the amount 
of element end force sensitivity information contained in the approxi- 
mate problem is increased. Indeed, in the case where the element end 
forces are assumed to be dependent on all of the design variables (run 
3) most of the intermediate designs are feasible. Again, this is not 
surprising since the structure clearly has two competing load paths. 

Due to the lack of significant convergence improvement for either 
the local or global element end force variation options, the remaining 
runs were made assuming that the element end forces are invariant during 
the solution of an approximate problem. Hun 4 (option 4(P)) is the same 
as run 1 except a mixed variable (hybrid) approximation is used for the 
stress constraints and the move limits are increased to 50% 
(d 1 = 0., d 2 = .5) on the RSP's. Comparison of the iteration histories 
for runs 1 and 4 shows a slight improvement in convergence rate and a 
significant improvement in terms of maximum constraint violation for 
intermediate designs, especially considering the more liberal move lim- 
its. This improvement can be attributed to the conservativeness of the 
mixed variable approximation as compared to the pure linear approxima- 
tion used in run 1. It is also interesting to note that the iteration 
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history for run 4 compares quite favorably with that of ran 3 indicating 
that, at least for some problems, the mixed variable approximation may 
be used successfully in place of higher levels of element end force sen- 
sitivity information (without the associated computational expense). 

In runs 5 and 6 (options 7(P) and 10(P)) the design is carried out 
in CSD space vith 25% move limits on the CSD's (d^ = .25, dj = 1.0). 
Each approximate problem is solved using a primal solution method. In 
run 5 the stress constraints are approximated using a linear approxima- 
tion while in run 6 the mixed variable approximation is used. Com- 
parison of the iteration histories for these runs shows superior results 
in terms of both convergence rate and maximum constraint violation for 
run 6. However, run 6 does have a slightly higher final mass resulting 
from a different and slightly less efficient material distribution. 
Also, run 6 compares well with run 4 except, again, for the slightly 
higher final mass. 

Run 7 (option 10(D)) is the same as run 6 except that each approx- 
imate problem is solved using a dual solution method. Comparing the 
iteration histories for runs 6 and 7 reveals little significant differ- 
ence. For this problem the dual solution method appears to be more 
efficient even though the numbers of retained constraints and design 
variables are the same. 

The final run (option 1(PU)) is the same as run 1 except that the 
approximate problem update procedure is used. In this case the approxi- 
mate problem is reconstructed, without recourse to structural analysis 
or response quantity sensitivity calculations, once between each 
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complete approximate problem generation. Comparison of the iteration 
histories for rnns 1 and 8 shows an improvement in the convergence rate 
by a factor of two while maintaining comparable maximum constraint vio- 
lations for the intermediate designs. The number of structural analyses 
required for convergence in this case is 7, resulting in the best con- 
vergence rate of all solution options used for this problem. 

The final designs and critical* constraints for runs 1-8 are 
given, along with those of the reference solution (Ref. 29) in Tables 8 
and 9. All of the final designs represent the same (intuitively 
correct) design concept. The sizing variables of the longer member (1) 
are at their lower bounds, with most of the load being carried through 
the shorter member (2) and the critical stresses occurring at the fixed 
end of this member. While these designs are conceptually the same it is 
interesting to note that there are two competing means of carrying the 
load through member 2. In the final designs for runs 1-5 and 8 member 2 
achieves nearly the maximum allowable outer dimensions (B.H) and minimum 
thickness (t^» t^) . However, the final designs for the reference solu- 
tion and runs 6 and 7 have a significantly smaller base dimension (B) 
and a larger wall thickness (t^). While the first design concept is 
more efficient (and intuitively more satisfying), the final design 
masses corresponding to the second design concept are only slightly 
higher (1-3%). This is not too surprising since it has become well 
recognized that many structural design problems do exhibit practical 
local minima having relatively close values of the objective function. 
As in this case, these local minima are often associated with distinct 
design concepts. 
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£.3,. 2 Case B: Stress , Buckling and Side Constraints 

This case was solved using seven different solution options. The 
iteration history data for these runs is given in Table 10. The itera- 
tion history plots are shown in Figs. 34-40. As in Case A, the first 
three runs (options 1(P), 2(P) and 3(P)) are made using the RSP design 
space option and linear constraint approximations. Each approximate 
problem is solved using a primal solution method with 40% move limits on 
the RSP's (d^ =0., d^ = .4). Comparison of the iteration histories for 
runs 1-3 indicates no change in the convergence rate and only slight 
improvement in the amount of constraint violation for intermediate 
designs as the amount of element end force sensitivity information con- 
tained in the approximate problem is increased. 

Due to the lack of significant convergence improvement for either 
the local or global element end force variation options, the remaining 
runs were made assuming that the element end forces are invariant during 
the solution of an approximate problem. Run 4 (option 4(P)) is the same 
as run 1 except that a mixed variable (hybrid) approximation is used for 
the stress and local buckling constraints and the move limits on the 
RSP's are increased to 50% (d^ =0., dj = .5). Comparison of runs 1 and 
4 shows a slight improvement in convergence rate and a significant 
improvement in the maximum constraint violation history, especially con- 
sidering the larger move limits. 

In runs 5 and 6 (options 10(P) and 10(D)) the design is carried 
out using the CSD design space option with mixed variable stress and 
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local buckling constraint approximations and 25% move limits on the 
CSD's (d^ *= .25, ** 0*). The approximate problems are solved using a 
primal and a dual solution method, respectively. In this case, com- 
parison of the iteration histories with those of the designs obtained 
using the RSP space design option (runs 1-4) indicates only slightly 
better convergence rate but a significantly improved maximum constraint 
violation history. Indeed, in both runs 5 and 6 most of the intermedi- 
ate designs are feasible. However, it should be noted that both of 
these runs have final mass values slightly greater than runs 1-4, 
resulting from a somewhat less efficient material distribution (see 
Table 11). 

The final run (option 1(PU)) is the same as run 1 except that the 
approximate problem update procedure is used. As in run 8, Case A, the 
approximate problem is reconstructed once between each complete approxi- 
mate problem generation. Again, the use of this procedure results in a 
significant improvement in the convergence rate, although the maximum 
constraint violation for the intermediate designs are quite large. The 
number of structural analysis required for convergence in this case is 
7, the best of all solutions options used for this problem. 

The final designs and critical constraints for runs 1-7 are 
given, along with those of the reference solution (Ref. 29), in Tables 
11 and 12. Again, as in Case A, all of the final designs represent the 
same overall design concept with the load being carried through the 
shorter member (2). Also, as in Case A, the same two competing means of 
carrying the load through member 2 are present. They are represented by 
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the designs from runs 1-4, 7 and the reference solution and runs 5-6, 
respectively. It is interesting to note, however, that in this case 
these competing design concepts not only have different final masses but 
they also are associated with different sets of critical constraints 
(see Table 12). In runs 1-3 and 7, the large base dimension (B) and 
small base thickness (t^) lead to both critical stress and local wall 
buckling constraints. On the other hand, the narrower and thicker base 
wall dimensions (B, t^) for the final designs of the reference solution 
and runs 4-6 preclude criticality of the wall buckling constraint. 

8.4 Three Member Frame ( Problem 3.) 

Figure 41 depicts a three member planar frame subject to two 
simultaneous out of plane loads. This structure is constructed of three 
thin walled box section frame elements (type 15) each having three siz- 
ing variables (B, H, t). All members are made of the same material. 
The three member frame is designed for minimum material volume subject 
to constraints on the maximum allowable member stresses. A summary of 
the material properties, loading conditions, constraint allowables, ini- 
tial design and bounds on the element sizing variables is given in Table 
13. The design element is described in Appendix C. 

This problem was solved using seven different solution options. 
The iteration history data for these runs is given in Table 14. The 
iteration history plots are shown in Figs. 42-48. The first three runs 
(options 1(P), 2(P) and 3(P)) are made using the RSP design space option 
and linear constraint approximations. Each approximate problem is 
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solved via a primal solution method with 50% move limits on the RSP's 
(d^ = 0., dj “ .5). The differences between runs 1, 2 and 3 lie in the 
assumptions made regarding the element end force variations during the 
solution of each approximate problem. Comparison of the iteration his- 
tories for these runs reveals no difference in the convergence rate for 
the three solution options. However, there is a substantial difference 
between runs 1 and 3, and run 2 in terms of the maximum constraint vio- 
lation for several of the intermediate designs. Two interesting obser- 
vations can be made here. First, the performance of the solution 
options 1(P) and 3(P) are nearly identical even though it would appear 
that, for a statically indeterminate structure such as this, the global 
element end force variation option (3(P)) would yield superior results. 
Secondly, the local element end force variation option (2(P)) yields 
substantially poorer results in terms of maximum constraint violation 
than the invariant option (1(P)), particularly during the early stages 
of the design. This behavior can be attributed to the fact that, 
although the structure is statically indeterminate, the symmetry of the 
structure and loading leads to a synthesis problem in which the coupling 
member (2) tends to vanish, thereby reducing the problem to the design 
of two determinate cantilevered beams (members 1 and 3). It should now 
be recognized that the local element end force variation option may, in 
some cases, lead to constraint approximations which are of poorer qual- 
ity than those resulting from the assumption that the element end forces 
are invariant during the solution of the approximate problem. 
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Based on the results of runs 1-3 * the remaining runs for the prob- 
lem were made using the invariance assumption for the element end 
forces. Run 4 (option 4(P)) is identical to run 1 except that a mixed 
variable (hybrid) approximation is employed for the stress constraints 
and the move limits are increased to 60% on the RSP's 
(d^ = 0., d^ = .6). Comparing the iteration history for run 4 with 
those of runs 1-3 shows a moderate improvement in the convergence rate 
and* considering the more liberal move limits* a significant improvement 
in maximum constraint violation history. 

In runs 5 and 6 (options 10(P) and 10(D)) the design is carried 
out in the CSD design space with 30% move limits on the CSD's 
(d^ = .3* dj — 0 . ) . Mixed variable constraint approximations are 
employed for the stress constraints. For run 5 a primal method is used 
to solve each approximate problem* while a dual method is used in run 6. 
Both runs yield essentially the same iteration histories* comparing 
favorably with those generated via the RSP design space option (runs 1- 
4) • Note that the maximum constraint violations for the intermediate 
designs are quite small with nearly all of these designs being feasible. 
Also* it is interesting to note that the dual solution method proves to 
be more efficient than the primal even though the number of retained 
constraints is greater than the number of design variables (16 as com- 
pared to 9) • 

The final run (option 1(PU)) is the same as run 1 except that the 
approximate problem update procedure is employed. As in the previous 
problem (Two Member Frame) the approximate problem is reconstructed* 
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without recourse to structural analysis or response quantity sensitivity 
calculations, once between each complete approximate problem generation. 
Comparison of the iteration histories for runs 1 and 7 indicates an 
improvement in the convergence rate by a factor of two while maintaining 
comparable maximum constraint violations for the intermediate designs 
(except in stage 1 where the violated constraint was not retained). 
Again, convergence was achieved after only 7 structural analyses. 

The final designs and critical* constraints for runs 1-7 are 
shown, along with those of the reference solution (Ref. 52), in Tables 
15 and 16. All of these designs represent the same (intuitively 
correct) design concept in which member 2 achieves its minimum allowable 
dimensions and the loads are carried through to the supports by members 
1 and 3. All of the designs have essentially the same final material 
volume and material distribution. The somewhat smaller material volume 
(less than 1%) of the reference solution can be attributed to a slightly 
different stress constraint formulation (compare Ref. 52 and Appendix 
C). 
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8.,5. Seven Member Frame ( Problem 4) 

A seven member planar frame structure subject to two independent 
in-plane loading conditions is shown in Fig. 49. This structure is 
modelled with seven thin walled box section frame elements (type 15) 
having three sizing variables (B, H, t). All members are made of the 
same material. The framework is designed for minimum mass subject to 
limits on the vertical displacements at node 3 and the allowable 
stresses at the ends of the members. It should be noted that the refer- 
ence solution (Ref. 28) also included constraints on the lowest funda- 
mental frequency, however this constraint did not participate in the 
design process. A summary of the material properties, loading condi- 
tions, constraint allowables, initial design and bounds on the element 
sizing variables is given in Table 17. The design element is described 
in Appendix C. 

The seven member frame problem was run using seven different solu- 
tion options. The iteration history data for these runs is given in 
Table 18. The iteration history plots for runs 2-7 are shown in Figs. 
50-55. In the first three runs (options 1(P)), 2(P) and 3(P)) the RSP 
design space option and linear constraint approximations are employed. 
Each approximate problem is solved via a primal solution method using 
move limits of 40% on the CSD's and 100% on the RSP's 
(d^ = .4, d2 = 1.0.). The differences between runs 1, 2 and 3 lie in 
the assumptions made regarding the element end force variations during 
the solution of each approximate problem. In run 1 the element end 
forces were assumed to be invariant during the solution of the approxi- 
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mate problem. In this case, convergence was not attained within 20 
design stages. The strong coupling among the members in the structure 
was not adequately represented by the constraint approximations. For 
runs 2 and 3, where the local and global element end force variation 
options are used, convergence is achieved after 16 structural analyses. 
Comparison of the iteration histories for runs 2 and 3 shows signifi- 
cantly improved maximum constraint violations for run 3, especially dur- 
ing the early stages of the design, as well as a slightly improved final 
design mass. 

Based on the results of the first three runs, the remaining runs 
were made using the assumption that the element end forces are dependent 
on all design variables in the structure. Run 4 (option 6(P)) is the 
same as run 3 except that mixed variable (hybrid) constraint approxima- 
tions are employed. In this case there is no convergence rate improve- 
ment over run 3 and only a small improvement in the maximum constraint 
violation for intermediate designs. 

The CSD design space option and mixed variable constraint approxi- 
mations are used in runs 5 and 6 (options 12(P) and 12(D)). For run 5, 
40% move limits are placed on the element CSD's (d^ = .4, d 2 « 0.) and 
each approximate problem is solved via a primal solution method. In run 
6 the move limits are 30% on the CSD's (d^ = .3, dj - 0.) and a dual 
solution method is employed. Comparison of the iteration histories for 
runs 4 and 5 shows little difference in performance, between the RSP and 
CSD design space options. There is, however, an improvement in the con- 
vergence rate in run 6. It should be noted, however, that while the 
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convergence rate is improved the solution time is significantly greater 
due to the fact that there are over tvice as many retained constraints 
as design variables (and therefore the dimensionality of the dual space 
is much larger than the primal space). 

The final run (option 3(PU)) is identical to run 3 except that the 
approximate problem update procedure is employed. As in the previous 
two problems, the approximate problem is updated once between each com- 
plete approximate problem generation. Comparing the iteration histories 
for runs 3 and 7 one can see that run 7 requires 50% fewer structural 
analyses for convergence while maintaining comparable performance in 
terms of the maximum constraint violation for intermediate designs. 

The final designs and critical* constraints for all runs are 

* 

shown, along with those of the reference solution (Ref. 28), in Tables 
19 and 20. While all of the designs have nearly the same final mass 
values there is considerable difference in the material distributions 
and a slight variation in the critical constraint sets. This type of 
behavior is not uncommon and was also observed in Ref. 28. Aside from 
this, two other observations can be made here. First, even though the 
material distributions are quite different, members 5 and 6 are con- 
sistently small for all cases. Secondly, the final design which most 
closely matches the reference solution (run 2) utilizes the same element 
end force variation assumption used to generate the reference solution. 
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£.6 Portal Frame ( Problem 5 ) 


Figure 56 depicts a three member planar frame subject to two 
independent loading conditions* All of the members are made of the same 
material and have the same cross section (symmetric I section with six 
sizing variables (B^, B 2 , H, t^, t 2 * * 3 ))* The structure is designed 
for minimum material volume subject to constraints on the lateral dis- 
placement and in-plane rotation at node number 3, stress constraints at 
the ends of each member and local buckling constraints for the web and 
flanges of the members* A summary of the material properties, loading 
conditions, constraint allowables, initial design and bounds on the siz- 
ing variables is given in Table 21* The design element is described in 
Appendix C. 

This problem was solved using five solution options. The itera- 
tion history data for these runs is given in Table 22* The correspond- 
ing iteration history plots are shown in Figs* 57-61. In the first 
three runs (options 1(P), 2(P) and 3(P)) the design is carried out in 
the RSP design space (note that since the cross section has six CSD r s 
the actual design variables include 4 RSP' s and 2 CSD's) using linear 
approximations of the constraints. Each approximate problem is solved 
via a primal solution method with move limits of 40% on the CSD's and 
60% on the RSP's (d^ = .4, ^2 = * 6 )* Comparing the iteration histories 
for runs 1-3 one can clearly see that the superior results (in terms of 
the number of analyses required for convergence) are obtained in run 3 
where the element end forces are assumed to be dependent on all of the 
structural design variables. This is not unexpected since the structure 
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clearly has two competing load paths. While run 3 exhibits the best 
convergence , it does result in a final design which has a material 
volume 4% greater than that obtained in run 1. Again, this is not too 
surprising as this problem is known to possess several local minimum 
solutions (Ref, 13), 

Based on the results of runs 1-3, runs 4 and 5 (options 12(P) and 
12(D)) were made using the global element end force variation option. 
The CSD design space option and mixed variable (hybrid) behavior con- 
straint approximations are employed. In both cases 40% move limits are 
placed on the element CSD 's (d^ = .4, d^ = 0,). Each approximate prob- 
lem is solved via a primal solution method in run 4 and a dual solution 
method in run 5. Comparison of the iteration histories for these runs 
reveals relatively little difference in convergence rate and overall 
maximum constraint violation for the intermediate designs. The final 
design material volume for run 5 is slightly less than that of run 4 but 
the design is also slightly infeasible. Comparing results of run 4-5 
with runs 1-3 reveals a significant difference in the final design 
material volumes between those obtained from the CSD design space option 
and those obtained from RSP space (10.2% - 13.6%), While, again, this 
is not unexpected for this problem it is interesting to note that alter- 
native design was obtained as a result of a change in the design space. 

The final designs and critical* constraints for runs 1-5 are 
given, along with those of the reference solution (Ref. 13), in Tables 
23 and 24. Comparison of these final designs, both in terms of final 
material volume and material distribution, clearly reveals the existence 
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of at least two local minimum solutions associated with distinct design 
concepts. In the case of the reference solution and runs 1-3 both pri- 
mary load paths (members 1-2, member 3) contribute significantly to the 
load carrying capacity of the structure. On the other hand, in runs 4 
and 5 the load path through member 3 is clearly abandoned in favor of 
the apparently more efficient path through members 1 and 2. 

8.7 One Bay / Two Story Frame ( Problem 6) 

A one bay/ two story frame structure subject to two independent 
loading conditions is shown in Fig. 62. This structure is modelled with 
eight thin walled tube elements (type 14), each haying two sizing vari- 
ables (R, t) . All of the members are made of the same material. This 
structure is designed for minimum weight subject to two independent sets 
of constraints. The first set (Case A) includes stress constraints at 
both ends of each member, side constraints on the element sizing vari- 
ables and constraints against both local and column buckling of each 
member. The local buckling constraints are applied in the form of upper 
bounds on the member R/t ratios. The second set of constraints (Case B) 
includes all of the constraints in Case A with the addition of con- 
straints on the lateral displacements at node numbers 2—7. It should be 
noted that while these two sets of constraints are quite similar to 
those used in generating the reference solution (Ref. 21), they are not 
identical. Additional constraints on the stresses and displacements at 
intermediate points along the members were included in Ref. 21, however, 
none of these constraints were reported as being critical in the final 
design. Also, the effective column length parameters used in the column 
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buckling constraint calculation were periodically updated in Ref. 21 
while these parameters were held constant (equal to the values 
corresponding to the final design of Ref. 21) in this work. A summary 
of the material properties, loading conditions, constraint allowables, 
initial design and bounds on the element sizing variables is given in 
Table 25. The design element is described in Appendix C. 

£.2«i Case A: Stress . Buckling and Side Constraints 

The solution for Case A was obtained using five different solution 
options. The iteration history data for these runs is given in Table 
26. The corresponding iteration history plots are shown in Figs. 63-67. 
The first two runs (options 1(P) and 3(P)) for this problem are made 
using the RSP design space option and linear constraint approximations. 
The approximate problems are solved via a primal solution method with 
40% move limits on the RSP's (d^ =0., d 2 = .4). Comparison of the 
iteration histories for runs 1 and 2 shows a slight improvement in the 
maximum constraint violation for the intermediate designs when the glo- 
bal element end force variation option is employed (run 2). It is, how- 
ever, worthwhile to note the considerable increase in the solution time 
resulting from the required element end force response quantity calcula- 
tions. 


Run 3 (option 6(P)), in this case, is the same as run 2 except 
that mixed variable (hybrid) approximations are used for the behavior 
constraints and the move limits are increased to 50% on the RSP's 
= 0., d 2 = .5). Comparing these results with run 2 indicates little 
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significant difference in terms of either convergence rate or maximum 
constraint violation history. 

The final two runs for this problem utilize the CSD design space 
option and mixed variable behavior constraint approximations. The 
approximate problems are solved via a primal method in run 4 and a dual 
method in run 5 with move limits of 30% (d^ = .3, dj c OJ and 15% 
(d^ = # 15, dj s 0.) on the element CSD's, respectively. Comparing the 
iteration histories for these two runs one can observe significantly 
poorer performance in the case where the dual solution method was used. 
This can basically be attributed to large numbers of retained con- 
straints (as many as 5-6 times the number of design variables) resulting 
in poor convergence of the optimizer. As a result of the optimizer con- 
vergence problems, relatively small move limits were required leading to 
an increase in the number of design stages needed for overall problem 
convergence. Comparison of runs 3 and 4 reveals that, although the 
design spaces are difference, there is little difference in the conver- 
gence rate. There is, however, significant difference in the maximum 
constraint violation histories, especially during the early stages of 
the design process. 

♦ 

The final designs and critical constraints for this case are 
given, along with those of the reference solution (Ref. 21), in Tables 
27 and 28. All of the designs have essentially the same final weight 
with the largest difference (between strictly feasible designs) occur- 
ring between the reference solution and run 2 (less than 1.4%). The 
critical constraints are also nearly identical. It is interesting to 
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note that the material distributions for all of the final designs are 
remarkably similar considering the considerable variations that one 
often encounters for frame type problems. This behavior is most likely 
due to the fact that the R/t constraints are critical for all of the 
members in the structure thereby effectively reducing the design freedom 
to one variable per member. 

8.7.2 Case B: Displacement . Stress . Buckling and Side Constraints 

This case was solved using five different solution options. The 
iteration history data is given in Table 29 and the corresponding itera- 
tion history plots are shown in Figs. 68-72. Runs 1 and 2 (options 1(P) 
and 3(P)) are made using the RSP design space option and linear con- 
straint approximations. Each approximate problem is solved via a primal 
solution method with move limits of 40% on the RSP's (d^ = 0 . , dj = .4). 
The results for these runs are essentially the same, both in terms of 
convergence rate and maximum constraint violation for the intermediate 
designs. It is important to note that, contrary to Case A where an 
improvement in the maximum constraint violation history was observed 
when the global element end force variation option was employed (compare 
runs 1 and 2, Table 26), here there is little constraint violation 
improvement for run 2 as compared to run 1. This is due to the fact that 
the column buckling constraints play only a small role in the design 
process (i.e., the design is essentially displacement and R/t con- 
strained, see Table 31). 
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Based on the resalts of runs 1 and 2, the remaining runs for this 
case are made using the invariance assumption for the element end 
forces. Run 3 (option 4(P)) is the same as run 1 except that the 
behavior constraints are approximated via a mixed variable approximation 
and the move limits are increased to 50% on the RSP's 
(d^ = 0., dj ■ .5). Comparison of the iteration histories for runs 1 
and 3 reveals an improvement in the maximum constraint violation history 
for run 3 (considering the more liberal move limits) but no substantial 
difference in the convergence rate. 

The final two solutions (options 10(P) and 10(D)) are obtained 
using the CSD design space option and mixed variable behavior constraint 
approximations. The approximate problems are solved via a primal method 
in run 4 and a dual method in run 5 with move limits of 30% 
(d^ = .3, dj e 0.) and 15% (d^ = .15, = 0.) on the element CSD's. As 
in Case A, the use of the dual solution method yields poorer results due 
to the large numbers of retained constraints (3-4 times the number of 
design variables). Also, as in Case A, run 4 compares well with run 3 
in terms of convergence rate but has a less attractive maximum con- 
straint violation history. 

The final designs and critical* constraints for runs 1-5 are 
given, along with those of the reference solution (Ref. 21), in Tables 
30 and 31. Again, as in Case A, all of these designs have essentially 
the same final weight, material distribution and critical constraint 
sets . 
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8.8 2x5 Grillage ( Problem 7) 

Figure 73 depicts a 2 x 5 grillage subject to a single loading 
condition. This structure is modelled with thin walled box section 
frame elements (type 11) each having four sizing variables 
(B, H, t fc , t^) . All members are made of the same material. Since both 
the structure and its loading are symmetric a half model can be used in 
solving the design problem. The grillage is designed for minimum 
material volume subject to two independent sets of constraints. In the 
first case (Case A) constraints are imposed on the vertical displace- 
ments at nodes numbers 4, 7 and 10 and side constraints are placed on 
the element sizing variables. The second set of constraints (Case B) 
includes all of those in Case A with the addition of member stress and 
local buckling constraints. A summary of the material properties, load- 
ing conditions, constraint allowables, initial design and bounds on the 
element sizing variables is given in Table 32. The design element is 
described in Appendix C. 

8,. 8.1. Case A: Displacement and Side Constraints 

This problem was solved using five different solution options. 
Since there are no stress or buckling constraints considered for this 
case all of the solution options used here utilize the element end force 
invariance assumption. The iteration history data for these runs is 
given in Table 33. The corresponding iteration history plots are shown 
in Figs. 74-78. The first two runs (options 1(P) and 4(P) ) for this 
case are made using the RSP design space option. In run 1 linear con- 
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straint approximations are employed while in rnn 2 mixed variable 
(hybrid) approximations are constructed for the displacement con- 
straints. Each approximate problem is solved via a primal solution 
method with 40% move limits on the RSP's (d^ = 0., d^ = ,4). Comparison 
of the iteration history data for these runs reveals a slight improve- 
ment in the convergence rate for run 2 but no significant improvement in 
the maximum constraint violation for the intermediate designs. 

In runs 3 and 4 (options 10(P) and 10(D)) the design is carried 
out in the CSD design space using mixed variable approximations for the 
displacement constraints. The approximate problems are solved via a 
primal method in run 3 and a dual method in run 4 with 40% move limits 
on the CSD's. The results for these runs compare well with those of 
runs 1 and 2 both in terms of convergence rate and maximum constraint 
violation. The final design material volume is slightly smaller (.3 - 
1.0%) for the designs generated using the CSD space option. The final 
run (option 1(PU)) for this problem is the same as run 1 except that the 
approximate problem update procedure is employed. Here, the approximate 
problem is updated, without recourse to structural analysis or response 
quantity sensitivity calculations, twice between each complete approxi- 
mate problem generation. Use of the update procedure leads to an 
improvement in the convergence rate by a factor of 2-2.7 while at the 
same time yielding a slightly improved constraint violation history and 
a final design material volume only 2.9% greater than run 1. It should 
be noted here that, due to the absence of stress and buckling con- 
straints, it is possible, in this case, to update the approximate prob- 
lem more than once between complete approximate problem generations 
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thereby making greater use of the high quality nodal displacement sensi- 
tivities. 

The final designs and critical constraints for this case are 
given, along with those of the reference solution (Ref. 29) in Tables 34 
and 35. In all runs the critical constraint set at the final design 
includes the vertical displacements at nodes 7 and 10. In runs 3 and 4 
the vertical displacement at node 4 is also critical. The material 
volume for the various final designs varies only slightly (the maximum 
variation is less than 3%), however, the corresponding material distri- 
butions are significantly different. Again, this result is not unex- 
pected and can be attributed to a "flatness" of the design space in the 
neighborhood of the optimum design which has been observed to exist for 
many statically indeterminate problems. 

8.. 8. 2 Case B: Displacement , Stress , Buckling and Side Constraints 

This problem was solved using four different solution options. 
Since it is expected that the local buckling constraints will play an 
important role in the design process and that the element end forces are 
strongly dependent on all of the problem design variables, all of the 
solution options employ the global force variation option. The itera- 
tion history data for these runs is given in Table 36. The correspond- 
ing iteration history plots are shown in Figs. 79-82. In the first two 
runs (options 3(P) and 6(P)) the design is carried out in the RSP design 
space using linear and mixed variable (hybrid) behavior constraint 
approximations, respectively. Each approximate problem is solved via a 
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primal solution method. In run 1 move limits of 20% and 40% 

= .2, &2 = .4) are placed on the CSD's and RSP's, respectively, 
while in run 2 the move limits are 20% and 50% (d = .2, d. = .5). Com- 
paring the results for runs 1 and 2 shows little difference in the con- 
vergence rate and only a small improvement in the maximum constraint 
violation history when the mixed variable approximation is used (run 2). 

In run 3 (option 12(P)) the design is carried out using the CSD 
design space option and mixed variable approximation for the behavior 
constraints. The approximate problems are solved via a primal solution 
method with move limits of 30% on the element CSD's (d^ *= .3, d 2 - 0.). 
The dual solution option is not used here due to the large number of 
retained constraints. The iteration history data for this run compares 
well with that of run 2 both in terms of the convergence rate and 
overall maximum constraint violation. 

The final run (option 3(PU)) for this problem is the same as run 1 
except that the approximate problem update procedure is used to update 
the approximate problem once between each complete approximate problem 
generation. The convergence rate for this run is improved over that for 
runs 1-3 , however, the maximum constraint violation history is slightly 
less attractive. 

The final designs and critical constraints for this case are 
given, along with those of the reference solution (Ref. 29), in Tables 
37 and 38. Again, as in Case A, the material volumes corresponding to 
the final designs are quite close, with the maximum variation being 
approximately 5.6%. Also, as was the previous case, there are 
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significant variations in the material distribution and slight differ- 
ences in the critical constraint sets. 

8.9 Two Bay / Six Story Frame ( Problem 8.) 

A two bay/six story planar frame structure is shown in Fig. 83. 
This structure is subjected to two independent loading conditions con- 
sisting of both concentrated nodal loading and uniform loading of the 
structural members. The structure is modelled with 30 thin walled tube 
elements (type 14), each having two sizing variables (R,t). All of the 
members are made of the same material. This structure is designed for 
minimum weight subject to two independent sets of constraints. The 
first set (Case A) includes stress constraints at both ends of each 
member, side constraints on the element sizing variables and constraints 
against both local and column buckling of each member. The local buck- 
ling constraints are applied as upper bounds on the member R/t ratios. 
The second set of constraints (Case B) includes all of the constraints 
in Case A with the addition of constraints on the lateral displacements 
at node numbers 1-18. As was mentioned in the previous discussion of 
Problem 6 , these two sets of constraints are not identical to those used 
in generating the reference solution (Ref. 21). Additional constraints 
on stresses and displacements at intermediate points along the members 
were included in Ref. 21, however, none of these constraints were 
reported as being critical in the final design. Also, the effective 
column length parameters used in the column buckling constraint calcula- 
tions were periodically updated in Ref. 21, while, in this work, these 
parameters were held constant (equal to the values corresponding to the 
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final design of Ref. 21). A summary of the material properties, load- 
ing conditions, constraint allowables, initial design and bounds on the 
element sizing variables is given in Table 39. The design element is 
described in Appendix C. 

8.9.1 Case A: Stress . Buckling and Side Constraints 

The solution for Case A was obtained using five different solution 
options. The iteration history data for these runs is given in Table 
40. The corresponding iteration history plots are shown in Figs. 84-88. 
The first two runs (options 1 (P ) and 3(P)) for this problem are made 
using the RSP design space option and linear constraint approximations. 
Each approximate problem is solved via a primal solution method with 40% 
move limits on the element RSP's (d^ = 0., dj a .4). Comparison of the 
iteration histories for runs 1 and 2 shows a moderate improvement in 
both the convergence rate and maximum constraint violation for the 
intermediate designs when the global element end force variation option 
is used (run 2). However, it is important to note the considerable 
increase in solution time resulting from the required element end force 
sensitivity calculations. 

Run 3 (option 6(P)) is the same as run 2 except that mixed vari- 
able (hybrid) approximations are used for the constraints and the move 
limits are increased to 50% on the RSP's (d 1 =0., d- = .5). Comparing 
these results with run 2 reveals a decrease in the number of analyses 
required for convergence and an improvement in the overall maximum con- 
straint violation history. However, the final design weight for run 3 
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is 1.4% greater than that of run 2. 


The fourth run (option 12(P)) utilizes the CSD design space option 
and mixed variable approximations for the behavior constraints. The 
approximate problems are solved via a primal solution method with 40% 
move limits on the CSD's (d^ = .4, dj s 0.). Comparison of the itera- 
tion history data for this run with that of run 3 reveals an increase in 
the number of analyses required for convergence and a slight deteriora- 
tion in the maximum constraint violation history. The final design 
weight is, however, approximately 1% lower than that of run 3. 

The omparing approximate problem update procedure is used. Here, 
the approximate problem is updated once between each complete approxi- 
mate problem generation. In this case the use of the update procedure 
results in only a slightly improved convergence rate. 

The final designs and critical* constraints for this case are 
given, along with those of the reference solution (Ref. 21), in Tables 
41 and 42, All of the final designs have essentially the same weight 
with the largest difference occurring between runs 3 and 5 (approxi- 
mately 1.6%). The material distributions and critical constraint sets 
are also very similar. 
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8.. 9. 2 Case B: Displacement . Stress . Buckling and Side Constraints 

Case B was solved using five different solution options. The 
iteration history data for these runs is given in Table 43. The 
corresponding iteration history plots are shown in Figs. 89-93. The 
first two rnns (options 1(P) and 3(P)) for this problem are made nsing 
the RSP design option and linear constraint approximations. Each 
approximate problem is solved via a primal solution method with 40% move 
limits on the ESP's (d^ = 0.. dj = .4). Comparison of the iteration 
histories for rnns 1 and 2 reveals a slight increase in the number of 
analyses required for convergence and a small improvement in the maximum 
constraint violation for the intermediate designs when the global ele- 
ment end force variation option is used (run 2). It is important to 
note that, contrary to Case A where the overall design process was 
improved through the use of the global element end force variation 
option (compare runs 1 and 2, Table 40), here there is little signifi- 
cant improvement. This is basically due to the fact that in Case B the 
displacement constraints (whose approximations do not depend on the 
choice of the element end force variation option) play and important 
role in the design process (see Table 45). 

Based on the results of runs 1 and 2, the remaining runs for this 
case are made using the element end force invariance assumption. Run 3 
(option 4(P)) is the same as run 1 except that mixed variable (hybrid) 
approximations are utilized for the behavior constraints and the move 
limits are increased to 50% on the element RSP's (d^ « 0., d^ = .5). 
Comparison of the iteration histories for runs 1 and 3 reveals a signi- 
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ficant improvement in the maximum constraint violation history for run 3 
(especially when considering the more liberal move limits) but little 
difference in the convergence rate. 

The fourth run (option 10(P)) utilizes the CSD design space option 
and mixed variable approximations for the behavior constraints. The 
approximate problems are solved via a primal solution method with 40% 
move limits on the element CSD's (d^ = .4, d^ ** 0) . Comparison of the 
iteration history data for run 4 with that for runs 1-3 reveals a signi- 
ficant improvement in the convergence rate for run 4 while maintaining a 
comparable maximum constraint violation history. 

The final run for Case B is the same as run 1 except that the 
approximate problem update procedure is employed. Here, the approximate 
problem is updated once between each complete approximate problem gen- 
eration. Comparing the results of run 5 with runs 1-3 reveals a sub- 
stantial improvement in the convergence rate for run 5 with only a small 
increase in the maximum constraint violation for the intermediate 
designs. 

The final designs and critical* constraints for this case are 
given, along with those of the reference solution (Ref. 21), in Tables 
44 and 45. As in Case A, all of the final designs have essentially the 
same weight with the largest difference being less than 1.6%. The 
material distributions and critical constraints for the final designs 
are also quite similar. 
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8 .10 Helicopter Tail Boom ( Problem 9) 

Figure 94 depicts a space frame idealization of a helicopter tail 
boom subject to a single loading condition. All members are made of the 
same material and have the same cross section (thin-walled tube with two 
sizing variables (R, t)). This structure is designed for minimum weight 
subject to constraints on the nodal displacements in the y and z direc- 
tions at node numbers 5-28 and side constraints on the element sizing 
variables. Stress constraints are also imposed at both ends of each 
member along with column buckling and local wall buckling constraints 
(in the form of R/t constraints). A summary of the material properties, 
constraint allowables, loading conditions, initial design and bounds on 
the element sizing variables is given in Table 46. The design element 
(type 14) is described in Appendix C. 

This problem was solved using five different solution options. 
Since this problem is expected to be essentially displacement and R/t 
constrained all of the solution options chosen make use of the element 
end force invariance assumption. The iteration history data is given in 
Table 47. The corresponding iteration history plots are shown in Figs. 
95-99. In the first two runs (options 1(P) and 4(P)) the design is car- 
ried out in the RSP design space using linear and mixed variable 
(hybrid) behavior constraint approximations, respectively. Each approx- 
imate problem is solved via a primal solution method with move limits on 
the RSP's of 60% (d x = 0., d 2 = .6) for run 1 and 70% = 0., d 2 - .7) 
for run 2. Comparison of the iteration histories for these runs reveals 
no difference in the convergence rate and only a slight change in the 
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maximum constraint violation for the intermediate designs. It is 
interesting, however, to observe that, even though the initial design is 
highly infeasible (211.6%), near feasible designs are achieved after 
only three design stages. 

In runs 3 and 4 (options 10(P) and 10(D)) the design is performed 
using the CSD design space option and mixed variable behavior constraint 
approximations. The approximate problems are solved via a primal method 
in run 3 and a dual solution method in run 4 with 40% move limits on the 
CSD's (d^ = .4, * 0.). There is little difference in convergence rate 
or maximum constraint violation for these runs, however, the dual solu- 
tion method is less efficient here due to the large number of retained 
constraints (twice the number of design variables). 

The final run (option 1(PU)) for this problem is the same as run 1 
except that the approximate problem update procedure is employed. Here, 
the approximate problem is updated, without recourse to structural 
analysis or response quantity sensitivity calculations, once between 
each complete approximate problem generation. Comparison of the itera- 
tion history for this run with those of runs 1-4 reveals a dramatic 
improvement in both the convergence rate and the maximum constraint vio- 
lation history. Also, it is important to note that the total solution 
time is improved here due to the fact that the increase in optimization 
time is more than offset by the decrease in analysis time (which 
includes the time required for the approximate problem generation). 
This improvement is expected to be much more dramatic for large practi- 
cal problems where the solution of the structural analysis problem 
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represents the primary computational bnrden. 


The final designs and critical* constraints for this problem, 
along with those of the reference solution (Ref. 21), are given in 
Tables 48 and 49. In all cases, the critical constraint set at the 
final design includes the displacement constraints at nodes 25 and 27 
and the local wall buckling (R/t) constraints for nearly all of the 
members. Also, the final design weight and material distribution is 
nearly the same for all runs. Finally, it is interesting to observe the 
intuitively satisfying result whereby the lighter designs contain 
slightly larger members at the base (fixed) end of the structure. 
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CHAPTER IX 


Conclusions and Recommendations 


9.1 Conclusions 

A synthesis methodology for the design of frame-truss structures 
subject to multiple static loading conditions has been developed. This 
methodology has been implemented in the COMPASS computer program and has 
subsequently been used to solve a variety of frame-truss synthesis prob- 
lems. The numerical results presented here illustrate the feasibility 
of obtaining near optimum designs after only 5-10 structural analyses 
for most frame-truss structures. 

While it is believed that the primary goal of this study has been 
achievedi it is important to realize that attaining this goal has 
required not only the introduction of a full gamut of approximation con- 
cept$! but also the proper application of these techniques to the prob- 
lem at hand. Unlike the truss-membrane synthesis methodology, which 
tends to employ a rather standard set of approximation concepts to the 
solution of most problems with uniform success, the key to the efficient 
solution of a frame-truss design problem, in many cases, lies in the 
thoughtful selection of appropriate approximation techniques and 
mathematical programming methods from a set of available options. 

Engineering insight and experience has long been a part of tradi- 
tional design methods and it is not surprising to find that somewhat 
analogous insights and experiences must now become a part of an 
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efficient frame-truss synthesis methodology. The considerable body of 
computational experience reported here is intended to provide the design 
engineer with the foundation needed to build expertise in the use of 
frame synthesis methodology. While most of the points that will be dis- 
cussed in the sequel have been previously mentioned in Chapter VIII, it 
is useful to summarize these results here in the form of guidelines for 
solving frame design problems. 

By now it should be apparent that the efficient solution of a 
frame-truss synthesis problem involves making important decisions about 
the construction and solution of the approximate design problems. The 
items which must be considered include the types of approximations 
employed, the choice of design space, the optimization method used to 
solve the mathematical programming problems, the choice of appropriate 
move limits and even the frequency of performing the structural analyses 
and response quantity sensitivity calculations. While all of these 
items are important, probably the most important decisions to make are 
those related to the construction of the behavior constraint approxima- 
tions, This involves deciding which design variables the constraint is 
believed to be dependent on (i.e. choosing the appropriate element force 
variation option) and selecting the basic character of the approximation 
to be used (i.e, linear vs, mixed variable). The following guidelines 
are offered: 

1, the element end force invariance option is recommended if 

a, the structure is statically determinate. 
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b. the structure is statically indeterminate but the coupling 
between the structural members is believed to be weak (see 
Problems 2 and 3, Chapter VIII). 

c. the design problem is believed to be essentially displace- 
ment constrainted (see Problems 6 and 8 (Case B) and Problem 
9, Chapter VIII). 

2. the local element end force variation option can be useful in 
cases where the invariance assumption is inadequate and the 
expense associated with the global force variance option is prohi- 
bitive. This option should be used with care as it can, in some 
cases, yield approximations which are inferior to the force 
invariance option (see Problem 3. Chapter VIII). 

3. the global element end force variation option is recommended if 
the structural behavior is strongly coupled and the design problem 
is believed to be essentially strength critical (see Problems 1, 5 
and 6 , Chapter VIII). 

4. the linear constraint approximations are recommended if the con- 
straints controlling the design process are believed to be linear 
or nearly linear functions of the design variables. 

5. the mixed variable (hybrid) approximation generally yields equal 
or superior performance as compared to the linear approximations 
(except as noted in item 4) with little additional computational 
expense . 
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the mixed variable (hybrid) approximations may, in some cases, be 
nsed in lien of higher levels of element end force sensitivity 
information vithont the added computational expense (see Problem 
2, Chapter VIII). 

The guidelines for selecting the appropriate design space option 
(i.e., either reciprocal section property (RSP) or cross sectional 
dimension (CSD) space) are not as specific as those for selecting the 
constraint approximations. For the problems solved in this study the 
overall performance of the tvo design spaces has been nearly the same, 
although distinctly different designs may be generated depending on the 
design space chosen (see Problem 5, Chapter VIII). Part of the reason 
that this has been the case can be attributed to the use of the approxi- 
mate sizing variable recovery transformation, described in Chapter VI, 
which compromises the quality of the constraint approximations in the 
RSP design space. It is possible that a more exact recovery method 
would lead to superior results for the RSP design space option. For the 
current implementation of the frame-truss synthesis methodology one may 
wish to consider the following when making the design space selection: 

1. the use of the RSP design space generally requires the addition of 
sizing variable side constraint approximations to the approximate 
design problem. Vhile this has not been a computational burden 
for the problems solved here it may be burdensome for problems 
involving large numbers of sizing variables. 

2. the selection of move limits in the CSD design space may be more 
physically meaningful than in RSP space. 
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The selection of the mathematical programming method (i.e., primal 
or dual) for solving the approximate problems is relatively straight 
forward. Generally, if the number of retained constraints is expected 
to be less than 1.5 times the number of design variables the dual method 
will be at least as, if not more, efficient than the primal method (see 
Problems 2, 3 and 5, Chapter VIII). For problems where the number of 
retained constraints is expected to be large compared to the number of 
design variables the primal solution method should be used (see Problems 
4, 6 and 9, Chapter VIII). Two additional comments are appropriate 
here. First, the development and implementation of a specialized dual 
method (in which the dimensionality of the dual space does not exceed 
the number of truly critical constraints) could make the dual method 
more efficient even when the number of retained constraints is large. 
Lastly, it should be recognized that for large problems, where the 
structural analysis and approximate problem generator can be expected to 
represent the main computational effort in the design process, the 
selection of the optimization method is not expected to seriously effect 
the efficiency of the overall problem solution. 

The selection of appropriate move limits on the structural design 
variables can be quite difficult in the absence of prior computational 
experience. Even with experience the selection process is often problem 
dependent and may require several attempts. Move limits which are 
either too tight or too loose can seriously effect the problem conver- 
gence rate or, in some cases, preclude convergence altogether. For the 
problems solved in this study the move limits ranged from 20% to 70% on 
the RSP's and 15% to 50% on the CSD's, with the most frequently used 
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values being 40-50% for the RSP's and 30-40% for the CSD's. The follow- 
ing guidelines are offered for making an initial choice of move limits: 

1. problems which are expected to have a large number of critical 
buckling constraints will generally require tighter move limits 
than those which are essentially stress and/or displacement con- 
strained. 

2. the use of the mixed variable (hybrid) constraint approximation 
option generally allows for more liberal move limits, especially 
when the problem is essentially strength critical. 

The final decision involves the use of the approximate problem 
update procedure. Based on the numerical examples presented here it can 
be concluded that the use of this option generally results in an 
improvement in the convergence rate of the design problem (i.e. the 
number of structural analyses required for convergence is reduced). For 
small problems, however, even though the convergence rate is improved 
the total solution time can increase when this option is employed (see 
Problem 2, Case B, Chapter VIII). For larger problems of practical 
interest the use of this option does result in savings in solution time 
as well as an improved convergence rate (see Problem 9, Chapter VIII). 
The solution time savings are expected to be more dramatic as the 
analysis problem size increases. Aside from the fact that the approxi- 
mate problem update option is best used for large problems two other 
comments are appropriate here: 

1. the approximate problem can be updated more frequently without 
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reanalysis if the design problem is essentially displacement crit- 
ical (see Problem 7, Case A, Chapter VIII). 

2. in the case vhere the design problem is essentially strength crit- 
ical. the approximate problem may be npdated more frequently 
without reanalysis when the structural behavior is weakly coupled 
than when it is strongly coupled. 

g.,2 Be commendations for Future Work 

While the frame-truss synthesis methodology presented in this work 
can be used to efficiently solve a significant class of structural 
design problems, several areas of future work and investigation can be 
identified which will broaden the applicability of the method to include 
a larger class of problems and/or lead to increased solution efficiency. 
These areas are summarized below: 

1. the expansion of the design element library to include a greater 
variety of cross sectional shapes. 

2. the addition of a modal analysis capability for the inclusion of 
frequency constraints in the design problem. 

3. the addition of an optimum design sensitivity capability. 

4. the extension of the approximate problem update procedure to 
include the updating of the stress and buckling constraint partial 
derivatives which are explicit functions of the CSD's and RSP's. 

5. investigate the use of a cumulative constraint formation (or some 
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other method) as a means of reducing the numbers of stress and 
buckling constraints retained for each structural member* 

6. develop a specialized, efficient dual solution method for solving 
all forms of the approximate problems* 

7* investigate the possibility of replacing the approximate recovery 
method with an efficient nonlinear recovery technique. 

The suggested extensions offered above can be divided into two 
groups based on the probability that the work can be successfully com- 
pleted. It is believed that items 1-4 pose relatively little risk in 
that the work basically entails the implementation of proven concepts* 
Items 5-7, however, may require considerable investigation and if com- 
pleted may not produce the desired results* 
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APPENDIX A 


Stiffness Matrices and Load Vectors for Structural Elements 


A.l Prismatic Frame Element 


The space frame element, shown in Fig. Al, is a two node, twelve 
degree of freedom element oriented with the longitudinal axis in the 
local x coordinate direction and the cross section principal axes in the 
local y and z coordinate directions. The element is assumed to have 
linear axial and torsional displacement states given by 


i(x> = [N 1 (x),N 2 (x)] jujj 


(A-l) 


e x (x) = l 6 x 


and cubic bending displacement states of the form 


v(x) = 


[N 3 (x),N 4 (x),N 5 (x).N 6 (x)]/ Vl 


e 


(A-2) 


(A-3 ) 


and 
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r(x) = [Ng(x) ,N^(x) ,Nj(x) ,Ng(x)] j Wj 

e_ 


e 

V y 2 


(A-4) 

where u, v and w are the displacements in the local x, y and z coordi- 
nate directions, 0 , 6 and 6 are the rotations about the x,y and z 
x y z 

axes and where the displacement shape functions are given by 


N^x) 



N 3 (x) = 1 - 3(J) 2 + 2<J) 3 

N 4 (x) = L[(J) - 2<|) 2 + (J) 3 ] 


(A-5) 


N 5 (x) = 3(J) 2 - 2(J) 3 


N 6 (x) = L[ - (J) 2 + (J) 3 ] 


The assumed displacement states (Eqs. (A-l) - (A-4)) lead to the follow- 
ing strain-displacement relations 


Axial strain: 


it - tB *i 


hi 

(”2) 
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where 


The total 


30 x / 

Torsional strain: = [B^] I 0^ 


Cnrvatnre: 


iiz 

dx 2 


[B lU ^ 


(A - 6 ) 



[B l] 


r _ i 


L V 


I 1 

L J 


[ V - [ - i (1 -T ) 'T (2 -r ) '-i ( ^ 1 ) '-l ( 1 -t>] 


(A-7) 


strain energy for the frame element can now be written as 
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< A— 8) 


„ EA h f 3u,T r 3u, - . GJ L . f fle x.T 

u = T y W fe }dx + T y { 3x } 
o o 


ae 

'I7 ld * + 


El L ,2 ^ .2 

^ y {^f) T {^f)dx + 

2 o 3x 2 3x 2 


El L -2 T -2 
-T 1 7 ( fi -f} T { L f}dx 
2 o dx 2 9x 2 


- \ {q} T [K] e {q} 

where 


E 

G 

A 

I 

I 


and 


material modulus of elasticity 

material shear modulus 

cross sectional area 

cross sectional principal moment of 
inertia about the y axis 

cross sectional principal moment of 
inertia about the z axis 


U) ■= (aj Vj »j 6 


V 9 I, °2 v 2 \ 9 x, 


e 


0 } 
y 2 Z 2 


Evaluation of the integrals of Eq. (A— 8 ) yields the following form of 
the element stiffness matrix in local coordinates 
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(A-9) 


Similarity, the external work expression for the space frame element 
subject to uniformly distributed loads can be written as 


L 

W = J 
o 



m 

x 


m 

y 
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( A— 1 0 ) 


T 

= CP}® (q) 


where p , p and p are forces per unit length in the x, y and z direc- 
x y z 

tions and m , m , m are moments per unit length about the x, y and z 
x y z 

axes. Evaluation of Eq. (A-10) using Eqs. ( ( A— 1 ) - (A-5)) and assuming 

p , p , p as well as m , m , m , are constants (uniformly distributed) 
x y z x y z 

leads to the following form for the work equivalent load vector 



m 





+ m 

7 



m L 
x 

2 ' 


m L 




12 ' 



(A-ll) 


A. 2 Prismatic Truss Element 

The space truss element, shown in Fig, A2, is a two node element 
oriented with the longitudinal axis in the local x direction. The ele- 
ment is assumed to have a linear displacement state of the form 


n(x) = [1 - f} 



(A-12) 


where u is the displacement in the x direction. This assumed 
ment state leads to the following strain-displacement relation 


displace- 


du(x) _ rZ l ll 
ax 1 L # L J 



The strain energy for the truss element can now be written as 


(A-13 ) 
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u =f 3 { fJ }T lfl }dl 

o 


- 2 { °1 S 11 * 1 ' ( tt lj 

( 2 ) (A-14) 
which leads to the following expression for the element stiffness matrix 
in local coordinates 


[K] e = 


= EA 1 1 

“ L L-i i. 


( A-15) 

Similarly, the external work due to a uniformly distributed axial load 

p can be written as 
x 


V=Jp. u ( x ) dx = {P} 


si 


which leads to the following work equivalent load vector 


(A-16) 


e T 

{P} 6 = 2 {1 11 


(A-17) 


Rewriting [K] e and {P } e in terms of the full twelve nodal degrees of 
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riw 


freedom gives 


and 



(A-18) 


T 


{PJ® 


P X L 

{1 0000010000 0 } 


(A-19) 
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APPENDIX B 


Generalized Inverse Transformation 


Consider the m-vector {X} and the n-vector {Y} related by the m x 
n linear transformation matrix [J] as follows: 


[JltYJ = {X} 


(B-l) 


For the case where m = n and [J] is non-singnlar the inverse transforma- 
tion is clearly given by 


{1} = [J] 1 {XJ 


(B-2) 


For the general case where m ^ n and the rows of [J] are not necessarily 

linearly independent, an inverse transformation of the form given by Eq. 

(B-2) does not generally exist. It is, however, possible to determine a 

T 

set of linearly independent (free) variables [{X^J, {Yp} ] and a set of 

T 

}] and to construct a linear 
transformation between them having the form 


dependent (basic) variables [{Xg}, {Y 



(B-3) 


The transformation matrix [H] can be constructed from [J] via two sets 
of pivot operations as shown in the following construction, which relies 
heavily on material contained in Ref. 54. 


Consider Eq. (B-l) and imagine that it is partitioned in terms of 
the free {Yp} and basic {Yg} variables as follows: 
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[B F] 


{I R } 

L { Vj 


= {X} 


(B-4) 

Pivoting to determine the basis inverse and allowing for the existence 
of linearly dependent rows gives 


[; H ■ ra » 


where 


(B-5) 


[B 1 ] = 


Eq. (B-5) yields the following two equations: 


{Y b J + [B X F] {Ip} = [B 1 ](X) 


( B— 6 ) 


(B-7) 


IB 2 ](X} = {0} 


(B-8) 


The existence of Eq, (B-8) indicates that the {X} variables are not 
linearly independent. Therefore, partitioning Eq. (B-8) in terms of the 
free and basic {X} variables gives 


k I %][£!] ■ 


(0) 


(B-9) 


Pivoting yields 


[ r i b " b %] [S3 ■ 


{0} 


(B-10) 


Solving Eq. (B-10) for {Xg} gives 
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(B-ll) 


-1 


tX B } " " [B 2 ®2 HX F } = t C HV 

B F 


We now have {Xg} expressed in terms of {Xp} . To express {Yg} in 
terms of the free variables let 


[B 3 = [B 1 B, ] 
B F 


(B-12) 


and 


tX} 


L {X F } J 


(B-13) 


Substituting Eqs. (B-12) and (B-13) into Eq. (B-7) and solving for {Y D } 


gives 


tYfiJ = [B ] {X B } + [B^CXpJ - [B _1 F]{Y f ) 


Substituting Eq. (B-ll) for {Xg} in Eq. (B-14) yields 

{Yg} = ( IB 1 ] + [B 1 HClHXp) - [B _1 F] {Y p ) 

F B 


Finally, writing Eqs. (B-ll) and (B-15) in matrix form yields 


(B-14) 


(B-15) 




L (y b'j 


- [C] I [0] 

[8 ] ♦ IB L] [C] | . [B -1 F] 

F D 



'«F>' 


;<v 


Rxp.fi rtxpji 

M ‘ [H1 M 


(B-16) 


(B-17) 

It should be noted that if Eq. (B-l) is written in terms of the variable 
perturbations (AX) and {AY} as 
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[ J] {AY} = {AX} 


( B— 18) 


where 


[j] = rin 

LJJ LdiJ 


then Eq. (B-17) becomes 


r{AX B }] 

L tAY B } J 


[H] {AZ} 


(B-19) 


(B-20) 


where 


{AZ} = [{AXp} {AYp} ] 


and 


[H] 


dZ 

!!b 

L az . 


(B-21) 


(B-22) 
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APPENDIX C 


Design Element Library 

The design element library contains the descriptions of the vari- 
ous design elements available for constructing the structural design 
model. Each element is completely described by the following: 1) a 

basic structural element type (e.g. frame or truss), 2) a description of 
the element cross section and 3) a set of element level constraints 
(e.g. stress, local buckling). The element end forces used in the ele- 
ment level constraint calculations are shown in Fig. Cl for both the 
frame and truss type elements. 

C.l Design Element Descriptions 

The following design elements are available and are described in 
this section. 


Element Type 

Truss with one sizing variable 1 
Box beam with four sizing variables 11 
I-beam with six sizing variables 12 
Square beam with one sizing variable 13 
Thin walled tube with two sizing variables 14 
Box beam with three sizing variables 15 
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Element: TRUSS 


Type : 1 

Description ; This element is a truss element with the cross 
described by the sizing variable A as shown in Fig. C2. One 
constraint is computed for this element. 

Reciprocal Section Property 
\ = 1/A 


Stress Constraint 


g 



l F X *1 I 

a 

a 


- 1 < 0 


where a 

a 


is the allowable stress in tension and compression. 


section 

stress 


(C-l) 


(C-2) 
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Element : BOX BEAM 


Tvne : 11 

Description : This element is a frame element with the cross section 
described by the sizing variables B, H, t^ and t^ as shown in Fig. C3 . 
Eight stress and four local buckling constraints are computed at each 
end of the element for a total of sixteen stress and eight local buck- 
ling constraints. The locations on the cross section at which these 
constraints are evaluated are shown in Fig. C3 for the first node (n^) 
end of the element. 


Cross Sectional Dimensions 


Y, = B 


Y, = H 


Y 3 * *h 


Y 4 = S 


(C-3) 


Rec iprocal Section Properties 


X l~ A~ YjTj - <y 2 


2Y 4 )(Y 1 - 2Y 3 ) 


_ 1 Y 2 Y 4 + Y 1 Y 3 

2 " J " 2Y 2 Y 2 Y Y 
ZI 1 J 2 *3 *4 


(C-4) 
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^ c JL = Afr 

3 Iy - (I 2 - 2 V (I 1 - 2I 3 )3 

x = 1_ m 12 

4 Ij *1*2 - ( V 2 V ,t 2 - 2 V 3 

Stress Constraints 

The stress constraints for this element are based on the Von Mises 
criterion and have the form 

g £ = -y (a\ + 3x\) - 1 i 0 ; 1 = 1.2, ...16 

a (C-5) 

where a is the allowable stress. The normal stresses (a.) are obtained 
a i 

from technical beam theory and are given by 

- V 2 M z Y 2 X 4 - V 2 M y Y 1 X 3 + F x X 2 ; i = 1 , 2 , 9,10 

V 2 M z Y 2 X 4 - V 2 MyY^ + ; i = 3 , 4 , 11,12 

V 2 M z Y 2 X 4 + ty 2 MyY^ + F x Xj ; i = 5,6,13,14 

- V 2 M z Y 2 X 4 + V 2 M y Y 1 X 3 + F x Xj ; i = 7,8,15.16 

(C-6) 

The shear stresses are calculated assuming that the box beam is thin 
walled (i.e. horizontal walls do not resist vertical shearing forces and 
vice-versa; shear stresses due to shearing forces are uniformly distri- 
buted over the appropriate wall areas; shear stresses due to the twist- 
ing moment is uniformly distributed over the cross section) and are 
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given by 


2 W4 2Y 1 Y 4 


M_ 


2 Y lV 3 2Y 2 Y 3 


M 


F. 


2Y 1 Y 2 Y 4 


2Y 1 Y 2 Y 3 


2Y Y 
X1 1 J 4 


2Y 2 Y 3 


; i = 1,8,9,16 

; i - 2,3,10,11 

; i - 4,5,12,13 

; i = 6,7,14,15 


(C-7) 


Local Beetling Constraints 

To protect against local buckling of the box beam walls each side 
of the member is conservatively modelled as a simply supported infin- 
itely long plate, subject to combined axial, bending and shear stresses 
as shown in Fig. C4. Defining the buckling stress ratios 


R = t/t 


cr 


*b = c b /a bcr 


R = a /a 
x x xcr 


(C-8) 


where, for an infinitely long plate of width b and thickness t (Ref, 55) 
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x - 5.35 S 
cr 

a, =23.8 S 
bcr 


o 


xcr 


4.0 S 


(C-9) 


S = E g — (t/b 2) 

12(1- v 2 ) 


- C (t/b) 2 

v 


and combining Eqs. (C-8) in one interaction formula leads to the follow- 
ing local buckling constraint equation 


g = R x + + E 2 - ! i 0 


(C-10) 


Eq. (C-10) is evaluated at four points on each end of the element for a 
total of eight constraints. The expressions for cr^, cr^, x and S for 
each of these constraints is given below: 



- V 2 M z T 2 X 4 - F A ; 

V 2 VjX,, - f A i 

Wl '.Vi ! 

- V 2 MyYjY, - F^ ; 


j = 1.5 
j = 2.6 

j = 3.7 
j - 4.8 


(C-ll) 

where a > 0 is a compressive stress. 

X j 
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(C-12) 


(C-13) 


r 


; j » 1.3. 5. 7 


; j = 2. 4, 6, 8 


(C-14) 
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Element : I-BEAM 


Type : 12 

Description : This element is a frame element with the cross section 
described by the sizing variables B^, tj, B 2 , tj. H, and t 3 as shown in 
Fig. C5. It is intended primarily for use in planar frame structures. 
Three stress and three local buckling constraints are computed at each 
end of the element for a total of six stress and six local buckling con- 
straints. The locations on the cross section at which these constraints 
are evaluated are shown in Fig. C5 for the first node (n^) end of the 
element. 



( C-l 5 ) 
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Reciprocal Section Properties 


x 1 _ 1 = 1 

1 * Y 1 Y 2 + Y 3 T 4 + T 5 Y 6 ^ + ^2 + ^3 

** 1 vl + V« + 

x = _L = 12 

3 Iy ♦ Vs - h*6 <c_16) 

x B _l_ = l 

4 ** * 1 ^ + A + h A * 126 


where 


Y Y 

B = V y " ^ )2 + V y ■ Y 2 " Y 5 " 2 t)2 + V y “ T 2 

_ A 1 (Y 2 /2) + ^(Yj + Y 5 + Y 4 /2) + Ag(Y 2 + Y 5 /2) 

y = A 



(C-17) 


Stress Constraints 


Constraints on normal stress for the flanges and on shear stress 
for the web are applied at each end of the element. These constraints 
have the form 


i 0 ii 


- 1 < 0 ; i = 1.2, 4, 5 


8 i = 


(C-18) 


1^1 


-110 ; i = 3 ,6 


where a and x are the normal and shear stress allowables. The ezpres- 
a a 
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sions for and x ^ are obtained from technical beam theory and are 

given by 


" M z (Y 2 + Y 4 + Y 5 ■ y)X 4 + F x X 1 ; i= 1,4 


°i = 


H yX. 
z 9 4 


; i = 2.5 


and 


T i * V4 


- A 

- Y v + — 
Y, 21, 2 y 2 


g T iV vj 


; i = 3,6 


(C-19) 


(C-20) 


It shoald be noted that since I-sections are generally not designed to 
carry twisting moments, bending moments about the vertical axis or hor- 
izontal shear forces, the stress computations do not include these 
effects. Therefore this element should not be used for structures in 
which these effects are expected to be significant. 


Local Buckling Constraints 

To protect against local buckling of the element flanges and web 
local buckling stress constraints are evaluated for both the flanges and 
the web at each end of the element. The flanges are assumed to experi- 
ence only normal compressive stresses and each half of the flange is 
conservatively modelled as an infinitely long plate, simply supported at 
the ends and along one side, and free along the remaining side. The web 
is assumed to be subject only to shear stresses and is conservatively 
modelled as an infinitely long simply supported plate. The buckling 
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constraints have the form 


— — - 1 i 0 ; i = 1,2,4, 5 

a 


(C-21) 


l^il 

— - 1 i 0 ; i = 3,6 

X 

cr i 


where o and v _ are given by (Ref. 56) 
cr. cr. 



where E is the material modulus of elasticity and is the Poissons 

ratio. 
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Element : SQUARE BEAM 


Type : 13 

Description : This element is a frame element vith the cross section 
described by the sizing variable B as shown in Fig. C6. Four normal 
stress and fonr shear stress constraints are calculated at each end of 
the element for a total of sixteen stress constraints. The locations on 
the cross section at which these contraints are evaluated are shown in 
Fig. C6 for the first node (n^) end of the element. 


Cross Sectional Dimensions 


Y 


1 


= B 


(C-22) 


Reciprocal Section Properties 



X, 



X 


4 



(C-23) 
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Stress Constraints 


The normal and shear stress constraints for this element have the 

form 


8 i - 


I 


-HO 


l T ,l 


- 1 i 0 


i = 1-4.9-13 


i = 5-8,14-16 


(C-24) 


where a and t ere the normal and shear stress allowables. The ezpres- 
a a 

sions for a ^ are obtained from technical beam theory and are given by 


- ^2 - l/ 2 M z T 1 X 4 + F^j ; 

i « 1,9 

- V 2 MyYjXj + 1/2 M 2ll X 4 + i 

i = 2,10 

V 2 MyY,X 3 + V 2 M i Y 1 X 4 + F 2 Xj ; 

i = 3,11 

MyYjX 3 - V 2 M i Y 1 I 4 * F lXl ; 

i = 4,12 



(C-25) 


The shear stress is obtained by the superposition of the parabolic 
horizontal shear stress distribution from technical beam theory and a 
torsional shear stress distribution for a solid square bar (Ref. 57) . 
The expressions for x^ are given by 
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(C-26) 



where c 


- c M Y- X- + 1/8 F Y?X. ; i = 5,13 

x 1 2 y l 4 

- c M Y 1 X- + 1/8 F Y?X, ; i = 6,14 

X 1 2 z 1 o 

c M I Y 1 X 2 + 1/8 F y Y^X 4 ; i = 7,15 

c M 1,1- + 1/8 F I?X, ; i = 8,16 

X X 2 Z 1 o 

= .676 
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Element: THIN TUBE 


Tvne : 14 

Description : This element is a frame element with the cross section 
described by the sizing variables R and t as shown in Fig. C7. Two 
stress and two local column buckling constraints are evaluated at each 
end of the element for a total of four stress and four local column 
buckling constraints. These constraints are evaluated at the points of 
maximum normal stress defined by the angle 0, as shown in Fig. C7, where 
6 is given by 


0 “ arctan (M /M ) 

Z 7 (C-27) 

It should be noted that 0 is assumed to be constant during the solution 

to each approximate design problem and is updated for each approximate 

problem. In addition, one local wall buckling constraint is evaluated 

for this element. 

Cross Sectional Dimensions 

Yi = R 

I 2 = t (C-28) 

Reciprocal Section Properties 

X = — = 

1 A 2nY 1 Y 2 

X = — = 

2 J n(R 4 -R 4 ) (C-29) 

o 1 
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4 


13 ' K = .<<-<) 


JL i 

*z n(R 4 -R 4 ) 
o i 


where 


R 


* T 2 /2 


®i ■ T 1 - V 2 


(C-30) 


Stress Constraints 

The stress constraints for this element are based on the Von Mises 
criterion and have the form 

g A = -j (®* + 3t*) - 1 i 0 ; i = 1.2. 3, 4 

* ( C-3 1 ) 

where a is the allowable stress. The normal stresses o . are obtained 
a i 

from technical beam theory and have the form 

M Y-X- cos0 + M Y.X. sin© +FI, ; i - 1,3 
y 1 o Z14 x 1 

-M Y,X a cos6 - M Y.X. sin0 +FL ; i = 2,4 (C-32) 

y 1 3 z 1 4 z 1 

The shear stresses are obtained from the snperposition of a parabolic 
horizontal shear stress distribution from technical beam theory and a 
uniform shear stress distribution due to the twisting moment from thin 
walled tube theory. The expressions for are given by 
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i = 1,3 


F Y?X, co 80 - F_ Y?X. sin0 + N I.X. 
y 1 3 z 1 4 zlZ 


T i = 


-F A, co 80 + F, Al sin0 + N Y 1 X. i 
y 1 o z l 4 xi* 


i = 2,4 


(C-33) 


Local Wall Bocklina Constraints 


To protect against local vail buckling of the thin tube element 
the following constraint between the wall thickness and mean radios is 
used (Ref , 58) . 


8 - * 1 i 0 


(C-34) 


with 


4 


1 .65x10 psi 


Y 1 


where a is the material yield stress in tension. 

y 


(C-35) 


Local Column Buckling Constraints 


To protect against column buckling of individual elements subject 
to axial compressive forces the following buckling stress constraint is 
used 


g 




-110 


(C-36) 


This interaction formula (Eq. (C-36)) is conservative (Ref. 59) when 
applied to elastically supported columns under combined loading. Equa- 
tion (C-34) is evaluated at two points on each end of the element as 
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shown in Fig. C7. The valnes of r and a are given by 

cr cr ° 


and 


a 

cr 




x > vr 



x l vT 



where 


FS 

X 

W 

a 

7 

E 

L 

k 


( 1.92 

( 1.67 + .265 X - .044 X 3 

.450 ft Jo /E 
1 V 7 


min 



561(L/2Y 1 ) * 
(2Y 1 /Y 2 )* 25 


51 


X > JT 
X < JT 


material yield strength in tension 
material modulus of elasticity 
column length 
effective length factor 


{ .0983 ; = 0 

.0629 ; 0 


The effective length factor k is dependent on the element end 


(C-37) 


(C-38) 


restraint 
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stiffnesses. Currently the value of k must be supplied as input data 
and is assumed to be constant during the design process. It should be 
noted that X and V are assumed to be constant during the solution of 
each approximate problem but they are updated for each approximate prob- 
lem. Additional details on the column buckling constraint formulation 
are available in Ref. 21. 
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Element: BOX BEAM 


Type : 15 

Description : This element is a frame element with the cross section 
described by the sizing variables B, H and t as shown in Fig. C8. Eight 
stress and four local buckling constraints are evaluated at each end of 
the element for a total of sixteen stress and eight local buckling con- 
straints. The locations on the cross section at which these constraints 
are evaluated are shown in Fig. C8 for the first node (n^) end of the 
element. 


Cross Sectional Dimensions 



( C-40 ) 


Reciprocal Section Properties 



1 

'lVWW 


* 2-1 


vw 

2<W 3 )2 


= _1_ 12 

3 Iy v5- ( V 2T 3 )(I r 2 v 3 


(C-41 ) 
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12 


X. 



Yi^-(ii-2Y 3 ) d 2 -n 3 ) 3 


Stress Constraints 

The stress constraints for this element are of the same form as 
those for element type 11 and they are given by Eqs. (C-5) - (C-7) with 
Y^ set equal to Yj. 

Local Buckling Constraints 

The local buckling constraints for this element are of the same 
form as those for element type 11 and they are given by Eqs. ( C— 8 ) - 
(C-14) with Y^ set equal to Y^. 
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APPENDIX D 


Derivatives of Structural Response Quantities 
with Respect to Element Reciprocal Section Properties 


D.l. Displacement Derivatives 


For the linear static structural analysis problem, the displace- 
ment derivatives are easily obtained through the implicit differentia- 
tion of the governing equilibrium equations with respect to the element 
reciprocal section properties (RSP's). In general, differentiation of 
Eq. (3-1) with respect to the j-th RSP of the i-th element (x^) yields 


dill 

3X ij 


(u) k + [K] 


3(h). 


dx 


ij 


d{P) k 

■^T ! 


k <= 1,2,. ..K 


(D-l) 


Under the assumption that the external loading is independent of the 

d{PL 


element RSP's (i.e. 


9l ij 


= 0) Eq. (D-l) becomes 


dtu} k _ 

IK] 77 — = V ijk ; k = 1 ,2 , . . .X 


(D-2) 


where the pseudo load vector ^ is given by 


dm 


ijk dx 


{»), 


ij 


(D-3) 


Writing the system stiffness matrix [K] as 


IK] = I [p.] T ([T i ] T [K i ] e [T i ])[p i ] 
i=l 


(D-4) 


where [K^] e is the element stiffness matrix in local coordinates, [T^] 
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is the element coordinate transformation matrix, [JK] is the element 
local to global degree of freedom transformation matrix and I is the 
total number of structural elements; substitution into Eq. (D-3) gives 

T T 

v Uk-“Wi ] IT i3 “aTT 

1J (D-5) 

Finally. Eq. (D-5) may be rewritten as 

V... = 1/x* ([p.J T [T ] T [K ] e [T ][p ] {u}.) 

J J J (D-6) 

▼here it is recognized that 

~ a 3 lK i ]e 

[ V - SuTTT 

J (D-7) 

is the unit element stiffness matrix formed by assigning the j-th sec- 
tion property a value of unity while the remaining section properties 
are set to zero. 

Using the expression for the pseudo load vector given by Eq. (D- 

6) , Eq. (D-2) can be solved for the unknown displacement derivatives 
d {u}^ 

“T ; k - 1,2,...K via the same procedure used to solve the equili- 

8i ij 

brium equations (Eq. (3-1)). Solving Eq. (D-2) directly yields the 
derivative values for all of the displacement degrees of freedom. For 
the case where the number of displacement degrees of freedom associated 
with the retained constraint set is fewer than the number of pseudo load 
vectors associated with Eq. (D-2) it is computationally more efficient 
to solve Eq. (D-2) using a partial inverse technique represented by the 
equation 
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(D-8) 


3 {u} k 

3x ij 


= [C]V 


ijk 


where {n}^ represents the displacement degrees of freedom associated 

with the retained constraint set. The partial inverse matrix [C] is 

T 

constructed such that its n-th row contains the vector {c} obtained 

n 

from the solution of the equation 


[KHc} = {e} 
n n 


(D-9) 


where {e}^ is a unit vector corresponding to the n-th degree of freedom 
associated with the retained constraint set. It should be noted that 
the solution of either Eq. (D-2) or Eq. (D-9) requires only the back 
substitution of the vectors or {®) n ^ the decomposed stiffness 
matrix has been saved from the previous structural analysis. 


D.2 Element Force Derivatives 


The element force derivatives are obtained through the implicit 
differentiation of the element force-displacement relations. Rewriting 
Eq. (3-8) in terms of the global coordinate system gives 


{F r J h = lK r ]8 { Vk + < FEF r } k 


(D-10) 


where {F )?, [K ] 8 , {u }. and {FEF )? are the element force vector, 
r x r r x r x 

stiffness matrix, vector of nodal displacements and fixed end force vec- 
tor for the r-th element and k-th load set. Differentiation of Eq. (D- 
10) with respect to x^ gives 
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3{FEF.}® 

“^T 

and Eq. (D-ll) becomes 


9l ij " 

Rewriting [K^]® as 

[K ] g = [T ] T [K ]®[T ] 
r r r r 

and substituting into Eq. (D-13) yields 



(D-12) 


i = r 


(D-13) 


(D-14) 
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T 3[K r ] T . 3UJ k 

'V -^V'-r’k * tT r ] T [K r ]'[T t ] — 


; i = r 


T e 3U r } k 

iv ( v [ v “^7 ; 1 * 1 


(D-15) 


Introducing the unit element stiffness matrix [K^]° # Eq. (D-15) becomes 


1 t ~ * *r - 

-f- ([T r ] T [K tj ] C lT r ]){u r ) k + [T r ] X [K r ] e [T r ] ; i*r 


8(F r ){ 


T e 3{u r } k 
[T r ] T « r ]»tT r ] -J— ; it r 

1J 


(D-16) 


Finally, writing the element force derivatives in the local coordinate 
system gives 


d Jl£ _ rT , •_*£ 

3x ij r 8x ij 


(D-17) 
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f- !E rjl‘ tT t H»r) k * [K r ]'[T t ] 

X ij 14 


3{u r } k 


atF r }® 


dx 


ij 


e 3 {u r } k 
[K ] e [T ] — 


r dx 


i i 4 r 


ij 


since, by orthogonality of [T^] , 


[T r ][T r ] A = IT r ][T r ] 1 = [I] 


i = r 


(D-18) 


(D-19) 
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APPENDIX E 


Data Command Descriptions 
E.JL Analysis Data Commands 

The synthesis problem structural analysis model may be defined 
using the analysis data commands listed below. 

* ANAL YSIS DATA * 

' BEAM ELEMENTS' 

' BOUN DARY CONDITIONS' 

' COOR DINATES' 

• LOAD CONDITIONS' 

' MATE RIALS' 

' NODAL LOADS' 

' TRAN SFORMATIONS' 

' TRUS S ELEMENTS' 

' UNIFO RM BEAM LOADING' 

These commands are described in this section in alphabetical order. 
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Command : * ANAL YSIS DATA * 

Description : This command denotes the beginning of the analysis input 
data block and must precede the first occurrence of any other analysis 
data command. 
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Command: 'BEAM ELEMENTS' 


Description : This command allows the user to define the beam element 
data for the structural analysis model. 


Data : n^ nl^ n2j [A J IYT IZZ IYZ] ^ Cj 


| 


n k nl k n2 fc [A J IYT IZZ IYZ] k 


n3 k 

°k 

(x,y,z) 


Entry 

Definition 

Note 

Variable 

n i 

Beam element number 


IDBEM(I,1) 

°i 

Material specification number 

(1) 

IBEAM(I,1) 


First end node number 

(2) 

IBEAM(I,2) 

^i 

Second end node number 

(2) 

IBEAM(I,3) 

A i 

Cross sectional area 

(3) 

PBEAM(I.l) 

J i 

Torsional constant 

(3) 

PBEAM(I,2) 

IYY. 

Moment of inertia about local 
y axis 

(3) 

PBEAM(I,3) 

izz A 

Moment of inertia about local 
z axis 

(3) 

PBEAM(I,4) 

IYZ . 
1 

Product of inertia 

(3) 

PBEAM(I,5) 

C i 

Beam orientation specification code 

(4) 

PBEAM(I,6) 

n3 i 

Beam orientation node number 

(4) 

PBEAM(I,7) 


Beam orientation angle 

(4) 

PBEAM(I,7) 
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FBEAM(I,7) 

Components of beam orientation vector (4) PBEAM(I,8) 

PBEAM(I,9) 

Total number of beam elements NBE 


Notes : 


1. The corresponding material must be defined by the ' MATER IAL ' com- 
mand • 


2. The end nodes must be defined by the 9 COOR DINATES 9 command, 

3, The beam section properties need not be given if the beam is to be 
associated with a design element. For the case where the beam is 
associated with a design element, any section properties supplied 
by the user will be replaced by section properties calculated from 
the corresponding design element information. 


The orientation of the beam element local y axis may be defined by 
supplying a node or vector which lies in the beam's local x-y 
plane (Fig. El) or an angle a (Fig. E2) which is a measure of the 
angle between the local x-y plane and the x i oca i ~ ^global P* ane 
(counter-clockwise looking down the local x-axis). For the spe- 
cial case where the local x-axis and global y-axis are coincident 
a is the measure of the angle between the local x-y plane and the 

X local ~ z global P* ane (Fig. E3) , The means by which the orienta- 
tion information is supplied is determined from c^ as follows: 
c^ = 1 - node point n3 ^ 
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c i = 2 - angle 

Cj = 3 - vector (x.y.zjj 
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Command; ’BOUNDARY CONDITIONS' 


Description : The boundary condition command allows the user to specify 
the nodal degrees of freedom which are to be considered restrained for 
the purposes of analysis. 

Data : ['SET' ] 




['SET* kj 


n„ SPC« 


\ SPC “l 
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Entry 

Definition 

Note 

Variable 

k. 

J 

Boundary condition set number 

(1) 

IDBCS( J, 1) 

% 

SPC 

j 

m. 

J 

Node number 

(2) 

ISPCS(I,1,J) 

Degree of freedom code 

(3) 

ISPCS(I,2.J) 

Number of restrained nodes 
for boundary condition set k^ 

(4) 

NCN 

L 

Number of boundary condition 
sets 


NBS 


Notes : 

1. The absence of the 'SET' descriptor and its associated boundary 
condition number will cause the following boundary condition 
specifications to be applied to all boundary condition sets, 

2. The use of the ' ALL ' descriptor .in place of the node number will 
cause the associated boundary condition specification to be 
applied to all nodes, 

3. Any combination of the values 1-6 may be used to specify res- 
traints on the nodal degrees of freedom. The values 1,2 ,3 and 
4,5,6 correspond to the x,y,z translations and rotations, respec- 
tively. A negative sign preceeding any individual value will 
cause a previously set restraint to be removed, 

4. The maximum number of restrained nodes for all boundary condition 

sets is stored in the variable NCN (i,e, NCN = max m.) 

j-l.L. J 
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Command: * COORDINATES ' 


Description : This command allows the user to specify the problem node 
point data. 

Data.: n x 1 y 1 [k 1 ] 


n x y z [k ] 
m m ■'m m m 


Entry 

Definition 

Note 

Variable 

n. 

l 

Node number 


IDGRD(I.l) 

x i 

X-coordinate 


PGRID (1,1) 

y i 

Y-coordinate 


PGRID(I,2) 

z i 

Z-coordinate 


PGRID (1,3) 

k. 

1 

User coordinate system number 

(i) 

IGRID(I) 

m 

Number of nodes 


NND 


Notes : 

1, The user coordinate system number k^ refers to a user supplied 
coordinate transformation given by the 'TRANSFORMATION' command. 
All quantities associated with the node except its location are 
described in the user coordinate system. If the transformation 
specification includes an origin specification then the node loca- 
tion is also described in the user coordinate system. 
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Command : ' LOAD CONDITIONS' 

Description : The load condition command allows the user to specify the 
load sets and their corresponding analysis types and boundary condition 
numbers. 

Data : n^ 'type'^ [bc^] 


n 'type' [be ] 
m m m 


Entry 

Definition 

Note 

Variable 

n i 

Load set number 


IDL0D(I,2) 

' type ' i 

Analysis type 

(I) 

LOADT (1,1) 

b°i 

Boundary condition number 

(2) 

LOADT(I,2) 

m 

Number of load sets 


NLS 


Notes ; 

1. The descriptor 'type' defines the type of analysis to be performed 
for the associated load set. Currently only linear static 
analysis may be performed and therefore 1 STAT ICS 9 is the only 
valid descriptor. 

2. The absence of the boundary condition number bc^ will cause the 
load set to be associated with boundary condition number 1. 
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Command : ’MATERIALS* 


Description : This command allovs the user to specify material properties 
to be used for analysis and design pnrposes. 

Data : Ej Gj t ^ [C^] 


« P« E „ V n G „ ° ■n T „ l C ,J 
n nnnnnn n 


Entry 

Definition 

Note 

Variable 

m i 

Material number 


IDMAT (1,1) 

p i 

Mass density 


PMATE(I,1) 

E i 

Modulus of elasticity 


PMATE(I,2) 

V. 

1 

Poissons ratio 


PMATE(I,3) 

G i 

Shear modulus 

(1) 

PMATE(I,4) 

*i 

Allowable stress in tension 
and compression 


PMATEU.8) 

*i 

Allowable stress in shear 


PMATE(I,9) 

c. 

1 

Additional material constant 

(2) 

PMATE(I.IO) 

n 

Total number of materials 


NMT 

Notes : 





1. If a valne of 0. is supplied for G the actual value of G will be 
calculated from 


2(l+v) 
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2. The additional material constant may be used to supply any 
additional data reqnired for the element strength constraint cal- 
culations (e.g. effective member length for buckling constraint 
evaluation, factor of safety) . 
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Command: 'NODAL LOADS' 


Description : This command allows the nser to describe discrete loads to 
be applied to specific nodal degrees of freedom. 


Data: 


['SET' k x ] 


n DOFj P 1 


n ml D0F ml P ml 


['SET' kj^] 


n mj+l D0F mj+l P mj+1 


n mL D0F mL P mL 


Entry 


Definition 


Note Variable 


“j 

DOF. 


P. 

J 


tnT . 


Load set number (1) 

Node number 

Degree of freedom specification (2) 

Magnitude of load 

Number of load sets 

Number of applied nodal loads 


NLOAD(J.l) 

NL0AD(J,2) 

NLOAD(J.3) 

FLOAD(J) 

NLS 

NAF 
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Notes: 


1. The absence of the ' SET 1 descriptor and its associated load set 
number will cause the following loading data to be associated with 
load set number 1. The load set number must correspond to a load 
set defined by the ' LOAD CONDITIONS' command. 

2. The degree of freedom specification may be any value between 1 and 
6. The values 1, 2, 3 and 4, 5* 6 correspond to the x,y, z trans- 
lations and rotations respectively. 
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Command: ‘TRANSFORMATIONS'* 


Description : This command may be used to to create user defined coordi- 
nate systems to facilitate data input and/or to impose boundary condi- 
tions and displacement constraints in directions other than the global 
coordinate directions. 

Data : n^ * axisl xl^ yl^ zl^ 'axis2'^ x2^ y2^ z2^ [*0* xo^ yo^ zo^] 


n 

m 


'axisl' 

m 


xl m yl m zl m 


axis2 ' n x2 - y 2 ™ z2 » i’ 0 ' x0 ™ y°» 


m 


Entry 


Definition 


Note Variable 


*i 

’axisl' 


xl. 

1 

yl i 


zl 


i 

axis2 ' 


x2 


y2 


i 

z2 i ) 


Transformation number 
First user axis designation 

Components of user 'axisl* 

Second user axis designation 

Components of user *axis2’ 

Indicates that user coordinate 
system origin specification follows 


( 1 ) 

( 2 ) 

( 1 ) 

( 2 ) 

(3) 


IDTRN(I,1) 

CTRAN(l.I) 
CTRAN(2 , I) 
CTRAN(3 , I) 

CTRAN(4,I) 

CTRAN(5,I) 

CTRAN(6,I) 
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Location of user coordinate 
system origin 


CTRAN(10,I) 



Notes; 


(3) CTRAN(11,I) 
CTRAN(12,I) 


1. The user axis designations 'axisl' and 'axis2' may be given as 

T» or 'Z' and need not be given in cyclic permutation 

order. 

2. The components of the nser coordinate axes mnst be given in terms 
of the global reference system. The axes 'axisl' and 'axis2' need 
not be perpendicular. The third user axis ('axis3') is calculated 
from the cross product of 'axisl' and *axis2' and then 'axis2' is 
recomputed from the cross product of 'axisl' and 'axis3'. 

3. If the origin specification 'O' is given then the user coordinate 
system origin (xo,yo,zo) must be given in terms of the global 
reference system. In this case the node point locations supplied 
via the ' COORD INATES ' command are assumed to be measured with 
respect to the user coordinate system origin. 


* The capabilities available through the use of this command are still 
under development at this time and have not been fully tested. 
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Command : * TRUSS ELEMENTS' 


Description : This command allows the user to define the truss element 
data for the structural analysis model. 

Data : n^ nlj n2j [A]j 

a 

nl k n2 k [A]^ 


Entry 

Definition 

Note 

Variable 

°i 

Truss element number 


IDTRS(I.l) 

"i 

Material specification number 

(1) 

ITRUS(I.l) 

nl . 

i 

First end node number 

(2) 

ITRUS(I,2) 

n2i 

Second end node number 

(2) 

ITRUS(I,3) 

A i 

Cross sectional area 

(3) 

PTRUS(I,1) 

k 

Total number of truss elements 


NTE 


Notes : 

1. The corresponding material must be defined by the 9 MATE RIAL 9 com- 
mand. 

2. The end nodes must be defined by the 9 COOR DINATES 9 command. 

3. The truss cross sectional area need not be given if the truss is 
to be associated with a design element. For the case where the 
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trass is associated with a design element any cross sectional area 
supplied by the user will be replaced by the cross sectional area 
calculated from the corresponding design element information. 
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Command: ' UNIFORM BEAM LOADING' 


Description : This command allows the user to specify uniformly distri- 
buted loading for any beam element in the structural model. This load- 
ing is applied in the form of work equivalent nodal loading. 

Data : ['SET' kj] 

n DIRj Pj 


“ml DIR ml P ml 


['SET' 1^] 

n .... DIR P . , . 
oj+1 mj+1 mj+1 


“mL WR mL P mL 
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Entry 

Definition 

Note 

Variable 

k i 

Load set number 

(1) 

IUL0D(J,1) 

“j 

Beam element number 


IUL0D(J,2) 

DiRj 

Direction specification 

(2) 

IULOD( J,3) 

P. 

J 

Magnitude of loading in 
load per unit length 


DLOAD(J) 

L 

Number of load sets 


NLS 

mT . 

Number of applied uniform 
beam loads 


NULB 


Notes: 


1. The absence of the ' SET * descriptor and its associated load set 
number will cause the following loading data to be associated with 
load set number 1. The load set number must correspond to a load 
set defined by the 1 LOAD CONDITIONS' command. 

2. The direction specification defines the direction of the applied 
loading in the beam element local coordinate system. The specifi- 
cation may be any value between 1 and 6. The values 1 , 2,3 and 
4,5,6 correspond to the directions of local z,y, z translations and 
rotations, respectively. 


175 



E.2 Design Data Commands 


The synthesis problem design model may be defined nsing the design 
data commands listed below. 

DESIG N DATA * 

* DESIG N ELEMENTS' 

' DISPL ACEMENT CONSTRAINTS' 

' ELEM ENT GEOMETRY' 

' LOCAL BUCKLING CONSTRAINTS' 

' STRE SS CONSTRAINTS' 

These commands are described in this section in alphabetical order. 
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Command : ' DESIG N DATA ' 

Description : This command denotes the beginning of design inpnt 
block and most precede the first occurrence of any other design 
command. 


data 

data 
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Command: * DESIGN ELEMENTS' 


Description : This command allows the user to specify the structural ele- 
ments which are to be considered as design elements during the design 
process. 

Data : ml^ [m2^] 


ml^ n k [m2 k ] 


Entry 

Definition 

Note 

Variable 

m!i 

Analysis element number 


IDESG(I.l) 

n i 

Geometry number 

(1) 

IDESG (1,2) 

m2 ^ 

Master element number 

(2) 

IDESG(I,4) 

k 

Number of design elements 


NDE 


Notes: 


1. The geometry number n^ must correspond to a design element 
geometry number specified by the * ELEM ENT GEOMETRY* command. 

2. The master element number m2^ allows the user to link several 
analysis elements together for design purposes. The master ele- 
ment number must appear in the design element list and must not 
itself refer to another master element. The absence of a user 
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specified master element causes m2^ to be set equal to ml^. 
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Command; DISPLACEMENT CONSTRAINTS' 


Description : This command allows the user to specify allowable upper and 
lower bounds on nodal translations and rotations. 

Data : Uj jj Cj LBj 


n j c LB^ DB^ 


Entry 

Definition 

Note 

Variable 

m i 

Load set number 


IDISC(I.l) 

n i 

Node number 


JDISC(I,2) 

h 

Degree of freedom number 

(1) 

IDISC(I,3) 

c i 

Constraint shift value 

(2) 

RDISC(I ,1) 

“i 

Lower bound value 

(2) 

RDISC(I,2) 

ra i 

Upper bound value 

(2) 

RDISC(I,3) 

k 

Number of displacement constraints 


NDC 

Notes : 

1 . 

The degree of freedom specification may 

be 

any value between 


6. The values 1, 2, 3, and 4, 5, 6 correspond to the x, y, z 

translations and rotations, respectively. 

2. The displacement constraint is written for strictly negative lower 
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bound and strictly positive upper bound values. For the case 
where the actual displacement bounds are of the same sign the con- 
straint shift value c^ may be used to shift the constraint value 
such that the upper and lower bounds will have the correct signs. 
The form of the displacement constraint is given by 

(u-c)/u a £ 1 

where u is the displacement value and u is the shifted displace- 

ft 

ment allowable (i.e. LB ^ or FB^) . 
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Command: * ELEMENT GEOMETRY' 


Description : This command allows the user to specify data pertaining to 
the design element sizing variables. 


Data: 

"i 

"1 





Y 1 

• 

TO 1 


• 

* 

Y kl 

• 

iu tl 


• 

"l 

n L 





Y ki+1 

• 

TO kl*l 



Y kL 

• 

^kL 

™kL 

Entry 


Definition 

Note 


m i 

Design element geometry number 


IDG0M(I,1) 

a i 

Design element type 

(i) 

IGEOM(I.l) 

y j 

Initial value of design element 
sizing variable 


DGE0M(J,1) 


Lower bound value of design 
element sizing variable 


DGE0M(J,2) 

"j 

Upper bound value of design 
element sizing variable 


DGE0M(J,3) 

L 

Number of design element 
geometries 


NEG 
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kL 

Notes 

1 . 


Total number of design element NCD 

sizing variables 


The design element type must correspond to an element type avail- 
able in the design element library. 
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Command: 'LOCAL BUCKLING CONSTRAINTS' 


Description : This command allovs the nser to request the calculation of 
local buckling constraint values for the specified load sets. Local 
buckling constraint values are calculated for all elements whose 
corresponding element type has an associated local buckling computation 
as described in the element library. 

Data : k^ 


k 

n 

Entry Definition Note Variable 

k^ Load set number IBLOD(I) 

n Number of load sets for which NBL 

buckling constraints are to be calculated 
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Command: * STRESS CONSTRAINTS' 


Description : This command allows the user to request the calculation of 
design element stress constraint values for the specified load sets. 
Stress constraint values are calculated for all elements whose 
corresponding element type has an associated stress computation as 
described in the element library. 


Data: 

k i 

• 



• 

• 

k 



n 


Entrv 

Definition 

Note Variable 

k. 

X 

Load set number 

ISLOD(I) 

n 

Number of load sets for vhich stress 
constraints are to be calculated 

NSL 
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E. 3 Control Data Commands 


The solution of the structural synthesis problem may be controlled 
using the commands listed below. 

' ANAL YSIS* 

' CHEC KPOINT* 

' CONM IN' 

* CONT ROL DATA' 

* CSD ' 

' DUAL * 

' FORC E VARIANCE' 

' ITER ATIONS' 

' MIXE D* 

' MOVE LIMITS' 

' OPTI MIZATION' 

' PRIN T' 

' REST ART* 

' SCAL E' 

' SENS ITIVITIES' 
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' SETUP' 


'UPDATE' 


These commands are described in this section in alphabetical order. 
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Command: 'ANALYSIS' 


Description : This program function control command causes the program to 
terminate after the structural analysis has been completed. In this 
case the program control variable IPCTL is set to 1. 
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Command: ' CHECK POINT' [n] 


Description : This command causes the program analysis and design data to 
be written on external file nnmber n upon successful termination of a 
design run. The external file number is stored in the variable ICKFL. 
The absence of the file number specification will cause ICKFL to be set 
to 1 . 
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Command: 'CONMIN' [c] 


Description : This command allows the nser to specify that the solution 
to the approximate design problem is to be performed using a primal 
mathematical programming formulation. The mass minimization problem is 
solved using the CONMIN optimization program. The parameter c controls 
the constraint push off factor. If c > 0. then the displacement, stress 
and local buckling constraints are treated as being nonlinear and c is 
used as the push off factor. If c «* 0. then all constraints are con- 
sidered to be linear. In the absence of a user specified value c is 
given a value of 1.0. Specification of the 'CONMIN' command causes the 
variable IOPTY to be set to 3. The value of c is stored in the variable 
OPTPRM(l) . 
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Command: 'CONTROL DATA' 


Description : This command denotes the beginning of the program control 
data block and must precede the first occurrence of any other control 
data command. 
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Command: 'CSD 


* 


Description : This command allows the user to specify that the element 
cross sectional dimensions are to be used as the design variables during 
the solution of each approximate design problem. In the absence of this 
command either the element reciprocal section properties or a combina- 
tion of element reciprocal section properties and cross sectional dimen- 
sions are chosen as the design variables. Specification of this command 
causes the variable ICSD to be set to 1. 
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Command : ’ DUAL * 


Description : This command allows the user to specify that the solution 
to the approximate design problem is to be performed using a dual 
mathematical programming formulation. The dual function maximization 
problem is solved using the CONMIN optimization program. This feature 
is currently operational only when used in conjunction with the ' CSD ' 
command and will automatically specify that a mixed variable approxima- 
tion is to be used for all behavior constraints. Design variable scal- 
ing is also activated automatically for this feature. Specification of 
this command causes the variables IOPTY, IMIX, and ISCALE to be set to 
2, 1 and 1, respectively. 
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Command: 'FORCE VARIANCE' n 


Description : This command allows the users to specify that the variation 
of the design element forces with respect to the problem design vari- 
ables is to be included in the stress and local buckling constraint 
approximations during the approximate problem generation. If this com- 
mand is not specified the element forces are assumed to be invariant 
during any design step. If n = 1 then the element force sensitivities 
are calculated only with respect to the design variables associated with 
that element. If n = 2 then the element force sensitivities are calcu- 
lated with respect to all design variables. The value of n is stored in 
the variable IFVAR. 
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Command : ' ITER ATIONS ' n 


Description : This command allows the nser to specify the number of 
design steps which are to be performed. Each design step consists of 
the generation and solution of one approximate problem. The number of 


iterations n is stored in the variable NSTEP 



Command : * MIXE D APPROXIMATIONS' 

Description : This command allows the user to specify that a mixed vari- 
able approximation will be used for all behavior constraints. Design 
variable scaling is activated automatically for this feature. Specifi- 
cation of this command causes the variables IHIX and ISCALE to be set to 
1 . 
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Command : ' MOVE LIMITS' d ± & 2 


Description : This command allows the user to specify the allowable 
changes in the problem variables during any single design step. The 
move limits d^ and dj are applied to the element cross sectional dimen- 
sions and reciprocal section properties, respectively. The limiting 
values (Y 1 , Y U , X L , X U ) are given by the following equations: 

Y L = Y - Yd x 

Y D = Y + Yd 1 

X L = X - Xd 2 

X U = X + Xd 2 

If zero is given for either move limit then the limiting values for the 
corresponding variables are determined from the overall limits on the 
element sizing variables supplied via the * GEOME TRY command. The values 
of d^ and d 2 are stored in the variables DMOVE(l) and DM0VE(2) , respec- 
tively. 
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Command: * OPTI MIZATION * [c] 


Description : This program function control command causes all major pro- 
gram functions (data processing, analysis, approximate problem genera- 
tion, optimization and design recovery) to be performed. In this case 
the program control variable IPCTL is set to 3. The value c controls 
the diminishing returns convergence criterion on structural mass. The 
design process is terminated if the relative change in structural mass 
is less than c for three consecutive design steps. The value of c is 
stored in the variable DELOBJ. In the absence of a user specified value 
for c DELOBJ is set to .01. 
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Command : ' PRIN T* n^ n^ n^ n^ ng n^^ 

Description : This command allows the user to control the program print- 
ing options. If the value of n^ is set to 0 then the corresponding 
printing option is disabled. A value of n^ greater than 0 enables the 
printing options as described below. The value of n^ is stored in the 
variable IPS(I). 

Design data is printed at the beginning of each 
design step. 

Analysis data is printed at the beginning of each 
design step. 

Not used at this time. 

Structural analysis results are printed for each 
design step. 

Design results are printed for each design step. 

CPU timing summary is printed for each design step. 

Not used at this time. 

Not used at this time. 

All program data storage is printed at completion 
of specified program function. 

Controls optimizer print option. 


♦ Currently n^ is used to control the CONMIN optimizer output. A value 
of n^g = 1 is suggested for normal program operation with higher values 
(up to = 6) being used for debug operations. 
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Command ; » REST ART* [n] 


Description : This command allows the nser to begin the design process 
from the termination point of a previous design run by causing the pre- 
viously written analysis and design data to be read from external file 
number n. The external file number is stored in the variable IRSFL. 
The absence of the external file number specification will cause IRSFL 
to be set to 1. It should be noted that at this time only control data 
modifications are allowable for restarted design runs. 
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Command : * SCAL E* 


Description : This command will cause the design variables to be scaled 
to unity at the beginning of each approximate problem stage. Specifica- 
tion of this command causes the variable ISCALE to be set to 1. The 
default value of ISCALE is 0 and scaling is not performed. 
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Command: ' SENSITIVITIES ' 


Description : This program function control command causes the program to 
terminate after the approximate problem generation has been completed. 
In this case the program control variable IPCTL is set to 2. 
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Command : * SETUP 9 


Description : This program function control command causes the program to 
terminate after the input data processing has been completed. In this 
case the program control variable IPCTL is set to the default value 0. 
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Command : * UPDA TE ' [n] 


Description : This command will canse the approximate design problem to 
be partially reconstructed and re-solved n times between each complete 
structural analysis and approximate problem generation. The value of n 
is stored in the variable NUPDAT. If n is not specified NUPDAT is set 
to 1. In the absence of the 'UPDATE* command NUPDAT is set to 0. 
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Fig. 1 - Structural Element Orientation 
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Fig. 2 


“fF 
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D.O.F.’s independent of changing 
boundary conditions 


D.O.F.'s associated with changing 
boundary conditions 


- Conceptual Matrix Partitioning for Multiple Boundary Conditions 
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Fig. 3 - Approximate Problem Options in RSP Design Space 
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Fig, 4 - Approximate Problem Options in CSD Design Space 
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Fig. 5 - Available Optimization Options 
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Fig. 6 - Program Organization 























lommand 


INITIALIZE 

Initialize data management 
dictionary and void area table. 

INSERT 

Insert new variable in the data vector 

LOCATE 

Locate a variable in the data vector. 

DELETE 

Delete a variable from the data vector 

EXPAND 

Increase the storage available for 
a variable. 

CONTRACT 

Decrease the storage avaiable for 
a variable. 

COMPRESS 

Remove voids from data vector. 

CHANGE NAME 

Change the name of a variable 
in storage. 

DEBDG 

Print data vector, dictionary and void 
area table. 

QUERY 

Return maximum storage used. 

CHECKPOINT 

Write data vector, dictionary and 
void area table on external file. 

RESTART 

Read data vector, dictionary and void 
area table from external file. 


Fig. 8. Storage Management Commands 
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MAIN — f- CONTRL 


PSTPRO 
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Fig. 9 - Main Routine Flow Diagram 
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-READA 
-READD 
— READC 
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SCARD 

SCARD 
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SORTI 
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SRCHI 

SHFTI 
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VCPRD 

VSPRD 


H VSPRD 
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CSDPRO 
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U RESTRT 


Fig. 10 - Pre-processor Flow Diagram 
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ANCTL 


APPCTL 


OPTCTL 


RLZCTL 


CONVRG 


PRTCTL 


Fig. 11 - Design Control Flow Diagram 
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r-BLDLBN 


I— BLDDFT 
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{ ROTATE H 
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r 
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TRSTRN 

ASSEMB 


MTMLT1 


MTMLT1 
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Fig. 12 - Analysis Control Flow Diagram 
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| — SKYFCT 


SKYFAC 



VIPSS 
SKY RDM 
SKYSOL 



VIPSS 

SKYMUL 
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— TRNMAT 
— INTBIN 
— BEMFEM 
— LODTRN 


MTMLTl 


H SKYSOL 



VIPSS 

SKYMUL 


•— DISPL INTBIN 
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“TRNMAT MTMLTl 
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MTMLTl 


Fig. 12 - Analysis Control Flow Diagram (cont.). 
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pAPGDAT 

-CONUPD — UPDATE 
-OBJCXY 


I— DISPLB 


l-CONEVL 
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STRS12 

STRS13 

STRS14 


- LBUK11 
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L INTBIN 
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Fig. 13 - Approximate Problem Generation Control Flow Diagram 
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h- DFORBM 


H SENSIT H 
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-BEMTRN 
-MTMLT4 
-MTMLT5 

- RDSPGR 
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Fig. 13 - Approximate Problem Generation Control Flow Diagram (cont.) 
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Fig. 13 - Approximate Problem Generation Control Flow Diagram 


(cont . ) 
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-SCALDV 

-SCALOB 

-SCALGR 


HDULINT H 
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Fig. 14 - Optimization Control Flow Diagram 
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Fig. 15 - Design Recovery Control Flow Diagram 
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rDDPRNT 
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Fig. 16 - Printing Control Flow Diagram 
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rCLENUP 


PSTPRO H 


‘-CHECKP 


17 - Post-processor Flow Diagram 


COMMAND' 


data [data] 


Form 2 

COMMA ND* 


data 


data 


Fora 3. 


COMMAND* 


* SUB- COMMAND* 


data 


data 


[data] 


[data] 


[data] 


[data] 


Fig. 18 - Data Command Forms 



</> <n -c/> — <r> <s> <s> <r> -</> <r> -c/> </> <r> </> <s> -</> <r> </> </> - <r> -<s> <r> </> <r> <s> <s> <s> 




I — u 




FRAME OPTIMIZATION PROBLEM - ONE BAY / TWO STORY FRAME 


DESIGN DATA BLOCK 
DESIGN DATA' 

DESIGN ELEMENTS' 

ANALYSIS ELEMENT NO. 


1 

2 

3 

5 
8 
7 

4 

6 


GEOMETRY NO. 

1 

1 

1 

1 

1 

1 

1 

1 


MASTER ELEMENT NO. 


ELEMENT GEOMETRIES' 

GEOMETRY NO. GEOMETRY TYPE 

1 14 

INITIAL VALUE LOWER BOUND 

16.61 1.00 

.45 .01 

STRESS CONSTRAINTS' 

LOAD SET NO. 

1 

2 

LOCAL BUCKLING CONSTRAINTS’ 

LOAD SET NO. 

1 

2 

$ 

'DISPLACEMENT CONSTRAINTS' 

$ 


1 

2 

3 

5 


UPPER BOUND 

25. 

5. 


Fig. 19 - Sample Program Input Data 
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</> </> -c /> --c/> -t/> </> -c/> - ~c/> -c/> </> </> — </> - <r> </> <s> 


LOAD SET NO. 

NODE NO. 

DOF 

SHIFT 

LOVER BOUND 

UPPER BOUND 

2 

2 

1 

0.0 

-.36 

.36 

2 

3 

1 

0.0 

-.72 

.72 

2 

4 

1 

0.0 

-.36 

.36 

2 

5 

1 

0.0 

-.72 

.72 

2 

6 

1 

0.0 

-.36 

.36 

2 

7 

1 

0.0 

-.72 

.72 


ANALYSIS DATA BLOCK 


ANALYSIS DATA 


f 


BOUNDARY CONDITIONS' 

NODE NO. DEGREE OF 

'ALL' 

2 

3 

4 

5 

6 
7 


FREEDOM SPECIFICATION 

123456 

-126 

-126 

-126 

-126 

-126 

-126 


MATERIAL NO. 

RHO 

E 

NU 

G 

SIGY 


FS 

K-EFF 

1 

.2836 

30.00E6 

0.30 

0.785E7 

3.6E04 


1.51 

1.00 

2 

.2836 

30.00E6 

0.30 

0.785E7 

3.6E04 


1.51 

1.14 

3 

.2836 

30.00E6 

0.30 

0.785E7 

3.6E04 


1.51 

1.38 

4 

.2836 

30.00E6 

0.30 

0.785E7 

3.6E04 


1.51 

1.37 

BEAMS ' 









ELEMENT NO. 

MATERIAL NO. I 

NODE 

1 NODE 

2 ORIENTATION DATA 

1 

1 


1 

2 


3 

1 . 

0 . 0 . 

2 

2 


2 

3 


3 

1 . 

0 . 0 . 

3 

3 


3 

5 


3 

0 . 

1 . 0 . 

4 

3 


5 

7 


3 

0 . 

1 . 0 . 

5 

4 


2 

4 


3 

0 . 

1 . 0 . 

6 

4 


4 

6 


3 

0 . 

1 . 0 . 

7 

2 


6 

7 


3 

1 . 

0 . 0 . 

8 

1 


6 

8 


3 

1 . 

0 . 0 . 


$ 

'COORDINATES’ 

$ 


Fig. 19 - Sample Program Input Data (cont.) 


227 


<r> <r> </> <n <r> <r> </> <n <r> - </> </> </> </> </> <r> </> - -c/> -c/> -c/> </> 


NODE NO. 

X 

Y 

Z 

1 

0.0 

0.0 

0.0 

2 

0.0 

180.0 

0.0 

3 

0.0 

360.0 

0.0 

4 

120.0 

180.0 

0.0 

5 

120.0 

360.0 

0.0 

6 

240.0 

180.0 

0.0 

7 

240.0 

360.0 

0.0 

8 

240.0 

0.0 

0.0 


$ 

'LOAD CONDITIONS' 


LOAD SET NO. 

TYPE 

BOUNDARY CONDITION NO 

1 

'STATICS' 

1 

2 

'STATICS' 

1 


NODAL LOADS 


I 


LOAD SET NO. 

'SET' 2 

NODE NO. 

2 
3 

UNIFORM LOADS' 


DIRECTION 

1 

1 


MAGNITUDE 

45000. 

45000. 


LOAD SET NO. 



'SET' 1 

BEAM NO. 

DIRECTION 

MAGNITUDE 

3 

2 

-500. 

4 

2 

-500. 

5 

2 

-500. 

6 

2 

-500. 


CONTROL DATA BLOCK 


'CONTROL DATA' 

? 

'ITERATIONS' 15 
'MOVE LIMIT' .0 .5 

'CONMIN' 


Fig. 19 - Sample Program Input Data (cont.) 


228 



</></>*/> 


'MIXED' 

' SCALE ' 

'OPTIMIZATION' .001 
'PRINT' 0 0 0 1 1 1 0 0 0 1 


END OF DATA 


Fig. 19 - Sample Program Input Data (cont.) 
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WEIGHT (lb) 


800.0 


750.0 


700.0 


650.0 


600.0 


550.0 



500.0 


450.0 I I I I I I I I 

0 1 2 3 4 5 6 7 

ANALYSIS NUMBER 


Fig. 21 - Iteration History for Problem 1, Run 1 (Option 1(P)) 
Tied Cantilevered Beam 
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WEIGHT (lb) 



ANALYSIS NUMBER 


Fig. 24 - Iteration History for Problem 1, Run 4 (Option 6(P)) 
Tied Cantilevered Beam 
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MASS (kg) 



Fig. 26 - Iteration History for Problem 2, Case A, Run 1 (Option 1(P)) 
Two Member Frame 
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MASS (kg) 



Fig. 27 - Iteration History for Problem 2, Case A, Run 2 (Option 2(P)) 
Two Member Frame 
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MASS (kg) 


O CONSTRAINT VIOLATION < 1% 
□ CONSTRAINT VIOLATION < 10% 
• CONSTRAINT VIOLATION > 10% 



Fig. 28 - Iteration History for Problem 2, Case A, Run 3 (Option 3(P)) 
Two Member Frame 


MASS (kg) 



0 5 10 15 

ANALYSIS NUMBER 


Fig. 29 - Iteration History for Problem 2, Case A, Run 4 (Option 4(P)) 
Two Member Frame 


MASS (kg) 



Fig. 30 - Iteration History for Problem 2, Case A, Run 5 (Option 7(P)) 
Two Member Frame 
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MASS (kg) 



Fig, 31 - Iteration History for Problem 2, Case A, Run 6 (Option 10(P)) 
Two Member Frame 
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MASS (kg) 



Fig. 32 - Iteration History for Problem 2, Case A, Run 7 (Option 10(D)) 
Two Member Frame 
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MASS (kg) 



Fig. 33 - Iteration History for Problem 2, Case A, Run 8 (Option 1(PU)) 
Two Member Frame 


243 


MASS (kg) 


800.0 


700.0 


600.0 


500.0 


400.0 


300.0 


200.0 


100.0 



Fig. 34 - Iteration History for Problem 2, Case B, Run 1 (Option HD) 
Two Member Frame 


2 


MASS (kg) 



Fig. 35 - Iteration History for Problem 2, Case B, Run 2 (Option 2(F)) 
Two Member Frame 
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MASS (kg) 



Fig. 36 - Iteration History for Problem 2, Case B, Run 3 (Option 3(P)) 
Two Member Frame 
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MASS (kg) 



Fig. 37 - Iteration History for Problem 2, Case B, Run 4 (Option 4(P)) 
Two Member Frame 
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MASS (kg) 



ANALYSIS NUMBER 


Fig. 38 - Iteration History for Problem 2, Case B, Run 5 (Option 10(P)) 
Two Member Frame 
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MASS (kg) 



ANALYSIS NUMBER 


Fig. 39 - Iteration History for Problem 2, Case B, Run 6 (Option 10(D)) 
Two Member Frame 
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MASS (kg) 



ANALYSIS NUMBER 


Fig. AO - Iteration History for Problem 2, Case B, Run 7 .(Option 1(PU)) 
Two Member Frame 
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VOLUME (in' 5 ) 



ANALYSIS NUMBER 


Fig. 42 - Iteration History for Problem 3, Run 1 (Option 1(P)) 
Three Member Frame 
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VOLUME (in J ) 



Fig. A3 - Iteration History for Problem 3, Run 2 (Option 2(P)) 
Three Member Frame 
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VOLUME (in' 5 ) 



Fig, 44 - Iteration History for Problem 3, Run 3 (Option 3(P)) 
Three Member Frame 
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VOLUME (in' 5 ) 



Fig. 45 - Iteration History for Problem 3, Run 4 (Option 4(P)) 
Three Member Frame 
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VOLUME (in J ) 



Fig. 46 - Iteration History for Problem 3, Run 5 (Option 10(P)) 
Three Member Frame 
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VOLUME (in 3 ! 



Fig. 47 - Iteration History for Problem 3, Run 6 (Option 10(D)) 
Three Member Frame 
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VOLUME (in'*) 



Fig* 48 - Iteration History for Problem 3, Run 7 (Option 1(PU)) 
Three Member Frame 
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MAbb (kg) 


O CONSTRAINT VIOLATION < 1% 
□ CONSTRAINT VIOLATION < 10% 
• CONSTRAINT VIOLATION > 10% 


ANALYSIS NUMBER 


Fig, 51 - Iteration History for Problem 4, Run 3 (Option 3(P)) 
Seven Member Frame 













MASS (kg) 



Fig. 55 - Iteration History for Problem 4, Run 7 (Option 3(PU)) 
Seven Member Frame 



VOLUME (cm J ) 



Fig. 57 - Iteration History for Problem 5, Run 1 (Option 1(P)) 
Portal Frame 
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O CONSTRAINT VIOLATION < 1% 
□ CONSTRAINT VIOLATION < 10% 
• CONSTRAINT VIOLATION > 10% 





VOLUME (cm 



Fig. 60 - Iteration History for Problem 5, Run 4 (Option 12(P)) 
Portal Frame 
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VOLUME (cm J ) 



Pig. 61 - Iteration History for Problem 5, Run 5 (Option 12(D)) 
Portal Frame 
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500 lb/in (L.C. 1) 



240" 


Fig. 62 - One Bay / Two Story Frame (Problem 6) 


2 




WEIGHT (lb) 



ANALYSIS NUMBER 


Fig, 63 - Iteration History for Problem 6, Case A, Run 1 (Option 1(F)) 
One Bay / Two Story Frame 



ANALYSIS NUMBER 


Iteration History for Problem 6, Case A, Hnn 2 (Option 3(P)) 
One Bay / Two Story Frame 


WEIGHT (lb) 



ANALYSIS NUMBER 


Fig. 65 - Iteration History for Problem 6, Case A, Run 3 (Option 6(P)) 
One Bay / Two Stroy Frame 


WEIGHT (lb) 


Mg 


O CONSTRAINT VIOLATION < 1% 
□ CONSTRAINT VIOLATION < 10% 
• CONSTRAINT VIOLATION > 10% 



ANALYSIS NUMBER 


[•!*] 


Iteration History for Problem 6, Case A, Run 4 (Option 12(P)) 
One Bay / Two Story Frame 





WEIGHT (lb) 



Fig. 67 - Iteration History for Problem 6, Case A, Run 5 (Option 12(D)) 
One Bay / Two Story Frame 
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WEIGHT (lb) 



ANALYSIS NUMBER 


Fig. 68 - Iteration History for Problem 6, Case B, Run 1 (Option 1(P)) 
One Bay / Two Story Frame 
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WEIGHT (lb) 



ANALYSIS NUMBER 


Fig. 69 - Iteration History for Problem 6, Case B, Run 2 (Option 3(P)) 
One Bay / Two Story Frame 
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WEIGHT (lb) 



Fig, 72 - Iteration History for Problem 6, Case B, Run 5 (Option 10(D)) 
One Bay / Two Story Frame 
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Fig. 73 - 2x5 Grillage (Problem 7) 




VOLUME (in 3 ) 



Fig. 74 - Iteration History for Problem 7, Case A, Run 1 (Option 1(P)) 
2x5 Grillage 
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VOLUME (in 



Fig. 76 - Iteration History for Problem 7, Case A, Run 3 (Option 10(P)) 
2x5 Grillage 


(c u !) awmoA 



Fig. 77 - Iteration History for Problem 7, Case A, Run 4 (Option 10(D)) 
2*5 Grillage 
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VOLUME (in d ) 
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15,000 \- 
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ANALYSIS NUMBER 



Fig. 78 - Iteration History for Problem 7, Case A, Run 5 (Option 1(PU)) 
2x5 Grillage 


VOLUME (in 



Fig. 79 - Iteration History for Problem 7, Case B, Run 1 (Option 3(P)) 
2x5 Grillage 


( P u n awmoA 



Fig. 80 - Iteration History for Problem 7, Case B, Run 2 (Option 6 (P) ) 
2x5 Grillage 
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(p u !) 3WmOA 



Fig. 81 - Iteration History for Problem 7, Case B, Run 3 (Option 12(P)) 
2x5 Grillage 
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20,000 
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O CONSTRAINT VIOLATION < 1% 
□ CONSTRAINT VIOLATION < 10% 
• CONSTRAINT VIOLATION > 10% 


ANALYSIS NUMBER 


Fig. 82 - Iteration History for Problem 7, Case B, Run 4 (Option 3(PU)) 
2x5 Grillage 


r 



Fig. 83 - Two Bay / Six Story Frame (Problem 8) 
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WEIGHT (lb) 



ANALYSIS NUMBER 

Fig. 84 - Iteration History for Problem 8, Case A, Run 1 (Option 1(P)) 
Two Bay / Six Story Frame 
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Fig. 85 - Iteration History for Problem 0, Case A, Run 2 (Option 3(P)) 
Two Bay / Six Story Frame 
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Fig, 86 - Iteration History for Problem 8, Case A, Run 3 (Option 6(P)) 
Two Bay / Six Story Frame 
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Fig. 87 - Iteration History for Problem 8, Case A, Run 4 (Option 12 (P)) 
Two Bay / Six Story Frame 
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Fig. 88 - Iteration History for Problem S, Case A, Run 5 (Option 3(PU)) 
Two Bay J Six Story Frame 
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Fig. 90 - Iteration History for Problem 8, Case B, Run 2 (Option 3(P)) 
Two Bay /Six Story Frame 
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Fig. 91 - Iteration History for Problem 8, Case B, Run 3 (Option 4(P)) 

Two Bay / Six Story Frame 
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Fig. 92 - Iteration History for Problem 8, Case B, Run 4 (Option 10(P)) 
Two Bay / Six Story Frame 
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Fig. 93 - Iteration History for Problem 8, Case B, Run 5 (Option 1(PU)) 
Two Bay / Six Story Frame 
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Fig. 94 - Helicopter Tail Boom (Problem 9) 
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Fig. 97 - Iteration History for Problem 9, Run 3 (Option 10(P)) 
Helicopter Tail Boom 
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1 



Fig. 98 - Iteration History for Problem 9, Run 4 (Option 10(D)) 
Helicopter Tail Boom 



WEIGHT (lb) 


O CONSTRAINT VIOLATION < 1% 
□ CONSTRAINT VIOLATION < 10% 
• CONSTRAINT VIOLATION > 10% 



ANALYSIS NUMBER 

Iteration for Problem 9, Run 5 (Option 1(PU)) 
Helicopter Tail Boom 
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Fig. C4 - Buckling Stresses on a Long Simply Supported Plate 



3 
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Stress and Local Buckling Constraint 
Evaluation Points 





i y 



Cross Sectional Dimensions 


y 



Stress and Buckling Constraint Evaluation Points 


Fig. C7 - Cross Section for Element Type 14 
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Fig. E3 - Angle Representation of Beam Element Orientation 
(Special Case) 
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Table 1* Descriptions of Design Problem Solution Options 


Option* 

Design Space 

Approximat ions 

Element Force Variation 



(objective function/behavior 
constraints/side constraints) 


1 

RSP 

None/Lin ear /Linear 

Invariant 

2 

RSP 

None/Linear /Lin ear 

Local 

3 

RSP 

None/ Lin ear /Linear 

Global 

4 

RSP 

None/Hybrid/Linear 

Invariant 

5 

RSP 

None/Hybrid/Linear 

Local 

6 

RSP 

None/Hybrid/Linear 

Global 

7 

CSD 

Lin ear /Lin ear /None 

Invariant 

8 

CSD 

Linear /Linear /None 

Local 

9 

CSD 

Linear/ Linear /None 

Global 

10 

CSD 

Hybrid/Hybrid /None 

Invariant 

11 

CSD 

Hybrid/Hybrid /None 

Local 

12 

CSD 

Hybrid/Hybrid/None 

Global 


* The letter designation following the option number indicates the type of optimization 
method used to solve the approximate problem ((P)=Primal, (D)=Dual) • Also, the desig- 
nation (U) indicates that the approximate problem update feature was employed. 


Table 2. Definition of Problem 1 
Tied Cantilevered Beam 


Material Properties 
Young's Modulus 
Shear Modulus 
Poisson's Ratio 
Weight Density 

Allowable Normal Stress (Truss) 
Allowable Normal Stress (Beam) 
Allowable Shear Stress (Beam) 


E = 30.0 x 10 6 PSI 

G = 11.5 x 10 6 PSI 

V = .3 

p = .284 lb/in 3 

0 = 120,000 PSI 

a 

a = 20,000 PSI 

a 

T = 10,000 PSI 

a 


Nodal Loading 


Load 

Case 

Node 

No. 

Loading Components (lb, in-lb) 

F 

X 

F y 

F 

z 

M x 

M 

y 

M 

z 

1 

2 

0 . 

-10,000 

0 . 

0 . 

0 . 

0 . 


Uniform Loading 


Load 

Case 

Member 

No, 

Loading Components (lb/in, in-lb/in) 

P x 

p y __ 

P z 

m 

X 

n 

m 

z 

2 

i 

0 . 

-83.33 

0 . 

0 . 

0 . 

0 . 


Initial Design and Side Constraints 


Member 

No. 

Sizing 

Variable 

Initial 
Value2 
(in, in ) 

Lower 
Bound2 
(in, in ) 

Upper 
Bound2 
(in, in ) 

i 

B 

5.00 

1.00 

10.00 

2 

A 

o 

CM 

.01 

1.00 
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Table 3. Iteration History Data for Problem 1 
Tied Cantilevered Beam 



Weight (lb) [Maximum Constraint Violation (%)] 

Analysis 

Run 1 

Run 2 

Run 3 

Run 4 

No 

• 

Option* 1(P) 

Option* 2(P) 

Option* 3(P) 

Option* 6(P) 

0 


861.64 [0] 

861.64 [0] 

861.64 [0] 

861.64 [0] 

1 


697.00 [0] 

697.00 [0] 

697.00 [0] 

697.00 [0] 

2 


563.97 [1.2] 

563.92 [2.2] 

563.98 [1.0] 

563.98 [1.0] 

3 


531.34 [1.4] 

519.21 [1.8] 

476.97 [3.1] 

476.97 [3.1] 

4 


523.20 [.5] 

520.14 [0] 

475.18 [0] 

475.18 [0] 

5 


520.14 [.2] 

520.14 [0] 

473.00 [0] 

473.00 [0] 

6 


520.15 [0] 

520.14 [0] 

473.04 [0] 

473.04 [0] 

7 


520.15 [0] 


473.04 [0] 

473.04 [0] 

CPU 

Tot. 

.800 

.797 

1.035 

1.033 

Time 

Anal . 

.131 

.196 

.249 

.249 

(sec) 

o 

rt 

• 

.166 

.142 

.192 

.192 


*See Table 1 



Table 4. Final Designs for Problem 1 
Tied Cantilevered Beam 



Final Design (in, in^) 

Linking 

Member 

Sizing 

Ref. 51 

Run 1 

Ref. 51 

Run 2 

Group 

Nos . 

Variable 

Method I 

Option* 1 (P) 

Method II-B 

Option* 2(P) 

1 

i 

B 

3.8850 

3.8874 

3.8819 

3.8874 

2 

2 

A 

.1061 

.1066 

.1062 

.1066 

Weight (lb) 


519.40 

520.15 

516.80 

520.14 

Number of Analyses 

— 

8 

— 

7 


*See Table 1 
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Table 4. Final Designs for Problem 1 

Tied Cantilevered Beam (cont.) 



Final 

Design (in, in^) 

Linking 

Member 

Sizing 

Ref. 51 

Run 3 

Run 4 

Group 

Nos. 

Variables 

Method IV -B 

Option* 3(P) 

Option* 6(P) 

i 

i 

B 

3.6394 

3.6421 

3.6421 

2 

2 

A 

.4571 

.4354 

.4354 

Weight (lb) 


473.40 

473.04 

473.04 

Number of Analyses 

— 

8 

8 


*See Table 1 
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Table 6. Definition of Problem 2 
Two Member Frame 


Material Properties 
Young's Modulus : 
Shear Modulus : 
Poisson's Ratio : 
Mass Density : 
Allowable Stress : 


E = 20.74 x 10 6 N/cm 2 
G = 7.97 x 10^ N/cm 2 

V = .3 

P = 2.77 x 10" 2 kg/cm 3 

4 2 

0 * 2.76 x 10 N/cm 

a 


Nodal Loading 


Load 

Case 

Node 

Nos. 

Loading Components (N, N-cm) 

F 

X 

F 

y 

F 

z 

M 

X 

mm 

M 

z 

1 

2 

0 . 

-44480. 

0 . 

0 . 

0 . 

0 . 


Initial Design and Side Constraints 


Member 

Sizing 

Initial 

Lower 

Upper 

No. 

Variable 

Value (cm) 

Bound (cm) 

Bound (cm) 


B 

15.20 

6.350 

25.40 


c b 

2.03 

.254 

2.54 

i 

H 

10.20 

6.350 

25.40 


C h 

2.29 

.229 

2.54 


B 

22.90 

6.350 

25.40 


C b 

2.03 

.254 

2.54 

2 

H 

20.30 

6.350 

25.40 


C h 

2.29 

.229 

2.54 
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Table 7. Iteration History Data for Problem 2, Case A 
Two Member Frame 



Mass (kg) [Maximum Constraint Violation (%) ] 

Analysis 

Run 1 

.Run 2 

Run 3 

Run 4 

No 

• 

Option* 1 (P) 

Option* 2(P) 

Option* 3 (P) 

Option* 4(P) 

0 


1220.78 [0] 

1220.78 [0] 

1220.78 [0] 

1220.78 [0] 

1 


531.31 [0] 

531.31 [0] 

531.31 [0] 

295.74 [0] 

2 


313.48 [0] 

313.48 [0] 

313.48 [0] 

188.09 [0] 

3 


204.35 [0] 

204.35 [0] 

204.35 [0] 

163.20 [3.5] 

4 


165.10 [24.6] 

165.23 [23.7] 

165.24 [23.6] 

151.29 [3.5] 

5 


161.08 [8.0] 

158.88 [6.3] 

159.99 [2.5] 

141.61 [.6] 

6 


162.68 [1.0] 

150.79 [2.8] 

152.69 [0] 

133.80 [.3] 

7 


151.29 [6.8] 

142.97 [2.2] 

143.70 [0] 

133.88 [0] 

8 


144.37 [2.4] 

136.46 [1.4] 

136.94 [0] 

130.20 [.3] 

9 


137.49 [1.3] 

132.14 [.8] 

136.88 [0] 

130.32 [0] 



132.09 [.3] 

130.20 [.3] 

132.35 [0] 

130.32 [0] 

11 


132.04 [0] 

130.20 [.3] 

130.47 [0] 


12 


130.25 [.1] 

130.20 [.3] 

130.46 [0] 


13 


130.25 [.1] 


130.46 [0] 


14 


130.25 [.1] 




msm 

Tot. 

2.186 

2.031 

2.192 

1.846 


Anal . 

.293 

.375 

.467 

.217 

Mill 

Opt . 

.999 

.863 

.875 

.894 


*See Table 1 






Table 7. 




1 Analysis 

Run 5 

No 

• 

Option* 7 (P) 

0 


1220.78 [0] 

1 


686.69 [0] 

2 


397.36 [0] 

3 


285.02 [33.5] 

4 


234.82 [25.0] 

5 


193.72 [20.7] 

6 


166.19 [0] 

7 


139.00 [5.9] 

8 


128.55 [8.2] 

9 


128.13 [8.4] 

10 


127.95 [8.2] 

11 


129-54 [2.3] 

12 


130.28 [.1] 

13 


130.26 [0] 

14 


130.26 [0] 

CPU 

mm 

1 .915 

Time 


.236 

(sec) 

1 

.882 


* See Table 1 

**Constraint was not Retained 


Iteration History Data for Problem 2, Case A 
Two Member Frame (cont.) 




Mass (kg) [Maximum Constraint Violation (%) ] 


Run 6 

Run 7 

Run 8 

Option* 10(P) 

Option* 10(D) 

Option* 1(PU) 

1220.78 [0] 

1220.78 [0] 

1220.78 [0] 

686.69 [0] 

686.69 [0] 

313.48 [0] ** 

397.36 [0] 

397.36 [0] 

165.10 [24.6] 

298.48 [3.6] 

298.52 [3.5] 

149.25 [8.1] 

247.78 [0] 

247.69 [0] 

137.26 [3.3] 

204.49 [0] 

204.45 [0] 

129.92 [1.2] 

168.88 [0] 

168.95 [0] 

130.40 [0] 

142.71 [0] 

142.69 [0] 


132.08 [.1] 

132.14 [0] 


132.06 [.1] 

132.09 [0] 


132.06 [.1] 

132.05 [0] 


1.451 

1.146 

2.078 

.177 

.173 

.449 

.596 

.329 

.884 
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Table 8. Final Designs for Problem 2, Case A 
Two Member Frame 


Group I Nos. I Variables 


Mass (kg) 

Number of Analyses 


*See Table 1 

■PSizing Variable at Upper Bound 
-Sizing Variable at Lower Bound 




Final Design ( 

cm) 


Ref. 29 

Run 1 

Option* 1(P) 

Run 2 

Option* 2(P) 

Run 3 

Option* 3(P) 

Run 4 

Option* 4(P) 






6.35" 

6.35“ 

6.35 

6.35 

6.35“ 






.229~ 

.229" 

.229“ 

.229“ 

.229" 

6.35- 

6.35" 

6.35" 

6.35" 

6.35" 

.254" 

.254“ 

.254“ 

.254" 

.254" 

6.35" 

25.36 

25.30 

25.34 

25.36 

1.14 

.248 

.248 

.249 

.249 

25.40 + 

25.38 

25.40 + 

25.40 + 

25.29 






.254" 

.254" 

.254“ 

.254“ 

.254“ 

133.70 

130.25 

130.20 

130.46 

130.32 

19 

15 

13 

14 

11 
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Table 8. Final Designs for Problem 2 , Case A 
Two Member Frame (cont.) 



Final Design (cm) 

- - - - V 

Linking 

Member 


Sizing 

Run 5 

Run 6 

Run 7 

Run 8 

Group 

Nos. 

Variables 

Option* 7 (P) 

Option* 10(P) 

Option* 10(D) 

Option* 1(PU) 




B 

6.35" 

6.35" 

6.35" 

6.35" 

i 

i 


c b 

.229" 

.229" 

.229" 

.229" 




H 

6.35" 

6.35" 

6.35~ 

6.35" 




c b 

.254" 

.254" 

.254" 

.254" 




B 

25.40 + 

9.48 

9.63 

25.31 




fc b 

.248 

.715 

.703 

.249 

2 

2 


H 

25.40 + 

25.40 + 

25.40 + 

25.40 + 




C h 

.254" 

.254“ 

.254" 

.254" 

Mass (kg) 



130.26 

132.06 

132.05 

130.42 

Number of Analyses 


15 

11 

11 

7 


*See Table 1 


+Sizing Variable at Upper Bound 
-Sizing Variable at Lower Bound 
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Table 9, Critical Constraints for Problem 2, Case A 
Two Member Frame 


Run 

No. 

Option 

M * 

No . 

Stress Constrained 
Members 

Ref. 29 

i 

i 

2 

1 

up) 

2 

2 

2(P) 

2 

3 

3(P) 

2 

4 

4(P) 

2 

5 

7(P) 

2 

6 

10(P) 

2 

7 

10(D) 

2 

8 

1(PU) 

2 


*See Table 1 
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Table 10. 


Analysis 

No. 

Run ] 
Option* 

1 (P) 

0 


1220.78 

[0] 

1 


531.31 

[0] 

2 


313.48 

[0] 

3 


204.35 

[42.2]** 

4 


182.90 

[0] 

5 


159.41 

[6.7] 

6 


153.69 

[2.6] 

7 


145.53 

[1.7] 

8 


138.31 

[1.4] 

9 


132.27 

[1.4] 

10 


130.22 

[3.6] 

11 


130.79 

[.1] 

12 


130.82 

[0] 

13 




14 




15 




CPU 

Time 

(sec) 

1 

2.051 

.276 

.940 


*See Table 1 


**Constra int 
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Table 10. Iteration History Data for Problem 2, Case B 
Two Member Frame (cont.) 



Mass (kg) [Maximum Constraint Violation (%)] 

Analysis 

Run 5 

Run 6 

Run 7 

No 


Option* 10(P) 

Option* 10(D) 

Option* 1(PU) 

0 


1220.78 [0] 

1220.78 [0] 

1220.78 [0] 

1 


686.69 [0) 

686.69 [0] 

313.48 [0] 

2 


397.37 [0] 

397.36 [0] 

165.10 [>100]** 

3 


298.48 [3.6] 

299.21 [2.0] 

149.55 [77.6] 

4 


247.78 [0] 

247.73 [0] 

135.39 [19.6] 

5 


204.49 [0] 

204.42 [0] 

130.35 [2.5] i 

6 


168.88 ro] 

168.95 [0] 

130.75 [0] 

7 


142.71 [0] 

142.73 [0] 


8 


132.08 [.1] 

132.15 [0] 


9 


132.06 [.1] 

132.10 [0] 


10 


132.06 £ .1] 

132.15 [0] 


11 





12 





13 





14 





15 





CPU 

Tot. 

1.526 

1.233 

2.274 

Time 

Anal . 

.182 

.177 

.465 

(sec) 

Opt . 

.602 

.350 

.977 


*See Table 1 **Constraint was not retained 








ns for Problem 2, Case B 
Frame 


Final Design (cm) 


Run 1 

Run 2 

Run 3 

Option* 1(P) 

Option* 2(P) 

Option* 3(P) 

6.35" 

6.35" 

6.35" 

.229" 

.229" 

.229" 

6.35" 

6.35“ 

6.35" 

.254" 

.254" 

.254" 

18.41 

18.36 

18.41 

.348 

.348 

.348 

25.40 + 

25.40 + 

25.40 + 

.254" 

.254“ 

.254" 


130.82 130.60 130.75 


13 


13 


13 
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Table 11. Final Designs for Problem 2, Case B 
Two Member Frame (cont.) 


Linking 

Group 

Member 

Nos. 

i 

i 

2 

2 


Mass (kg) 

Number of Analyses 


*See Table 1 

+Sizing Variable at Upper Bound 
-Sizing Variable at Lower Bound 


Final Design (cm) 

Run 4 

Option* 4(P) 

Run 5 

Option* 10(P) 

Run 6 

Option* 10(D) 

Run 7 

Option* 1(PU) 

6.35" 

6.35" 

6.35" 

6.35" 

.229" 

.229" 

.229" 

.229" 

6.35~ 

6.35" 

6.35" 

6.35" 

.254“ 

.254" 

.254" 

.254" 

16.74 

9.48 

9.64 

18.35 

.385 

.715 

.704 

.349 

25.38 

25.40 + 

25.40 + 

25.40 + 

.254" 

.254" 

.254" 

.254" 

130.81 

132.06 

132.15 

130.74 

10 

11 

i. - 

11 

7 





























Table 12, Critical Constraints for Problem 2, Case B 
Two Member Frame 


Run 

No. 

Option 

No.* 

Stress Constrained 
Members 

Local Buckling 
Constrained Members 

Ref. 29 

— 

2 

— 

1 

1 (P) 

2 

2 

2 

2(P) 

2 

2 

3 

3(P) 

2 

2 

4 

4(P) 

2 

— 

5 

10(P) 

2 

— 

6 

10(D) 

2 

— 

7 

1(PU) 

2 

2 


*See Table 1 





Table 13. Defintion of Problem 3 
Three Member Frame 


Material Properties 
Young’s Modulus : 
Shear Modulus : 
Poisson’s Ratio : 
Allowable Stress : 


E = 30.0 x 10 6 PSI 
G = 11.5 x 10 6 PSI 
V - .3 

a = 40,000 psi 

a 


Nodal Loading 


Load 

Case 

Node 

No. 


Loading Components (lb, 

in-lb) 


F 

X 

urn 

F 

z 

M 

X 

mm 

M 

z 


2 

0 . 

0 . 

-10000 

0 . 

0 . 

0 . 

i 









3 

0 . 

0 . 

-10000 

0 . 

0 . 

0 . 


Initial Design and Side Constraints 


Member 

Nos. 

Sizing 

Variable 

Initial 
Value (in) 

Lower 
Bound (in) 

Upper 
Bound (in) 


B 

9.0 

2.5 

10.0 

1-3 

H 

9.0 

2.5 

10.0 


t 

.9 

. 1 

1.0 
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Table 14. Iteration History Data for Problem 3 
Three Member Frame 



3 

Volume (in ) [Maximum Constraint Violation (%)] 

Analysis 

Run 1 

Run 2 

Run 3 

Run 4 

No. 


Option* 1(P) 

Option* 2 (P) 

Option* 3(P) 

Option* 4 (P) 

0 


8748.00 [0] 

8748.00 [0] 

8748,00 [0] 

8748,00 [0] 

1 


3138.34 [ 0 ] 

3138.34 [ 0 ] 

3138.34 [0] 

1976.35 [40.3]** 

2 


1808.55 [48.7] 

1600.01 [ 100 ] 

1808,51 [48.7] 

1895.92 [8.0] 

3 


1857.16 [10.1] 

1797.94 [19.1] 

1857.60 [10.1] 

1869.33 [0] 

4 


1867.32 [.7] 

1880.41 [ 0 ] 

1867.53 [.7] 

1806.43 [ 0 ] 

5 


1850.23 [0] 

1853.33 [ 0 ] 

1824.50 [0] 

1758.83 [ 0 ] 

6 


1779.38 [0] 

1774.28 [.7] 

1780.04 [ 0 ] 

1720.58 [ 0 ] 

7 


1744.35 [0] 

1742.52 [. 2 ] 

1746.63 [ 0 ] 

1690.07 [ 0 ] 

8 


1714.93 [ 0 ] 

1715.20 [0] 

1716.11 [0] 

1669.84 [ 0 ] 

9 


1690.33 [0] 

1690.44 [0] 

1691 39 [ 0 ] 

1668.23 [.1] 

10 


1673.16 [0] 

1682.79 [0] 

1682.85 [0] 

1667.82 [.1] 

11 


1670.60 [0] 

1666.43 [0] 

1667.96 [0] 


12 


1670.18 [0] 

1666.43 [ 0 ] 

1667.96 [0] 


13 


1670.18 [0] 

1666.43 [ 0 ] 

1667.96 [ 0 ] 


CPU 

Tot. 

3.189 

3,423 

3.884 

2.591 

Time 

Anal. 

.335 

.675 

.866 

.262 

(sec) 

Opt. 

1.935 

1 .820 

2.114 

1.551 


*See Table 1 **Constraint was not Retained 
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1 


Table 14, Iteration History Data for Problem 3 
Three Member Frame (cont.) 



3 

Volume (in ) [Maximum Constraint Violation (%) ] 

Analysis 

Run 5 

Run 6 

Run 7 

No 


Option* 10(P) 

Option* 10(D) 

Option* 1(PU) 

0 


8748.00 [0] 

8748.00 [0] 

8748.00 [0] 

1 


4286.52 [3.5]** 

4286.52 [3.5]** 

1590.87 |>100]** 

2 


3066.00 [0] 

3084.45 [0] 

1751.41 [16.7] 

3 


2312.58 [0] 

2311.07 [0] 

1773.40 [1.2] 

4 


1982.69 [0] 

1981.12 [0] 

1715.05 [0] 

5 


1846.23 [0] 

1841.54 [0] 

1674.29 [0] 

6 


1739.91 [0] 

1727.85 [0] 

1671.03 [0] 

7 


1680.03 [0] 

1673.63 [0] 


8 


1668.90 [0] 

1666.48 [0] 


9 


1666.65 [0] 

1666.66 [0] 


10 


1666.65 [0] 

1666.71 [0] 


11 


1666.65 [0] 



12 





13 





cru 

Tot. 

2.437 

1 .449 

3.321 

Time 

Anal . 

.243 

.221 

.507 

(sec) 

Opt . 

1.421 

.516 

1 .982 


*See Table 1 **Constraint was not Retained 
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Table 15. Final Designs for Problem 3 
Three Member Frame 



Final Design (in) 

Linking 

Group 

Member 

Nos. 

S iz ing 
Variables 

Ref. 52 

Run 1 

Option* 1(P) 

Run 2 

Option* 2(P) 

Run 3 

Option* 3(P) 



B 

10.00 + 

10.00 + 

9.99 

9.99 

1 

i 

H 

+ o 

o 

o 

r-H 

10.00 + 

10.00 + 

10.00 + 



t 

.199 

.201 

.200 

.201 



B 

2.50" 

2 .50“ 

2.50“ 

2.50“ 

2 

2 

H 

2.50" 

2.50" 

2.50" 

2.50“ 



t 

.100" 

.100" 

.100" 

.100“ 



B 

10.00 + 

10.00 + 

10.00 + 

10.00 + 

3 

3 

H 

10.00 + 

9.99 

10.00 + 

9.99 



t 

.199 

.201 

.200 

.200 

Volume (in^) 


1656.96 

1670.18 

1666.43 

1667.96 

Number of Analyses 

14 

14 

14 

14 


*See Table 1 

+Sizing Variable at Upper Bound 
-Sizing Variable at Lower Bound 
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Table 15. Final Designs for Problem 3 
Three Member Frame (cont.) 



Final Design (in) 

Linking 

Group 

Member 
Nos . 

Sizing 

Variables 

Run 4 

Option* 4 (P) 

Run 5 

Option* 10(P) 

Run 6 

Option* 10(D) 

Run 7 

Option* 1(PU) 




B 

10.00 + 

10.00 + 

10.00 + 

10.00 + 

i 

1 


H 

9.98 

10.00 + 

10.00 + 

9.98 




t 

.201 

.200 

.200 

.201 




B 

2.50" 

2.50" 

2.50" 

2.50" 

2 

2 


H 

2.50" 

2.50“ 

2.50" 

2.50" 




t 

.100" 

.100" 

.100" 

.100" 




B 

10.00 + 

10.00 + 

10.00 + 

10.00 + 

3 

3 


H 

9.98 

10.00 + 

10.00 + 

9.99 




t 

.201 

.200 

.200 

.201 

Volume (in^) 



1667.82 

1666.65 

1666.72 

1670.77 

Number of Analyses 


11 

12 

11 

7 


*See Table 1 

+Sizing Variable at Upper Bound 
-Sizing Variable at Lower Bound 
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Table 16. Critical Constraints for Problem 3 
Three Member Frame 


Run 

No. 

Option 

No. 

Stress Constrained 
Members 

Ref. 52 

— 

1,3 

1 

1(P) 

1,3 

2 

2(P) 

1,3 

3 

3(P) 

1,3 

4 

4 (P) 

1,3 

5 

10(P) 

1,3 

6 

10(D) 

1,3 

7 

1(PU) 

1,3 


*See Table 1 





Table 17. Definition of Problem 4 
Seven Member Frame 


Material Properties 
Young f s Modulus : 
Shear Modulus : 
Poisson’s Ratio : 
Mass Density : 
Allowable Stress : 


E = 20.74 x 10^ N/cm^ 

G = 7.85 x 10 6 N/cra 2 

V = .32 

p = 7.81 x 10 2 kg/cm 2 
2 

a = 20,000 N/cm 
a 


Nodal Loading 


Load 

Node 


Loading Components (N, 

N-cm) 


Case 

No. 

F 

F 

F 

M 

M 

M 



X 

y 

z 

X 

y 

z 

1 

3 

-40000 

-40000 

0 . 

0 . 

0 . 

0 . 


3 

-50000 

0 . 

0 . 

0 . 

0 . 

0 . 

2 

4 

0 . 

-50000 

0 . 

0 . 

0 . 

0 . 


Displacement Constraints 


Load 

Case 

Node 

No. 

Direction 

Lower Bound 

Upper Bound 

i 

3 

y 

-.2 cm 

. 04 cm 

2 

3 

y 

-.2 cm 

.04 cm 


Initial Design and Side Constraints 


Member 
Nos . 

Sizing 

Variable 

Initial 
Value (cm) 

Lower 
Bound (cm) 

Upper 
Bound (cm) 


B 

7.62 

1.00 

10.00 

1-7 

H 

7.62 

1.00 

10.00 


t 

.18 

.076 

.30 
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Table 18. Iteration History Data for Problem 4 
Seven Member Frame 



Mass (kg) [Maximum Constraint Violation (%) ] 

Analysis 

No. 

Run 1 

Option* 1 ( P) 

Run 2 

Option* 2(P) 

Run 3 

Option* 3(P) 

Run 4 

Option* 6(P) 

0 


16.26 

[11.9] 

16.26 

[11.9] 

16.26 

[11.9] 

16.26 

[11.9] 

1 


9.52 

[63.8]** 

9.75 

[76.7]** 

9.93 

[74.6]** 

10.07 

[66.3]** 

2 


9.25 

[20.7] 

8.59 

[41.9] 

8.40 

[20.7] 

8.59 

116.6] 

3 


8.75 

[62.5] 

8.37 

[21.2] 

8.77 

[1.5] 

8.85 

[-2] 

4 


8.49 

[100] 

8.49 

[5.1] 

8.56 

[.8] 

8.77 

1.3] 

5 


8.87 

[94.4] 

8.53 

[0] 

8.57 

[0] 

8.65 

[.1] 

6 


9.26 

[62.3] 

8.38 

[2.9] 

8.43 

[.1] 

8.55 

[0] 

7 


9.78 

[31.3] 

8.34 

[1.5] 

8.31 

[.2] 

8.52 

[0] 

8 


10.03. 

[15.1] 

8.36 

[0] 

8.27 

[0] 

8.40 

[0] 

9 


10.02 

[13.5] 

8.29 

[.5] 

8.23 

[0] 

8.31 

f.l] 

10 


1C. 22 

[7.8] 

8.23 

[3.7] 

8.12 

[.9] 

8.28 

[.1] 

11 


10.20 

[6.9] 

8.29 

[0] 

8.14 

[-2] 

8.16 

[.41 

12 


10.14 

[6.0] 

8.25 

[5.9] 

8.11 

[0] 

8.18 

[0] 

13 


10.06 

[5.7] 

8.12 

[3.4] 

8.10 

[0] 

8.12 

[-6] 

14 


10.00 

[4.9] 

8.19 

[0] 

8.08 

[0] 

8.12 

[0] 

15 


9.94 

[4.5] 

8.17 

[-7] 

8.08 

[0] 

8.10 

[.2] 

16 


9.86 

[4.7] 






17 


9.82 

[3.8] 







18 


9.78 

[2.9] 







19 


9.75 

[2.1] 







20 


9.72 

[1.3] 







CPU 

Time 

(sec) 

Tot. 
Anal . 
Opt . 

No Convergence 

10.627 

1.963 

6.369 

10.615 

3.388 

4.921 

10.929 

3.291 

5.322 


*See Table 1 


**Constraint was not Retained 


Table 13. Iteration History Data for Problem 4 
Seven Member Frame (cont.) 



Mass (kg) 

[Maximum Constraint Violation (%)] 

Analysis 

Run 5 

Run 6 

Run 7 

No. 

Option* 12(P) 

Option* 12(D) 

Option* 3 (PU) 

0 

16.26 [11.9] 

16.26 [11.9] 

16.26 [11.9] 

1 

9.70 [100]** 

10.60 [81.8]** 

8.62 [95.5]** 

2 

8.70 [39.6]** 

9.00 [10.8] 

8.47 [5.6] 

3 

8.97 [.1] 

8.65 [0] 

8.37 [0] 

4 

8.50 [.2] 

8.24 [1.3] 

8.30 [0] 

5 

8.30 [0] 

8.13 [1.4] 

8.23 [.1] 

6 

8.14 [0] 

8.09 [ .7] 

8.21 [0] 

7 

8.13 [0] 

8.07 [ .8] 

8.10 [.3] 

8 

8.12 [0] 

8.08 [.3] 

8.09 [0] 

9 

8.10 [0] 

8.09 [.1] 


10 

8.10 [0] 



11 




12 




13 




14 




15 




16 




17 




18 




19 




20 




CPU 

Tot. 

8.215 

17.090 

10.333 

Time 

Anal . 

2.175 

1.945 

2.589 

(sec) 

Opt . 

4.459 

13.730 

5.307 


*See Table 1 


**Constraint was not Retained 


Table 19. Final Designs for Problem 4 
Seven Member Frame 



Final Design (cm) 

Linking 

Group 

Member 

Nos. 

Size 
Var . 

Ref. 28 

Run 2 

Option* 2(P) 

Run 3 

Option* 3 (P) 

Run 4 

Option* 6(P) 



B 

10.00 + 

9.96 

8.39 

7.22 

1 

i 

H 

3.73 

2.73 

5.18 

5.09 



t 

.165 

.162 

.173 

.192 



B 

10.00 + 

10.00* 

7.91 

8.06 

2 

2 

H 

10.00* 

10.00 + 

2.10 

2.97 



t 

.096 

.099 

.199 

.180 



B 

7.81 

9.97 

9.99 

9.97 

3 

3 

H 

10.00 + 

10.00* 

1.73 

2.49 



t 

.110 

.096 

.158 

.149 



B 


10.00 + 

10 00* 

9.99 

4 

4 

H 

mm 

8.06 

3.08 

3.56 



t 

1 

.173 

.201 

.198 



B 

1.50 

1 .07 

1.76 

2.13 

5 

5 

H 

1.32 

1 .00' 

1.06 

1.29 



t 

.076' 

.076' 

.076" 

.076" 



B 

1.28 

1.00" 

1.00" 

1.29 

6 

6 

H 

1.28 

1 .00" 

1.00" 

1 .oo- 



t 

.076" 

.076" 

.076" 

.076" 



B 

5.31 

3.02 

2.91 

3.06 

7 

7 

H 

5.13 

7.56 

10.00* 

10.00* 



t 

• 076~ 

.087 

.081 

.079 

Mass (kg) 

Number of Analyses 

8.25 

8.17 

16 

i 

8.08 

16 

8.10 

16 


*See Table 1 +Sizing Variable at Upper Bound 

-Sizing Variable at Lower Bound 
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Table 19. Final Designs for Problem 4 
Seven Member Frame (cont.) 



Final Design (cm) 

Linking 

Group 

— 
Member 
Nos . 

Size 

Var. 

Run 5 

Opt ion* 12 (P) 

Run 6 

Option* 12(D) 

Run 7 

Option* 3(PU) 



B 

7.62 

8.13 

10.00 + 

1 

i 

H 

A. 64 

4.18 

5.52 



t 

.191 

.184 

.153 



B 

6.49 

6.88 

9.17 

2 

2 

H 

6.00 

5.56 

1.43 



t 

.155 

.158 

.188 



B 

7.72 

7.49 

9.36 

3 

3 

H 

3.98 

3.81 

1.28 



t 

.163 

.170 

.173 



B 

7.92 

8.88 

9.96 

4 

4 

H 

6.14 

5.40 

2.16 



t 

.202 

.199 

.214 



B 

1.18 

1.07 

1.80 

5 

5 

H 

1.17 

1.11 

1.12 



t 

.076" 

.076" 

.076" 



B 

1.00“ 

1 .00' 

1.15 

6 

6 

H 

1.00~ 

1 .00' 

1.00“ 



t 

.076" 

.076' 

.076“ 



B 

4.95 

4.52 

2.99 

7 

7 

H 

5.34 

4 .65 

10.00 + 



t 

.092 

.106 

.081 

Mass (kg) 

Number of Analyses 

8.10 

11 

8.09 

10 

8.09 

10 


*See Table 1 +Sizing Variable at Upper Bound 

-Sizing Variable at Lower Bound 
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Table 21 


Run 

Option 

No. 

No.* 



Ref. 28 

— 

2 

2(P) 

3 

3(P) 

4 

6(P) 

5 

12(P) 

6 

12(D) 

7 

3(PU) 


*See Table 1 



















Table 21 


Definition of Problem 5 
Portal Frame 


Material Properties 
Young’s Modulus 
Shear Modulus 
Poisson’s Ratio 
Allowable Normal Stress 
Allowable Shear Stress 


E = 7.0 x 10 6 N/cm 2 
G = 2.7 x 10 6 N/cm 2 
v = .3 

a = 2.0 x 10 4 N/cm 2 
a 

T = 1.16 x 10 4 N/cm 2 


Nodal Loading 


Load 

■ 

Loading Components (N, 

N-cm) 


Case 

1 

F 

mgm 


M 

■1 

M 



X 

mem 

■9 

X 

ma 

z 

i 

3 

4 

5 x 10 

0. 

0. 

0. 

0. 

0. 

2 

3 

0. 

0. 

0. 

0. 

0. 

1 

ro 

►— 

o 


Displacement Constraints 


Load 

Case 


Direction 

Lower Bound 

Upper Bound 

1 

3 

X 

-4.0 cm 

4.0 cm 

2 

3 

0 

Z 

-.015 rad 

.015 rad 
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Table 21. Definition of Problem 5 
Portal Frame (cont.) 


Initial Design and Side Constraints 



354 
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Table 22. Iteration History Data for Problem 5 
Portal Frame 



O 

Volume (cm ) [Maximum Constraint Violation (%) ] 

Analysis 

Run 1 

Run 2 

Run 3 

No 

• 

Option* 1 (P) 

Option* 2(P) 

Option* 3(P) 

0 


275,000 [0] 

275,000 [0] 

275,000 [0] 

1 


193,179 [0] 

193,179 [0] 

193,179 [0] 

2 


126,224 [84.0]** 

126,224 [84.0]** 

126,224 [84.0]** 

3 


119,299 [17.6] 

133,432 [10.9] 

121,745 [21.5] 

4 


102,977 [>100.0]** 

119,212 [12.5] 

102,563 [>100.0]** 

5 


101,275 [55.1] 

101,976 [>100]** 

101,018 [46.1] 

6 


104,134 [1.0] 

99,087 [70.3] 

97,564 [8.6] 

7 


99,506 [>100.0]** 

100,041 [12.9] 

97,444 [.3] 

8 


86,049 [>100.0]** 

100,527 [.8] 

97,459 [.2] 

9 


94,010 [>100.0]** 

100,547 [.1] 

97,460 [0] 

10 


91,575 [>100.0] 

99,083 [1.7] 


11 


93,696 [>100.0] 

99,110 [.1] 


12 


93,575 [37.2] 

97,880 [1.3] 


13 


93,955 [.2] 

97,934 [0] 


14 


94,091 [.1] 

97,560 [.6] 


15 


93,702 [.2] 

97,675 [0] 


16 


93,773 [.1] 

97,574 [0] 


17 


93,732 [.3] 

97,226 [.6] 


18 



97,365 [0] 


19 



97,239 [0] 


20 



97,240 [0] 


CPU 

Tot. 

5.339 

6.093 

3.759 

Time 

Anal . 

.956 

1.613 

.953 

(sec) 

Opt. 

2.825 

2.739 

1.804 


*See Table 1 **Constraint was not Retained 
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Table 22. Iteration History Data for Problem 5 
Portal Frame (cont.) 



3 

Volume (era ) [Maximum Constraint Violation (%)] 

Analysis 

Run A 

Run 5 

No. 


Option* 12(P) 

Option* 12(D) 

0 


275,000 [01 

275,000 [0] 

1 


136,028 [ 86 . 2 ]** 

108, A6A [>100]** 

2 


109,962 [>100]** 

95.0AA [31.1] 

3 


9A.857 [30.2] 

88,177 [ A 1 . 2 ]** 

A 


88,267 [A3. 8 ]** 

86,557 [0] 

5 


86,899 [0] 

8A.766 [0] 

6 


85.A69 [0] 

8A , 238 [. 8 ] 

7 


8A.82A [0] 

8A.109 [.7] 

8 


8A.A83 [91.0]** 

8A,022 [.5] 

9 


8A.265 [10.6] 

8A.056 [.7] 

10 


8A.215 [.1] 

8A,058 [.3] 

11 


8A ,272 [0] 


12 




CPU 

Tot. 

3.75A 

3.527 

Time 

Anal . 

1.0A2 

1.069 

(sec) 

Opt . 

1.672 

1.A36 


*See Table 1 


**Constraint was not Retained 
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Table 23. Final Designs for Problem 5 
Portal Frame 



Final Design (cm) 

Linking 

Member 

Size 

Ref. 53 

Run 1 

Run 2 

Run 3 

Run 4 

Run 5 

Group 

Nos. 

Var. 


Option*. 1 (P) 

Option* 2(P) 

Option* 3(P) 

Option* 1 2 (P) 

Option* 12(D) 



B 1 

13.00 

24.26 

22.90 

22.66 

11.55 

11.26 



C 1 

.450 

.624 

.567 

.550 

.415 

.410 



H 

74.90 

66.62 

66.83 

67.17 

77.86 

78.21 

1 

i 

t 3 

.497 

.475 

.461 

.460 

.523 

.523 



B 2 

12.10 

5.00" 

5 . 00~ 

5.00" 

10.40 

10.17 



C 2 

.487 

1.290 

1.496 

1 .610 

.463 

.456 



B 1 

11.40 

19.34 

18.26 

17.57 

11.63 

11.69 



h 

.404 

.463 

.463 

.455 

.410 

.417 



H 

89 .90 

87.81 

76.00 

72.92 

100.00 + 

99.47 

2 

2 

C 3 

.397 

.401 

.354 

.341 

.436 

.435 



B 2 

10.70 

5.00" 

5.00" 

5.00“ 

10.71 

10.94 



C 2 

.435 

1.236 

1.356 

1.430 

.446 

.447 
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Table 23. Final Designs for Problem 5 
Portal Frame (cont.) 


Group I Nos. I Var. 


Volume (cm ) 

Number of Analyses 


*See Table 1 

+Sizing Variable at Upper Bound 
-Sizing Variable at Lower Bound 




Final 

Design (cm) 



Ref. 53 

Run 1 

Option* 1 (P) 

Run 2 

Option* 2(P) 

Run 3 

Option* 3(P) 

Run 4 

Option* 12(P) 

Run 5 

Option* 12(D) 

7.50 

14.27 

5.00~ 

5.00" 

5 ,00~ 

5.00~ 

.268 

.356 

1.356 

1.354 

.142 

.143 

61.90 

36.35 

55.31 

59.21 

25.00" 

25.00“ 

.250 

.152 

.262 

.282 

.100" 

.100“ 

10.00" 

10.00" 

15.37 

15.24 

10.00“ 

10.00“ 

.369 

.619 

.535 

.541 

.276 

.276 

90,592 

93,732 

97,240 

97,460 

84,272 

84,058 

— 

18 

21 

10 

12 

11 



































Table 25. Definition of Problem 6 

One Bay / Two Story Frame 


Material Properties 


Young’s Modulus 


E = 

30.0 x 10 6 PSI 

Shear Modulus 


G = 

11.5 x 10 6 PSI 

Poisson’s Ratio 


V = 

.3 

Weight Density 


p = 

.2836 lb/in 3 

Yield Stress 


a = 
a 

; 36,000 PSI 

Factor of Safety 


FS = 

■ 1.51 


Nodal Loading 


Load 

Case 

Node 

No. 

Loading 

Components (1 

b, in-lb) 

F 

X 

mm 

F 

z 

M 

X 

mm 

M 

z 


2 

45000 

0 . 

0 . 

0 . 

0 . 

0 . 

2 









3 

45000 

0 . 

0 . 

0 . 

0 . 

0 . 


Uniform Loading 


Load 

Case 

Member 

No. 

Loading Components (‘ 

.b/in, in-lb/in) 

Px 

p y 


m 

X 

mm 

m 

z 


3 

0 . 

-500. 

0 . 

0 . 

0 . 

0 . 


4 

0 . 

-500. 

0 . 

0 . 

0 . 

0 . 

i 









5 

0 . 

-500. 

0 . 

0 . 

0 . 

0 . 


6 

0 . 

-500. 

0 . 

0 . 

0 . 

0 . 
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Table 25. Definition of Problem 6 

One Bay / Two Story Frame (cont.) 


Displacement Constraints 


Load 

Case 

Node 

No. 

Direction 

Lower Bound 

Upper Bound 


2 

X 

-.36 in 

.36 in 


3 

X 

-.72 in 

.72 in 


4 

X 

-.36 in 

.36 in 

2 

5 

X 

-.72 in 

.72 in 


6 

X 

-.36 in 

.36 in 


7 

X 

-.72 in 

.72 in 


Initial Design and Side Constraints 


Member 
Nos . 

Sizing 

Variable 

Initial 
Value (in) 

Lower 
Bound (in) 

Upper 
Bound (in) 

1-8 

R 

16.61 

1.00 

25.0 


t 

.45 

.01 

5.0 
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Table 26. Iteration History Data for Problem 6, Case A 
One Bay / Two Story Frame 




Weight (lb) [Maximum Constraint Violation (%)] 

Analysis 

Run 1 

Run 2 

Run 3 

No 


Option* 1(P) 

Option* 3(P) 

Option* 6(P) 

0 


15982.6 [0] 

15982.6 [0] 

15982.6 [0] 

1 


12210.8 [0] 

12210.8 [0] 

11906.4 [0] 

2 


10922.1 [2.7] 

10995.8 [0] 

10603.6 [0] 

3 


10011.0 [2.3] 

10117.5 [0] 

9936.6 [0] 

4 


9414.1 [2.2] 

9701.1 [0] 

9126.1 [0] 

5 


8939.0 [1.9] 

9075.9 [0] 

8869.6 [.3] 

6 


8830.8 [1.2] 

8857.4 [0] 

8885.9 [0] 

7 


8887.5 [0] 

8857.4 [0] 

8885.9 [0] 

8 

9 

10 


8887.1 [0] 
8888.5 [0] 

8857.4 [0] 

8885.9 [0] 

CPU 

Tot. 

2.615 

4.410 

4.755 

Time 

Anal . 

.418 

2.384 

3.000 

( sec) 

Opt . 

.906 

.848 

.902 


*See Table 1 
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Table 26. 



*See Table 1 


Iteration History Data for Problem 6, Case A 
One Bay / Two Story Frame (cont.) 


Weight (lb) [Maximum Constraint Violation (%) ] 


Run 4 

Run 5 

Option* 1 2 CP) 

Option* 12(D) 

15982.6 (0] 

15982.6 [0] 

10573.1 [8.3] 

11812.6 [5.7] 

9075.5 [36.4]** 

10130.0 [0] 

9007.6 [0] 

9116.0 [0] 

8886.3 [0] 

8915.5 [0] 

8886.3 [0] 

8799.3 [2.5] 

8886.3 [0] 

8791.9 [1.4] 


8826.9 [1.1] 


8848.5 [.6] 


8850.0 [.6] 


8845.6 [.6] 

4.150 

16.299 

2.203 

3.479 

.957 

11.494 


**Constraint was not Retained 
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Table 27. Final Designs for Problem 6, Case A 
One Bay / Two Story Frame 



Final Design (in) 

Linking 

Member 

Size 

Ref. 21 

Run 1 

Run 2 

Run 3 

Run 4 

Run 5 

Group 

Nos. 

Var . 


Option* 1(P) 

Option* 3(P) 

Option* 6(P) 

Option* 12(P) 

Option* 12(D) 



R 

16.099 

16.076 

16.008 

15.988 

15.951 

16.005 

i 

1,8 

t 

.3513 

.3508 

.3496 

.3492 

.3505 

.3493 



R 

11.825 

11.691 

11.726 

11.779 

11.711 

11.681 

2 

2,7 

t 

.2580 

.2553 

.2560 

.2570 

.2571 

.2549 



R 

11 .427 

11.432 

11.305 

11.293 

11.238 

11.361 

3 

3,4 

t 

.2493 

.2494 

.2468 

.2465 

.2499 

.2481 



R 

15.277 

15.134 

15.166 

15.243 

15.137 

15.158 

4 

5,6 

t 

.3333 

.3302 

.3314 

.3333 

.3348 

.3308 

Weight (lb) 


8980.7 

8888.5 

8857.4 

8885.9 

8886.3 

8845.6 

Number of Analyses 

16 

10 

9 

9 

7 

11 


*See Table 1 
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Table 29. Iteration History Data for Problem 6, Case B 
One Bay / Two Story Frame 




Weight (lb) [Maximum Constraint Violation (%) ] 

Analysis 

Run 1 

Run 2 

Run 3 

No 


Option* 1 (P) 

Option* 3(P) 

Option* 4(P) 

0 


15982.6 [0] 

15982.6 [0] 

15982.6 [0] 

1 


12210.8 [0] 

12210.8 [0] 

11889.0 [0] 

2 


10930.2 [2.4] 

11055.4 [0] 

10633.0 [1.4] 

3 


10462.8 [0] 

10491.1 [0] 

10265.0 [.7] 

4 


10212.0 [0] 

10219.2 [0] 

10231.1 [0] 

5 


10201.8 [0] 

10195.1 [0] 

10224.2 [0] 

6 

7 

8 
9 

10 

11 


10199.6 [0] 

10193.9 [0] 
10197.3 [0] 

10193.8 [0] 
10188.3 [0] 
10188.2 [0] 

CPU 

Tot. 

2.148 

3.056 

2.700 

Time 

Anal . 

.616 

1.335 

.815 

(sec) 

Opt . 

.506 

.596 

.656 


*See Table 1 
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Table 29. Iteration History Data for Problem 6, Case B 
One Bay / Two Story Frame (cont.) 



Weight (lb) [Maximum Constraint Violation (%) ] 

Analysis 

Run A 

Run 5 

No 

• 

Option* 10(P) 

Option* 10(D) 

0 


15982.6 [0] 

15982.6 [0] 

1 


10400.4 [27.4]** 

11817.5 [5.7] 

2 


10350.2 [5.9] 

10398.8 [2.9] 

3 


10252.0 [4.6] 

10364.2 [ .3] 

4 


10238.4 [.9] 

10166.8 [0] 

5 


10181.7 [0] 

10214.0 [ .7] 

6 


10181.6 [0] 

10156.4 [1.1] 

7 


10181.5 [0] 

10163.6 [.l] 

8 



10232.4 [0] 

9 



10253.8 [0] 

10 



10240.8 [0] 

11 



10220.7 [0] 

CPU 

Tot. 

2.613 

7.309 

Time 

Anal . 

.696 

.996 

(sec) 

Opt . 

.794 

4.873 


*See Table 1 


**Constraint not retained 
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Table 30. Final Designs for Problem 6, Case B 
One Bay / Two Story Frame 


Group Nos. Var. 


R 


Weight (lb) 

Number of Analyses 




Final 

Design (in) 



Ref. 21 

Run 1 

Option* 1 (P) 

Run 2 

Option* 3(P) 

Run 3 

Option* 4(P) 

Run 4 

Option* 10(P) 

Run 5 

Option* 10(D) 

15.658 

15.635 

15.599 

15.606 

15.699 

15.730 

.3416 

.3411 

.3404 

.3405 

.3426 

.3435 

13.206 

13.259 

13.269 

13.277 

13.227 

13.449 

.2881 

.2894 

.2895 

.2902 

.2901 

.2937 

13.072 

13.109 

13.015 

13.031 

12.934 

13.534 

.2852 

.2871 

.2850 

.2844 

.2822 

.2952 

17.011 

17.031 

17.133 

17.074 

17.032 

16.401 

.3711 

.3717 

.3740 

.3729 

.3719 

.3582 

10166.6 

10199.6 

10197.3 

10188.2 

10181.5 

10220.7 

22 

7 

8 

9 

8 

12 


*See Table 1 
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Table 32. Definition of Problem 7 
2x5 Grillage 


Material Properties 
Young’s Modulus : 
Shear Modulus : 
Poisson’s Ratio : 
Allowable Stress : 


E = 30.0 x 10 6 PSI 

G = 11.5 x 10 6 PSI 

v = .2963 

0 = 20,000 PSI 

a 


Nodal Loading 


Load 

Node 


Loading Components ( lb. 

in-lb) 


Case 

No. 

F 

F 

F 

M 

M 

M 



X 

y 

Z 

X 

y 

z 


i 

0 . 

0 . 

0 . 

-13330 

0 . 

0 . 


2 

0 . 

0 . 

0 . 

0 . 

37040 

0 . 


3 

0 . 

0 . 

-9000 

0 . 

-27780 

0 . 


4 

0 . 

0 . 

-3333 

0 . 

0 . 

0 . 


5 

0 . 

0 . 

0 . 

0 . 

37040 

0 . 


6 

0 . 

0 . 

-9000 

0 . 

-27780 

0 . 


7 

0 . 

0 . 

-3333 

0 . 

0 . 

0 . 

i 

8 

0 . 

0 . 

0 . 

0 . 

37040 

0 . 


9 

0 . 

0 . 

-9000 

0 . 

-27780 

0 . 


10 

0 . 

0 . 

-3333 

0 . 

0 . 

0 . 


11 

0 . 

0 . 

0 . 

0 . 

37040 

0 . 


12 

0 . 

0 . 

-9000 

0 . 

-27780 

0 . 


13 

0 . 

0 . 

-3333 

0 . 

0 . 

0 . 


14 

0 . 

0 . 

0 . 

0 . 

37040 

0 . 


15 

0 . 

0 . 

-9000 

0 . 

-27780 

0 . 


16 

0 . 

0 . 

-3333 

0 . 

0 . 

0 . 


17 

0 . 

0 . 

0 . 

13330 

0 . 

0 . 
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Table 32. Definition of Problem 7 
2x5 Grillage (cont.) 


Displacement Constraints 


Load 

Case 

Node 

No. 

Direction 

Lower Bound 

Upper Bound 


4 

z 

-0.1 in 

1 0 in 

i 

7 

z 

-0.1 in 

1.0 in 



z 

-0.1 in 

1.0 in 


I nitial Design and Side Constraints 


Member 

Nos. 

Sizing 

Variable 

Initial 
Value (in) 

Lower 
Bound (in) 

Upper 
Bound (in) 


B 

12.00 

1.00 

19.00 

1-16 

C b 

.95 

.045 

1.00 


H 

15.00 

1.00 

20.00 


C h 

.80 

.05 

.95 
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Table 33. Iteration History Data for Problem 7, Case A 
2x5 Grillage 



Volume (in^) 

[Maximum Constraint Violation (%)] 

Analysis 

Run 1 

Run 2 

Run 3 

No. 


Option* 1 (P) 

Option* 4(P) 

Option* 10(P) 

0 


32,382.4 [0] 

32,382.4 [0] 

32,382.4 [0] 

1 


15,687.5 [16.2] 

15,687.5 [16.2] 

19,946.8 [0] 

2 


9,737.9 [13.6] 

9,566.4 [14.9] 

11,727.0 [15.5]** 

3 


9,759.1 [1.9] 

9,702.2 [1.9] 

9,962.9 [.5] 

4 


8,853.3 [.2] 

8,829.2 [.1] 

8,415.9 [0] 

5 


7,961.3 [0] 

7,953.8 [0] 

7,440.8 [0] 

6 


7,418.6 [0] 

7,434.6 [0] 

6,898.1 [0] 

7 


7,123.0 [0] 

7,100.9 [0] 

6,828.7 [0] 

8 


7,087.6 [0] 

6,900.7 [0] 

6,812.7 [0] 

9 


6,872.6 [.2] 

6,887.7 [0] 

6,800.1 [0] 

10 


6,861.4 [0] 

6,840.0 [0] 

6,794.8 [0] 

11 


6,851.1 [0] 

6,771.4 [.7] 

6,786.7 [0] 

12 


6,774.9 [.5] 

6,807.4 [0] 

6,776.7 [0] 

13 


6,802.2 [0] 

6,801.5 [0] 

6,769.7 [0] 

14 


6,800.3 [0] 

6,796.4 [0] 

6,762.6 [.1] 

15 


6,794.9 [0] 


6,756.4 [0] 

CPU 

Tot . 

5.094 

5.626 

4.505 

Time 

Anal . 

1.695 

1.607 

1.642 

(sec) 

Opt . 

1.955 

2.361 

1.284 


*See Table I 


**Constraint was not Retained 
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Table 33. Iteration History Data for Problem 7, Case A 
2x5 Grillage (cont.) 



3 

Volume (in ) [Maximum Constraint Violation (%)] 

Analysis 

Run 4 

Run 5 

No 

• 

Option* 10(D) 

Option* 1(PU) 

0 


32,382.4 [0] 

32,382.4 [0] 

1 


19,943.8 [0] 

10,417.2 [7.1] 

2 


11,696.8 [16.8]** 

8,903.1 [0] 

3 


9,859.8 [.2] 

7,857.5 [0] 

4 


8,406.6 [0] 

7,166.5 [0] 

5 


7,419.1 [0] 

6,989.3 1 .11 

6 


6,874.5 [0] 


7 


6,822.5 [0] 


8 


6,800.2 [0] 


9 


6,793.7 [0] 


10 


6,783.0 [0] 


11 


6,778.7 [0] 


12 


6,772.8 [0] 


13 




14 




15 




CPU 

■SB 

3.304 


Time 


1.300 

1.105 

(sec) 

la 

.738 



*See Table 1 **Constraint was nett: Retained 
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Table 34. Final Designs for Problem 7, Case A 
2x5 Grillage 



Final Design (in) 

Linking 

Member 

Size 

Ref. 29 

Run 1 

Run 2 

Run 3 

Run 4 

Run 5 

Group 

Nos. 

Var. 


Option* 1(P) 

Option* 4(P) 

Option* 10(P) 

Option* 10(D) 

Option* 1(PU) 



B 

6.31 

8.93 

9.03 

3.35 

3.65 

19.00 + 



t b 

.045“ 

.047 

.046 

.230 

.236 

.049 

i 

1-6 

H 

18.90 

20.00 + 

19.97 

20.00 + 

20.00 + 

20.00 + 



m 

.050“ 

.050“ 

.050“ 

.050“ 

.050“ 

.050“ 



B 

6.62 

6.49 

5.54 

1.67 

1.00“ 

17.82 



mm 

.045“ 

.045“ 

.045“ 

.133 

.056 

.045“ 

2 

7-10 

mm 










15.7 

20.00 + 

19.99 

10.27 

12.10 

19.98 



m 

.050“ 

.050“ 

.050“ 

.050“ 

.050“ 

.050“ 



B 

13.50 

18.97 

19.00 

15.36 

13.09 

19.00 + 



am 

.822 

.482 

.500 

.558 

.639 

.381 

3 

13-16 

■1 









H 

20.00 + 

20.00 + 

20.00 + 

20.00 + 

20.00 + 

20.00 + 



B 

.050“ 

.050“ 

.050“ 

.050“ 

.050“ 

.050“ 

















Table 34. Final Designs for Problem 7, Case A 
2x5 Grillage (cont.) 





Final 

Design (in) 



Linking 

Member 

Size 

Ref. 29 

Run 1 

Run 2 

Run 3 

Run 4 

Run 5 

Group 

Nos. 

Yar. 


Option* 1 (P) 

Option* 4(P) 

Option* 10(P) 

Option* 10(D) 

Option* 1(PU) 



B 

4.89 

18.98 

19.00 + 

13.63 

13.45 

12.95 



t b 

.993 

.360 

.328 

.601 

.640 

.724 

4 

11-12 










H 

20.00 + 

20.00 + 

20.00 + 

20.00 + 

20.00 + 

20.00 + 



c h 

.050" 

.050' 

.050" 

.050' 

.050" 

.050' 

Volume (in^) 


6971.0 

6794.9 

6796.5 

6756.4 

6772.8 

6989.3 

Number of Analyses 

38 

16 

15 

16 

13 

6 


*See Table 1 

+Sizing Variable at Upper Bound 


-Sizing Variable at Lower Bound 


Table 35. Critical Constraints for Problem 7, Case A 
2x5 Grillage 


Run 

No. 

Option 

No.* 

Displacement 
Constrained Nodes 

Ref. 29 

— 

7,10 ' 

1 

1(P) 

7,10 

2 

A (P) 

7,10 

3 

10(P) 

4.7,10 

A 

10(D) 

4,7,10 

5 

1(PU) 

7,10 


*See Table 1 
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Table 36. Iteration History Data for Problem 7, Case B 
2x5 Grillage 





3 

Volume ( in ) [Maximum Constraint Violation (%) ] 


Analysis 

Run 1 

Run 2 

Run 3 

Run 4 

No. 


Option* 3(P) 

Option* 6(P) 

Option* 12(P) 

Option* 3(PU) 

0 


32,382.4 [0] 

32,382.4 [0] 

32,382.4 [0] 

32,382.4 [0] 

l 


25,259.7 [0] 

25,268.5 [0] 

21,663.8 [0] 

20,339.4 [01 

2 


19,997.0 [0] 

20,724.0 [0] 

13,969.4 [0] 

14,469.9 [0] 

3 


17,270.7 [0] 

17,578.1 [0] 

11,573.5 [0] 

11,099.0 [0] 

4 


14,898.8 [0] 

14,605.3 [0] 

9,712.1 [0] 

8,762.0 [0] 

5 


12,535.6 [0] 

12,916.8 [.6] 

8,644.9 [0] 

7,544.4 [100]** 

6 


11,092.8 [.2] 

11,463.6 [.7] 

7,807.2 [100]** 

7,523.9 [100] 

7 


10,032.6 [.4] 

10,630.3 [0] 

7,730.0 [20.6] 

7,526.1 [42.3]** 

8 


9,156.0 [.4] 

9,115.1 [1.8] 

7,717.4 [1.4] 

7,525.9 [4.3] 

9 


8,414.4 [79.1]** 

8,753.1 [77.4]** 

7,641.5 [0] 

7,510.4 [0] 

10 

11 

12 

13 

14 

15 

16 

17 

18 


7,881.9 [65.5]** 

7.886.5 [7.7] 
7,884.7 [9.5] 

7.887.6 [0] 

7,959.9 [51.9]** 
7,925.6 [22.2]** 

7.914.4 [9.3] 

7.917.5 [0] 

7.913.5 [0] 

7.581.5 [0] 
7,557.1 [1.0] 
7,550.4 [0] 

7.545.6 [0] 

7,505.2 [0] 


mm 

9.032 

10.459 

10.899 

14.065 

■ 


3.469 

3.945 

4.979 

4.959 


EBB 

2.486 

3.263 

2.981 

5.269 


*See Table 1 **Constraint was not retained 
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Table 37. Final Designs for Problem 7, Case B 
2x5 Grillage 



Final Design (in) 

Linking 

Member 

S iz ing 

Ref. 29 

Run 1 

Run 2 

Run 3 

Run A 

Group 

Nos. 

Variables 


Option* 3(P) 

Option* 6(P) 

Option* 12(P) 

Option* 3(PU) 



B 

6.10 

9.86 

10.60 

3.26 

4.27 



t b 

.159 

.096 

.105 

.224 

.234 

i 

1-6 

H 

20.00 + 

20.00* 

19.99 

20.00* 

20.00* 




.093 

.091 

.090 

.096 

.094 



B 

8.28 

11.29 

11.59 

'1.91 

1.84 



t b 

.074 

.093 

.095 

.148 

.131 

2 

7-10 









H 

15.20 

19.99 

19.98 

9.07 

9.14 



C h 

.064 

.074 

.074 

.055 

.056 



B 

6.33 

15.94 

16.85 

12.01 

11.64 



c b 

1.00 + 

.325 

.299 

.595 

.514 

3 

13-16 









H 

20.00* 

20.00* 

20.00* 

20.00* 

20.00* 




.098 

.095 

.094 

.103 

.098 
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Table 37. Final Designs for Problem 7, Case B 
2x5 Grillage (cont.) 



Final Design (in) 

Linking 

Member 

S iz ing 

Ref, 29 

Run 1 

Run 2 

Run 3 

Run A 

Group 

. Nos. 

Variables 


Option* 3(P) 

Option* 6(P) 

Option* 12(P) 

Option* 3(PU) 



B 

11.50 

17.62 

18.81 

14.20 

17.65 



% 

1.00 + 

.657 

.614 

.683 

.643 

4 

11-12 

H 

20.00 + 

19.99 

19.99 

20.00 + 

20.00 + 



C h 

.117 

.119 

.119 

.113 

,118 

Volume (in^) 


7927.0 

7887.6 

7913.5 

7545.6 

7505.3 

Number of Analyses 

41 

14 

15 

14 

11 


*See Table 1 

+Sizing Variable at Upper Bound 
-Sizing Variable at Lower Bound 
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Table 38. Critical Constraints for Problem 7, Case B 
2x5 Grillage 


Run 

No. 

Option 

„ * 
No. 

Displacement 
Constrained Nodes 

Stress Constrained 
Members 

Local Buckling 
Constrained Members 

Ref. 29 

— 

7,10 

— 

1,6,8,9,11,13,16 

1 

i(p) 

7,10 

— 

1,2,5,6,8,10,11,13,15 

2 

6(P) 

7,10 

— 

1,6,8,10,11,13,15 

3 

1 2 (P) 

4,7,10 

— 

1,6,8,10,11,13,15 

4 

3(PU) 

4,7,10 

8,10 

1,6,8,10,11,13,15 


*See Table 1 




Table 39, Definition of Problem 8 

Two Bay / Six Story Frame 


Material Properties 
Young’s Modulus : 
Shear Modulus : 
Poisson’s Ratio : 
Weight Density : 
Yield Stress : 
Factor of Safety : 


E = 30.0 x 10 6 PSI 
G = 11 .5 x 10 6 PSI 
V = .3 

p = .2836 lb/in 3 

O = 36,000 PSI 
a 

FS = 1.51 


Nodal Loading 


Load 

Node 


Loading Components (lb, 

in-lb) 


Case 

No. 

F 



F 

M 

■9H 

M 



X 



2 

X 

mm 

z 


i 

9000. 


0 . 

0 . 

0 . 

0 . 

0 . 


4 

9000. 


0 . 

0 . 

0 . 

0 . 

0 . 


7 

9000. 


0 . 

0 . 

0 . 

0 . 

0 . 

2 

10 

9000. 


0 . 

0 . 

0 . 

0 . 

0 . 


13 

9000. 


0 . 

0 . 

0 . 

0 . 

0 . 


16 

9000. 


0 . 

0 . 

0 . 

0 . 

0 . 
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Table 39. Definition of Problem 8 

Two Bay / Six Story Frame (cont.) 
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Table 39. Definition of Problem 8 

Two Bay / Six Story Frame (cont.) 


Displacement Constraints 


Load 

Case 

Node 
Nos . 

Direction 

Lower Bound 

Upper Bound 


H 

X 

-1.728 in 

1.728 in 


B 

X 

-1.440 in 

1.440 in 


B 

X 

-1.152 in 

1.152 in 

2 

B 

X 

- .864 in 

.864 in 


B 

X 

- .576 in 

.576 in 


B 

X 

- .288 in 

.288 in 


Initial Design and Side Constraints 


Member 
Nos . 

Sizing 

Variables 

Initial 
Value (in) 

Lower 
Bound (in) 

Upper 
Bound (in) 


R 

9.00 

1.00 

25.0 

1-10 






t 

.20 

.01 

5.0 


R 

11.00 

1.00 

25.0 

11-20 






t 

.24 

.01 

5.0 


R 

14.00 

1.00 

25.00 

21-30 






t 

.31 

.01 

5.0 
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Table 40. Iteration History Data for Problem 8, Case A 
Two Bay / Six Story Frame 



Weight (lb) [Maximum Constraint Violation (%) ] 

Analysis 

Run 1 

Run 2 

Run 3 

No 

• 

Option* 1 (P) 

Option* 3 (P) 

Option* 6(P) 

0 


23,141.4 [63.6] 

23,141.4 [63.6] 

23,141.4 [63.6] 

1 


23,503.0 [19.8] 

23,475.2 [14.7] 

23,792.0 [8.8] 

2 


23,188.3 [3.5] 

23,243.1 [1.0] 

23,793.4 [0] 

3 


23,240.5 [.5] 

23,257.5 [.1] 

23,536.4 [0] 

4 


22,945.8 [4.6] 

23,083.8 [1.1] 

23,117.4 [0] 

5 


22,804.6 [3.5] 

22,793.7 [.2] 

22,915.5 [0] 

6 


22,854.3 [.3] 

22,777.9 [0] 

22,906.4 [0] 

7 


22,668.5 [5.3] 

22,804.0 [.8] 

22,887.8 [0] 

8 


22,621.7 [2.9] 

22,608.7 [0] 


9 


22,691.1 [.7] 

22,603.1 [.2] 


10 


22,559.1 [4.7] 

22,579.8 [0] 


11 


22,645.2 [.7] 

22,566.8 [.1] 


12 


22,531.1 [2.0] 



13 


22,557.0 [.9] 



14 


22,545.7 [.5] 



15 


22,568.6 [.6] 



CPU 

BN® 

19.281 

58.239 

40.108 

Time 

E 

2.050 

43.969 

29.103 

(sec) 


12.201 

9.155 

8.043 


*See Table 1 
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Table 40. Iteration History Data f< 
Two Bay / Six Story Fraim 



Weight (lb) [Maximum 

Analysis 

Run 4 

No. 


Option* 12(P) 

0 


23,141.4 [63.6] 

1 


23,919.1 [8.7] 

2 


23,157.4 [.4] 

3 


22,999.2 [.1] 

4 


22,936.1 [0] 

5 


22,866.8 [0] 

6 


22,841.1 [0] 

7 


22,794.7 [.1] 

8 


22,747.5 [.1] 

9 


22,713.3 [.1] 

10 


22,693.4 f .1] 

11 


22,674.2 [.1] 

12 



13 



14 



CPU 

m 

56.484 

Time 


42.334 

(sec) 

EH 

10.120 


*See Table 1 


Constraint Violation (%) ] 


Run 5 

Option* 3(PU) 

23,141 .4 

[63.6] 

23,270.9 

[.8] 

23,014.6 

[.7] 

23,027.4 

[.3] 

22,982.3 

[.6] 

22,754.3 

[2.4] 

22,641.6 

[0] 

22,655.2 

[.2] 

22,534.2 

[.1] 

22,523.0 

[.1] 


54.005 

35.074 

14.229 
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Table 41. Final Designs for Problem 8 , Case A 
Two Bay / Six Story Frame 



Final Design (in) 

Linking 

Member 

Size 

Ref. 21 

Run 1 

Run 2 

Run 3 

Run 4 

Run 5 

Group 

Nos. 

Var. 


Option* 1 (P) 

Option* 3(P) 

Option* 6 (P) 

Option* 1 2 (P) 

Option* 3(PU) 

1 

1 ? 

R 

10.235 

10.250 

10.192 

9.985 

10.175 

10.215 



t 

.2233 

.2238 

.2225 

.2181 

.2224 

.2229 

2 

3,5 

R 

10.571 

10.462 

10.396 

8.732 

10.351 

10.422 


t 

.2306 

.2285 

.2276 

.3030 

.2286 

.2272 

3 

4 

R 

5.547 

7.063 

7.669 

8.565 

8.193 

7.397 



t 

.1210 

.1544 

.1681 

.1869 

.1798 

.1608 

4 

6 7 

R 

9.808 

9.813 

9.674 

9.712 

9.742 

9.688 


o, / 

t 

.2140 

.2143 

.2116 

.2122 

.2132 

.2115 


8,10 

R 

7.208 

7.389 

6.802 

8.015 

7.383 

6.786 

J 

t 

.1573 

.1616 

.1496 

.1751 

.1627 

.1493 

5 

Q 

R 

9.140 

8.837 

9.047 

8.648 

8.716 

9.096 


j 

t 

.1994 

.1937 

.1978 

.1890 

.1905 

.1992 

7 

11,12 

R 

10.082 

10.139 

10.092 

10.019 

10.096 

10.104 

t 

.2200 

.2212 

.2203 

.2192 

.2207 

.2208 
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Table 41, Final Designs for Problem 8, Case A 
Two Bay / Six Story Frame (cont.) 



Final Design (in) 

Linking 

Member 

Size 

Ref. 21 

Run 1 

Run 2 

Run 3 

Run 4 

Run 5 

Group 

Nos. 

Var. 


Option* 1(P) 

Option* 3(P) 

Option* 6(P) 

Option* 12 (p) 

Option* 3 (PU) 

Q 

13 IS 

R 

10.535 

10.356 

10.491 

8.864 

10.115 

10.524 

o 

i j j x j 

t 

.2298 

.2262 

.2293 

.2628 

.2332 

.2291 

Q 

14 

R 

9.233 

9.212 

9.541 

9.859 

9.688 

9.199 

-7 


t 

.2014 

.2013 

.2089 

.2151 

.2118 

.2004 

10 

16 17 

R 

10.395 

10.420 

10.432 

10.390 

10.422 

10.478 



t 

.2268 

.2273 

.2278 

.2269 

.2279 

.2293 

1 1 

18 20 

R 

8.973 

8.963 

8.814 

9.238 

8.943 

8.661 

X X 

1 vl | XU 

t 

.1958 

.1955 

.1933 

.2022 

.1963 

.1897 

12 

19 

R 

11.336 

11.150 

10.835 

9.918 

10.899 

11.355 


X J 

t 

.2473 

.2444 

.2596 

.2713 

.2490 

.2486 

13 

21 22 

R 

11.076 

11.176 

11.182 

11.182 

11.174 

11.182 



t 

.2417 

.2439 

.2441 

.2440 

.2442 

.2440 

14 

23 25 

R 

10.619 

10.299 

10.349 

10.203 

10.229 

10.444 


4 J j L J 

t 

.2317 

.2260 

.2258 

.2232 

.2525 

.2273 






Table 41. Final Designs for Problem 8, Case A 
Two Bay / Six Story Frame (cont.) 



Final Design (in) 

Linking 

Member 

Size 

Ref . 21 

Run 1 

Run 2 

Run 3 

Run 4 

Run 5 

Group 

Nos. 

Var. 


Option* 1 (P) 

Option* 3(P) 

Option* 6(P) 

Option* 12(P) 

Option* 3 (PU) 

15 

24 

R 

11.236 

10.319 

11.274 

11.588 

11.345 

11.187 



t 

.2452 

.2476 

.2468 

.2539 

.2490 

.2448 

16 

26,27 

R 

10.986 

11.014 

11.094 

10.933 

11.016 

11.116 



t 

.2397 

.2408 

.2422 

.2389 

.2408 

.2430 

17 

28,30 

R 

10.687 

11.095 

10.989 

11.810 

11.258 

10.790 



t 

.2332 

.2421 

.2412 

.2581 

.2467 

.2365 

18 

29 

R 

15.021 

14.141 

13.168 

12.193 

13.572 

13.352 



t 

.3277 

.3126 

.3455 

.3249 

.3221 

.3 508 

Weight (lb) 


22530.5 

22568.6 

22566.8 

22887.8 

22674.2 

22523.0 

Number of Analyses 

23 

16 

12 

8 

12 

10 


*See Table 1 
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Table 42. Critical 
Two Bay / 



*See Table 1 


nstraints for Problem 8, Case A 
ix Story Frame 


Buckling Constrained 
Members 

Load Case 1 

Load Case 2 

1,3,4,11, 

13,14,20,24 

9,16,19,21, 

22,24-26,29 

1,3,11,13, 

14,20,23 

9,16,19-22, 

24-26,29,30 

1,3,10,11, 

13,20,23,24 

9,10,16,19-22 

24-26,29,30 

1,3,11,13, 

20,23 

16,17,19-22, 

24-27,29,30 

1,3,11,13, 

20,23,24 

9,16,19-22, 

24-26,29-30 

1,3,10,11,13, 

14,20,23,24 

9,16,19-22, 

24-26,29-30 


R/ t Constrained 
Members 


1,2,4,6-12,14, 

16-18,20-28,30 


1-12,14,16-28,30 
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Table 43, Iteration History Data for Problem 8, Case B 
Two Bay / Six Story Frame 


Analysis 

No. 


CPU 

Time 

(sec) 



Weight (lb) [ 

Maximum Constraint \ 

Violation (%)] 


Run 

1 

Run 

2 

Run 

3 

Option* 

: 1(P) 

Option* 

3 (P) 

Option* 

4 (P) 

28 

,536.8 

[55.7] 

28,536.8 

[55.7] 

28,536.8 

[55.7] 

27 

,484.3 

[16.9] 

27,184.2 

[9.9] 

27,270.2 

[10.0] 

25 

,692.2 

[4.3] 

25,495.6 

[.7] 

25,546.8 

[1-8] 

25 

,404.1 

[.8] 

25,168.4 

[0] 

25,179.7 

[0] 

25 

,356.5 

[.2] 

25,056.6 

[1.5] 

25,128.6 

1.5] 

25 

,268.1 

[0] 

24,651.6 

[0] 

24,871.5 

[.2] 

24 

,939.8 

[.2] 

24,689.4 

[0] 

24,889.6 

[0] 

24 

,853.9 

[-4] 

24,488.8 

[0] 

24,818.5 

[0] 

24 

,865.1 

[1.0] 

24,530.4 

[0] 

24,855.5 

[0] 

24 

,672.7 

[0] 

24,339.1 

[0] 

24,808.1 

[0] 

24 

,661.9 

[.7] 

24,363.1 

[0] 

24,644.2 

[0] 

24 

,666.6 

[.4] 

24,314.3 

[0] 

24,598.2 

[0] 




24,278.8 

[0] 

24,618.5 

[0] 






24,608.8 

[0] 

17 

.205 


54.538 


23.801 


7 

.316 


39.959 


8.639 


5 

.702 


10.048 


10.435 



*See Table 1 















Table 43, Iteration History Data for Problem 8, Case b 
Two Bay / Six Story Frame (cont.) 



Weight (lb) [Maximum Constraint Violation (%)] 

Analysis 

Run 4 

Run 5 

No. 


Option* 10(P) 

Option* 1(PU) 

0 


28,536.8 [55.7] 

28,536.8 [55.7] 

1 


25,656.7 [19.1]** 

26,576.6 [10.0] 

2 


24,952.5 [2.5] 

24,952.4 [1.7] 

3 


24,654.2 [.6] 

24,889.6 [1.0] 

4 


24,609.4 [.2] 

24,890.1 [1.8] 

5 


24,624.1 [0] 

24,676.1 [1.4] 

6 


24,615.2 [0] 

24,b75.9 [.9] 

7 



24,665.6 [.3] 

8 




9 




10 




11 




12 




13 




14 




15 




CPU 

Tot . 

13.617 

19.436 

Time 

Anal . 

3.749 

6.483 

(sec) 

Opt. 

7.324 

8.018 


*See Table 1 
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Table 44. Final Designs for Problem 8, Case B 
Two Bay / Six Story Frame 



Final Design (in) 

Linking 

Member 

Size 

Ref. 21 

Run 1 

Run 2 

Run 3 

Run A 

Run 5 

Group 

Nos. 

Var. 


Option* 1(P) 

Option* 3(P) 

Option* 4(P) 

Option* 10(P) 

Option* 1(PU) 

i 

1 9 

R 

10.009 

10.071 

9.913 

10.110 

10.008 

10.079 

l 


t 

.2183 

.2203 

.2165 

.2210 

.2226 

.2201 

o 


R 

10.361 

10.298 

10.137 

10.342 

10.240 

10.296 

4 

J , 0 

t 

.2261 

.2252 

.2214 

.2266 

.2273 

.2250 


A 

R 

6.322 

7.871 

6.494 

6.628 

6.353 

7.062 

J 


t 

.1379 

.1732 

.1432 

.1446 

.1504 

.1709 

A 

A 7 

R 

9.682 

9.771 

9.521 

9.690 

9.647 

9.780 


o, / 

t 

.2112 

.2141 

.2078 

.2120 

.2148 

.2137 


8,10 

R 

7.825 

8.338 

7.746 

7.681 

7.655 

7.923 

_/ 

t 

.1707 

.1826 

.1705 

.1676 

.1706 

.1823 

A 

q 

R 

10.953 

9.748 

11.275 

11.179 

10.796 

10.345 

u 


t 

.2390 

.2128 

.2462 

.2439 

.2407 

.2259 

7 

11,12 

R 

10.639 

10.637 

10.813 

10.827 

10.719 

10.652 

/ 

t 

.2321 

.2322 

.2363 

.2363 

.2374 

.2324 


) 
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Table 44. Final Designs for Problem 8, Case B 
Two Bay / Six Story Frame (cont.) 





Final 

Design (in) 



Linking 

Member 

Size 

Ref. 21 

Run 1 

Run 2 

Run 3 

Run 4 

Run 5 

Group 

Nos • 

Var. 


Option* 1 (P) 

Option* 3(P) 

Option* 4(P) 

Option* 10(P) 

Option* 1 (PU) 

8 

13,15 

R 

9.891 

9.668 


9.781 

9.615 


t 

.2158 

.2127 


.2135 

.2132 



14 

R 

10.763 

11.211 

11.214 

10.532 

10.777 

11.330 

9 

t 

.2348 

.2448 

.2469 

.2298 

.2399 

.2474 

10 

16,17 

R 

11.406 

11.687 

11.585 

11.665 

11.642 

11.675 

t 

.2488 

.2557 

.2535 

.2556 

.2578 

.2554 

11 

18,20 

R 

8.687 

9.675 

8.782 

9.148 

9.188 

9.440 

t 

.1895 

.2116 

.1947 

.2004 

.2015 

.2106 

12 

19 

R 

13.104 

11.744 

12.826 

12.618 

12.377 

12.208 

t 

.2859 

.2565 

.2801 

.2754 

.2737 

.2666 

13 

21,22 

R 

12.114 

12.235 

12.108 

12.217 

12.183 

12.223 

t 

.2643 

.2678 

.2644 

.2674 

.2694 

.2668 

14 

23,25 

R 

10.690 

10.417 

10.456 

10.392 

10.309 

10.117 

t 

.2332 

.2287 

.2282 

.2267 

.2268 

.2271 














Table 44. Final Designs for Problem 8, Case B 
Two Bay / Six Story Frame (cont.)* 



Final Design (in) 

Linking 

Member 

Size 

Ref. 21 

Run 1 

Run 2 

Run 3 

Run 4 

Run 5 

Group 

Nos. 

Var. 


Option* 1(P) 

Option* 3(P) 

Option* 4(P) 

Option* 10(P) 

Option* 1(PU) 

15 

24 

R 

11.919 

12.441 

12.201 

12.018 

12.092 

12.414 

t 

.2601 

.2715 

.2696 

.2662 

.2694 

.2706 

16 

26,27 

R 

11 .707 

11.577 

11.225 

11.584 

11.576 

11.729 

t 

.2554 

.2533 

.2519 

.2531 

.2552 

.2559 

17 

28,30 

R 

11.413 

11.636 

9.480 

11.493 

11.385 

11.468 

t 

.2490 

.2545 

.2068 

.2513 

.2520 

.2533 

18 

29 

! 

R 

13.915 

13.067 

16.746 

13.801 

13.557 

13.155 

t 

.3036 

.2851 

.3654 

.3013 

.2987 

.2866 

Weight (lb) 


24405.4 

24666.6 

24278.8 

24608.8 

24615.2 

24665.6 

Number of Analyses 

27 

12 

13 

14 

7 

8 


*See Table 1 
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Table 


Run 

No. 

Option 

XT * 

No . 

Displacement 

Constrained 

Nodes 



Load Case 1 

Ref. 21 

— 

1-12 

1 

i (p> 

4-12 

2 

3 (P) 

1-12 

3 

4 (P) 

4-12 

4 

10(P) 

| 

1-12 

5 

1(PU) 

4-12 


*See Table 1 


Critical Constraints for Problem 8, Case B 
Two Bay / Six Story Frame 


Stress Constrained Buckling Constrained R/t 

Members Members Constrained 

Members 


Load Case 1 I Load Case 2 I Load Case 1 I Load Case 2 


1,3,4,12,20 25,29 1-30 
1,3,13 29,30 1-30 
1,3,4,13,30 25 1-30 
1,3,13 25,29 1-30 

25,29 1-3,5-30 


1,3,13 

1,3,13 


25,29,30 


1-3, 5-8 
10-30 




















Table 46. Definition of Problem 9 
Helicopter Tail Boom 

Material Properties 


Young's Modulus 

: E = 

10.5 x 10 6 PSI 

Shear Modulus 

: G = 

40.4 x 10 5 PSI 

Poisson's Ratio 

: v = 

.3 

Weight Density 

: p = 

.1 lb/in 3 


4 

Allowable Stress : a = 4.2 x 10 PSI 

a 


Factor of Safety : FS = 1.25 

Nodal Loading 


Load 

Case 

1 

Loading Components (lb, in-! 

lb) 


F 

X 

F 

y 

F 

z 

M 

X 

M 

y 

M 

z 


13-16 

0. 

0. 

-140.0 

0. 

0. 

0. 


25 

1490.3 

1691.8 

0. 

0. 

0. 

0. 

i 

26 

1490.3 

-1365.8 

0. 

0. 

0. 

0. 


27 

-1490.3 

1691.8 

0. 

0. 

0. 

0. 


28 

-1490.3 

-1365.8 

0. 

0. 

0. 

0. 


Displacement Constraints 


Load 

Case 

■ 

■BlSH 

Direction 

Lower Bound 

Upper Bound 


5-28 

y 

-.5 in 

. 5 in 

i 

5-28 

z 

-.5 in 

. 5 in 


Initial Design and Side Constraints 


Member 

Nos. 

Sizing 

Variable 

Initial 
Value (in) 

Lower 
Bound (in) 

Upper 
Bound (in) 


R 

2.0 

.25 

25.0 

1-44 




5.0 

t 

.051 

.001 


396 

































397 


) 


f 


Table 47. Iteration History Data for Problem 9 
Helicopter Tail Boom 


Weight (lb) [Maximum Constraint Violation (%)] 


Analysis 

No. 



Run 1 

Run 2 

Run 3 

Option* 1(P) 

Option* 4(P) 

Option* 10(P) 

69.11 [211.6] 

69.11 [211.6] 

69.11 [211.6] 

95.86 [34.7] 

93.97 [34.0] 

99.51 [18.6] 

107.94 [7.2] 

105.77 [7.2] 

109.45 11. 0] 

110.69 [.4] 

109.53 [.2] 

109.59 [.1] 

109.51 [.1] 

109.26 [0] 

109.52 [0] 

110.13 [0] 

109.59 [0] 

109.45 [0] 

109.33 [0] 

109.33 [0] 

109.40 [0] 

109.25 [0] 

109.12 [0] 

109.36 [0] 

108.83 [0] 

109.18 [0] 

109.27 [0] 

108.71 [0] 

108.60 [0] 

109.22 [0] 

108.80 [01 

108.66 [0] 


15.819 

16.337 

13.664 

9.368 

9.904 

8.426 

2.661 

2.652 

1.991 


*See Table 1 




Table 47. Iteration History Data for Problem 9 
Helicopter Tail Boom (cont.) 



Weight (lb) [Maximum Constraint Violation (%)] 

Analysis 

No. 

Run 4 

Option* 10(D) 

Run 5 

Option* 1(PU) 

0 


69.11 

[211.6] 

69.11 

[211.6] 

1 


97.87 

[20.3] 

105.44 

[7.3] 

2 


108.65 

[1.7] 

112.34 

[0] 

3 


109.62 

[0] 

110.20 

[0] 

4 


108.74 

[-4] 

109.15 

[0] 

5 


108.34 

[-5] 

108.90 

[0] 

6 


108.74 

[0] 

108.70 

[0] 

7 


108.66 

[0] 



8 


108.48 

[0] 



9 


108.52 

[0] 



10 


108.35 

[0] 



CPU 

Time 

(sec) 

Tot. 
Anal . 
Opt ♦ 

41.147 

9.367 

28.039 

13.498 

6.379 

3.235 


*See Table 1 
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Table 48. Final Designs for Problem 9 
Helicopter Tail Boom 



Final Design (in) 

Linking 

Member 

Size 

Ref. 21 

Run 1 

Run 2 

Run 3 

Run 4 

Run 5 

Group 

Nos. 

Var . 


Option* 1(P) 

Option* 4(P) 

Option* 10(P) 

Option* 10(D) 

Option* 1(PU) 

1 

1-4 

R 

2.6695 

3.1432 

3.0791 

3.1473 

3.0816 

3.1198 



t 

.0880 

.0801 

.0792 

.0811 

.0782 

.0794 

2 

5-8 

R 

1.9152 

1.2850 

1.3675 

1.4169 

1.3848 

1.3364 


t 

.0487 

.0355 

.0377 

.0380 

0.372 

.0386 

3 

9 -12 

R 

2.6530 

2.8813 

2.9242 

2.8725 

2.8850 

2.9110 


t 

.0829 

.0744 

.0744 

.0738 

.0736 

.0737 

4 

13-16 

R 

2.1035 

2.0071 

1.9918 

1.9786 

1.9927 

1.9975 


t 

.0535 

.0511 

.0509 

.0512 

.0507 

.0509 

5 

17-20 

R 

2.6784 

2.8101 

2.8255 

2.8328 

2.8086 

2.8239 


t 

.0753 

.0717 

.0719 

.0736 

.0716 

.0717 

6 

21-24 

R 

2.1488 

2.0703 

2.0709 

2.0573 

2.0656 

2.0756 


t 

.0547 

.0527 

.0527 

.0533 

.0526 

.0529 

7 

25-28 

R 

2.6238 

2.6724 

2.6681 

2.6071 

2.6848 

2.6672 


t 

.0673 

.0680 

.0679 

.0674 

.0685 

.0678 

8 

29-32 

R 

2.1569 

2.0965 

2.0959 

2.0862 

2.0884 

2.0983 

t 

.0549 

.0535 

.0533 

.0540 

.0532 

.0534 


400 


Table 48. Final Designs for Problem 9 
Helicopter Tail Boom (cont.) 



Final Design (in) 

Unking 

Group 

Member 

Nos. 

Size 

Var. 

Ref. 21 

Run 1 

Option* 1 (P) 

Run 2 

Option *4 ( P ) 

Run 3 

Option* 10(P) 

Run 4 

Option* 10(D) 

Run 5 

Option* 1(PU) 

9 

33-36 

R 

t 

2.5038 

.0637 

2.5179 

.0642 

2.5101 

.0639 

2.4670 

.0638 

2.4965 

.0637 

2.5103 

.0640 

10 

37-40 

R 

t 

2.1730 

.0553 

2.1475 

.0548 

2.1459 

.0546 

2.1347 

.0533 

2.1583 

.0550 

2.1480 

.0547 

11 

41-44 

R 

t 

2.3748 

.0604 

2.3703 

,0605 

2.3664 

.0602 

2.3654 

.0603 

2.4109 

.0615 

2.3618 

.0602 

12 

45-48 

R 

t 

1.9707 

.0502 

1.9487 

.0497 

1.9476 

.0496 

1.9300 

.0501 

1.9938 

.0508 

1.9485 

.0496 

Weight (lb) 


111 .20 

108.80 

108.66 

109.22 

108.35 

108.70 

Number of Analyses 

13 

11 

11 

10 

11 

7 


*See Table 1 
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Table 49. Critical Constraints for Problem 9 
Helicopter Tail Boom 


Run 

No. 

Option 

No.* 

Displacement 
Constrained Nodes 

R/ 1 Constrained 
Members 

Ref. 21 

— 

25,27 

5-8,13-16,21-24,29-48 

1 

1 (P) 

25,27 

1-4,9-48 

2 

4(P) 

25,27 

1-4,9-48 

3 

10(P) 

25,27 

1-4,9-48 

4 

10(D) 

25,27 

1-4,9-48 

5 

1(PU) 

25,27 

1-4,9-48 


*See Table 1 
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